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Abstract 

To  increase  the  performance  of  the  current  EiS  satellite  launeh  eapability,  new 
roeket  designs  must  be  undertaken.  One  eoneept  that  has  been  around  since  the  50s  but 
yet  to  be  utilized  on  a  launeh  platform  is  the  aerospike,  or  plug  nozzle.  The  aerospike 
nozzle  concept  demonstrates  globally  better  performance  compared  to  a  conventional  bell 
nozzle,  sinee  the  expansion  of  the  jet  is  not  bounded  by  a  wall  and  therefore  ean  adjust  to 
the  environment  by  changing  the  outer  jet  boundary.  A  dual-expander  aerospike  nozzle 
(DEAN)  roeket  concept  would  exceed  the  Integrated  High  Payoff  Roeket  Propulsion 
Teehnology  initiative  (IHPRPT)  phase  three  goals.  This  document  covers  the  design  of 
the  chamber  and  nozzle  of  the  DEAN.  The  validation  of  the  design  of  the  DEAN  are 
based  on  the  model  in  Numerieal  Propulsion  System  Simulation  (NPSS  added  with 
the  nozzle  design  from  Two-Dimensional  Kinematies  (TDK  04™).  The  result  is  a  rocket 
engine  that  produces  57,23 1  Ibf  (254.5  kN)  of  thrust  at  an  kp  of  472  s.  Additionally,  the 
oxygen  wall  is  made  of  silieon  carbide,  with  a  melting  point  of  5580  R  (3100  K),  and  has 
a  maximum  temperature  at  the  throat  of  1625  R  (902  K).  The  hydrogen  side  is  made  of 
copper,  with  a  melting  point  of  2444  R  (1358  K),  and  has  a  maximum  wall  temperature 
of  1224  R  (680  K)  at  the  throat.  Based  on  these  result,  future  investigation  into  this 
design  is  merited  since  it  has  the  potential  to  save  $19  million  in  the  eost  to  launch  a 
satellite.  NPSS  proved  to  be  a  powerful  tool  in  the  development  of  rocket  engines.  TDK, 
however,  was  left  wanting  in  the  area  of  aerospike  design. 
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COMPUTATIONAL  DESIGN  OF  UPPERSTAGE  CHAMBER,  AEROSPIKE,  & 
COOEING  JACKET  OF  DUAL-EXPANDER  ROCKET  ENGINE 

1.  Introduction 

Engineers  live  to  design  at  the  edge  of  what  is  possible.  Sinee  the  beginning  of 
rocketry,  the  bell  nozzle  has  dominated  the  rocket  motor  design.  Now  the  time  is  right  to 
revive  an  old  concept  and  couple  it  with  new  technology  to  break  the  bell’s  dominance. 
The  concept  is  the  aerospike  nozzle  and  the  dual-expander  aerospike  nozzle,  or  DEAN,  is 
the  design  that  will  usher  in  an  age  of  improved  performance  and  cost  in  the  space  launch 
arena. 

1.1.  Motivation 

The  Integrated  High  Payoff  Rocket  Propulsion  Technology  initiative  (IHPRPT)' 
is  a  joint  government  and  industry  effort  focused  on  developing  affordable  technologies 
for  reach  capability,  sustainable  strategic  missiles,  long  life  or  increased  maneuverability 
spacecraft  capability  and  high  performance  tactical  missile  capability.  The  objectives  for 
the  boost  and  orbit  transfer  part  is  to  increase  the  specific  impulse  (fip),  increase  the 
thrust-to-weight  ratio  or  the  mass  fraction,  reduce  the  failure  rate,  and  increase  the 
reusability.  All  of  this  will  lead  to  the  end  goal  of  a  reduction  in  the  cost  of  launching 
satellites.  The  goals  of  phase  III  are  to  increase  the  fip  by  26  seconds  and  improve  thrust- 
to-weight  by  100%.  The  result  of  the  program  should  be  an  increase  in  payload  of  22% 
and  95%  and  a  reduction  in  cost  of  33%  and  82%  for  expendable  and  reusable  launch 
vehicles,  respectively. 

One  way  to  determine  the  cost  of  a  launch  is  per  pound.  Table  1  shows  the 
estimated  average  cost  per  pound. 
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Table  1.  Average  Price  per  Pound  of  Launch  Vehicles^ 


Vehicle  Class 

LEO 

GTO 

Western 

Non- 

Western 

Western 

Non- 

Western 

Small 

$8,445 

$3,208 

$18,841 

N/A 

Medium 

$4,994 

$2,404 

$12,133 

$9,843 

Heavy 

$4,440 

$1,946 

$17,032 

$6,967 

In  a  simple  analysis  detailed  in  Appendix  D,  for  the  same  ehange  in  veloeity  (An) 
a  one-second  increase  in  hp  results  in  a  savings  of  134.5  Ibm  (61  kg).  Based  on  the 
numbers  from  Table  1  for  medium  western  vehicles,  this  results  in  a  savings  of  as  much 
as  $671,693  a  launch,  or  larger  and  more  capable  satellites.  Moreover,  that  only 
represents  a  single  second  increase  in  hp  so  the  actual  savings  could  be  dramatically 
greater. 

1,2.  Objectives 

To  meet  the  goals  of  IHPRPT,  a  Dual  Expander  Aerospike  Nozzle  engine 
(DEAN)  is  being  designed  to  provide  50,000  Ibf  (222.4  kN)  of  thrust  with  an  kp  of  464  s. 
This  research  effort  focuses  on  the  development  of  a  thrust  chamber,  nozzle,  and  cooling 
jacket  in  an  attempt  to  satisfy  the  following  three  goals. 

1 .  Determine  feasibility  of  meeting  the  IHPRPT  Phase  III  orbit  transfer  vehicle 
goals  with  the  DEAN  concept 

2.  Implement  and  improve  upon  a  design  process  focused  on  the  energy 
conversion  section  of  a  rocket  engine  (combustion  chamber,  nozzle) 
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3.  Perform  detailed  design  analysis  of  the  energy  transfer  eomponents  (eooling 
jaekets)  making  the  DEAN  possible 
1.2.1.  The  DEAN 

The  DEAN  will  utilize  liquid  hydrogen  and  oxygen  as  the  fuel  and  oxidizer.  Eaeh 
will  operate  in  their  own  expander  cyele  powered  by  the  heat  from  the  thrust  ehamber 
and  nozzle.  The  DEAN  is  also  designed  with  an  aerosptke  or  plug  nozzle.  The  DEAN 
has  the  potential  to  dramatieally  reduee  the  eost  of  launehing  a  satellite.  Eigure  1  show  a 
sehematie  of  the  Dean. 


Figure  1,  The  Dean  Schematic 


Eigure  1  shows  that  eaeh  of  the  fluid  flows  is  eontain  their  own  turbo-maehinery 
and  cool  a  separate  wall  (oxygen  cools  the  outer-wall  and  hydrogen  cools  the  inner-wall). 
Eour  system  design  choices  make  the  DEAN  a  revolutionary  engine.  Eirst,  the  Dean 
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utilizes  the  aerospike  nozzle  that  will  dramatically  reduce  the  weight  and  improved 
performance  over  a  bell  nozzle.  Since  the  fuel  and  oxidizer  are  separated  until  injection 
into  the  chamber  there  is  no  need  for  inter  propellant  seals  thus  reducing  a  critical  failure 
mode  increasing  reliability.  The  split  flow  on  the  fuel  side  reduces  the  required 
horsepower  that  increases  the  life  of  the  turbo-machinery.  Lastly,  the  design  has  the 
ability  to  be  throttled. 

In  order  to  analyze  the  design  of  the  DEAN,  two  computational  programs  will  be 
used.  Two-dimensional  kinematics  (TDK  04™)  will  be  used  to  design  the  contour  of  the 
nozzle.  Numerical  Propulsion  System  Simulation  (NPSS  ™)  will  be  used  to  asses  the 
performance  of  all  the  components  of  the  DEAN.  The  thrust  chamber,  nozzle,  and 
cooling  jacket  are  the  focused  of  this  document. 

1.2.2.  Thrust  Chamber 

The  thrust  chamber  embodies  the  essence  of  rocket  propulsion:  the  acceleration 
of  matter  and  the  reaction  imparting  propulsive  force  to  the  vehicle.  The  aim  is  to 
achieve  a  device  of  maximum  performance,  stability,  durability  while  minimizing  the 
size,  weight,  and  cost.  In  this  report,  the  thrust  chamber  consists  of  the  combustion 
chamber  and  the  injector.  The  combustion  chamber  provides  a  volume  for  proper  mixing 
of  the  propellants  and  length  for  complete  combustion.  The  injector  distributes  the 
prescribed  propellant  mass  flows  to  the  chamber.  The  key  variables  for  the  thrust 
chamber  are  the  geometry,  pressure,  and  temperature  of  the  chamber. 

1.2.2.  Nozzle 

The  nozzle  is  directly  connected  to  the  combustion  chamber.  The  nozzle  converts 
the  enthalpy  of  the  hot  combusted  gases  into  kinetic  energy  and  produces  the  thrust  of  the 
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engine.  The  key  variables  affeeted  by  the  design  of  the  nozzle  are  thrust,  kp,  nozzle 
length,  and  expansion  ratio  (s).  Maximizing  the  kp  is  benefieial  to  any  roeket  design  as 
shown  earlier.  In  the  design  of  the  nozzle,  the  ideal  length  is  often  quite  long,  adding 
weight  to  the  engine.  Therefore,  a  design  eonsideration  is  to  minimizing  the  length  of  the 
nozzle  while  still  maintaining  performanee  near  the  ideal  ease. 

1.2.4.  Cooling  Jacket 

The  eombustion  temperatures  in  the  thrust  ehamber  are  extremely  high. 
Additionally,  the  heat-transfer  rates  from  the  eombusted  gases  to  the  wall  are  high. 
Consequently,  the  eooling  jaeket  requires  major  design  eonsideration.  Not  only  does  the 
jaeket  need  to  keep  the  walls  eool  enough  to  maintain  their  struetural  integrity,  adequate 
heat  must  be  transferred  to  the  eooling  fluids  to  power  the  turbines.  The  eooling  jaeket 
will  be  analyzed  based  on  the  temperature  of  the  ehamber  wall  and  the  temperature 
ehange  in  the  eooling  fluid. 

1,3.  Preview 

The  remainder  of  this  report  begins  with  review  of  the  subjeets  relevant  to  the 
thrust  ehamber,  nozzle,  and  eooling  jaeket.  The  methodology  employed  to  eonduet  the 
work  within  this  report  is  presented  next.  The  results  obtained  from  the  methodology 
deseribed  is  presented  and  diseussed.  Finally,  the  eonelusions  inferred  based  on  the 
results  of  this  report  are  stated  and  any  reeommendations  based  on  the  work  done  are 
deelared. 
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2,  Background 


A  substantial  amount  of  work  has  already  been  completed  covering  the  wide 
range  of  topics  relating  to  the  design  in  this  thesis.  The  purpose  of  this  Chapter  is  to 
present  the  prior  knowledge  that  exists  which  can  set  benchmarks  and  help  in  the 
execution  of  the  design  of  the  DEAN.  Ideas  and  benchmarks  from  these  works  were 
considered  in  the  design  of  the  three  main  elements  in  this  thesis. 

2,1,  Orbit  Transfer  Engines  and  Mission  Requirements 

2.1.1.  Baseline  Engine 

The  DEAN  was  primarily  designed  to  replace  the  Pratt  and  Whitney  RE  10  rocket 
engine.  The  original  RE  10  was  designed  in  1959  as  an  upper-stage  liquid-oxygen  liquid- 
hydrogen  expander  cycle  rocket  engine.  The  most  current  inceptions  of  the  RE  10  is  the 
RLlOB-2.  The  RElOB-2  features  the  world’s  largest  carbon-carbon  extendible  nozzle. 
This  high-expansion  ratio  nozzle  enables  the  RE  1  OB-2  to  achieve  a  remarkable  465.5  sec 
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of  specific  impulse  and  lift  payloads  of  up  to  30,000  Ibm.  Figure  2  shows  an  image  of 
theRLlOB-2. 
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Figure  2.  RL10B-2(with  permission  from  Pratt  &  Whitney)^™^'  Bookmark  not  defined 

In  Figure  2,  the  RLlOB-2  is  shown  inside  the  expandable  skirt.  The  RLlOB-2 
eurrently  powers  the  upper  stage  of  the  medium  and  heavy-lift  versions  of  Boeing’s  Delta 
IV  for  Evolved  Expendable  Launeh  Vehiele  (EELV),  government  and  eommereial 
missions. 

The  RE  10  is  also  used  on  the  Centaur  upper  stage.  With  the  Titan  IVB  and  the 
Centaur  upper-stage,  they  are  able  to  insert  payloads  greater  than  12,700  pounds  direetly 
into  geosynehronous  orbit.  The  high-energy  Centaur  upper  stage  has  evolved  to  beeome  a 
very  versatile  vehiele.  Performing  a  three-bum  mission,  the  Centaur  aehieves  parking 
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orbit  with  the  first  bum,  boosts  itself  and  the  satellite  to  a  highly  elliptieal  orbit  with 
seeond  burn,  and  cireularizes  the  orbit  at  geosynchronous  altitude  with  the  third  burn"^. 

The  Centaur  propulsion  system  uses  2  RL10A-3-3A  Pratt  &  Whitney  engines. 
Each  engine  produces  16,500  Ibf  of  thmst,  a  444.4  sec  nominal  fip  at  5.0:1  mixture  ratio 
and  an  area  ratio  of  61:1."^  This  series  of  engines  has  been  used  successfully  since  1963 
on  the  Saturn  and  Atlas/Centaur  vehicles.  The  RL10A-3-3A  uses  an  expander  cycle, 
where  all  of  the  LH2  is  burned  in  the  combustion  chamber,  except  for  a  small  amount 
used  for  autogenously  pressurization  and  pump  bearing  cooling/gear  box  pressurization."^ 
The  turbine  working  fluid  is  the  supercritical  hydrogen  heated  in  the  regeneratively 
cooled  thmst  chamber. 

2.1.2.  Typical  Mission  Requirements 

One  of  the  simplest  and  most  common  methods  of  putting  a  spacecraft  into 
geostationary  orbit  involves  three  steps:  launch  to  low  earth  (or  parking)  orbit  with 
chemical  propulsion;  erect  a  GTO  (geostationary  transfer  orbit)  with  an  additional 
chemical  stage;  and  perform  a  simultaneous  circularization-and-plane-change  maneuver 
at  the  apogee  of  the  GTO.^ 

With  an  Atlas/Centaur  class  launch  vehicle,  the  conventional  path  to  the 
geostationary  orbit  (GEO)  is  a  two-burn  Centaur  stage.  The  first  burn  establishes  a 
slightly  elliptical  parking  orbit  and  the  second  bum  places  the  satellite  in  a  geostationary 
transfer  orbit.  ^ 
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2.1.3.  Expander  Cycles 

In  general  there  are  three  elassic  engine  system  eonfigurations;  gas-generator 
eyele,  expander  cyele,  and  staged  eombustion  eycle.  The  DEAN  utilizes  the  expander 
eycle.  The  expander  cyele  places  the  turbine  inline  with  the  thrust  chamber,  exhausting 
directly  into  the  chamber.^  Figure  3  shows  a  schematic  of  the  expander  cycles  of  the 
RLlOB-2. 


Figure  3,  RLlOB-2  Cycle  Schematic  (with  permission  from  Pratt  &  Whitney)^"^”' 

Bookmark  not  defined. 


Like  in  Figure  3,  most  expander  cycles  use  the  fuel  heated  through  cooling  of  the 
chamber  wall  as  the  working  fluid  for  the  turbine. 

The  turbine  powers  the  pumps  that  allow  the  propellants  to  be  stored  at  lower 
pressures,  their  by  reducing  the  structural  weight  of  the  tanks.  The  limiting  factor  to  the 
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performance  of  an  expander  cycle  is  the  turbine  inlet  temperature,  that  in  turn  limits  the 
attainable  chamber  pressure.  Consequently,  the  expander  engine  is  primarily  used  as  a 
space  engine  where  it  can  exhaust  to  a  vacuum  and  can  have  a  very  high  nozzle  area  ratio 
even  though  it  has  a  lower  chamber  pressure. 


2.2,  Thrust  Chamber 

2.2.1.  Rocket  Engine  Performance 

A  common  performance  parameter  used  to  define  a  rocket  engine  is  the  specific 
impulse  (fip).  The  fip  compares  the  thrust  of  the  engine  to  the  propellant  mass  flow  rate. 
Equation  1  show  the  fip  as  defined  by  Humble,  Henry  &  Larson:*’ 
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In  Equation  1 ,  E  is  the  thrust,  m  is  the  propellant  mass  flow  rate,  and  go  is  the 
acceleration  due  to  gravity.  The  units  for  fip  is  seconds  for  both  English  or  SI  units. 

In  order  to  describe  the  performance  of  each  component  of  the  thrust  chamber, 
two  coefficients  are  defined.  Eor  the  combustion  chamber,  the  characteristic  velocity  (c*) 
characterizes  the  influence  of  propellant  choice  through  absolute  maximum  temperature 
achievable.  Eor  the  nozzle,  the  thrust  coefficient  (Cj)  provides  the  conversion  of  the 
potential  energy  to  kinetic  energy  as  well  as  the  efficiency  of  the  nozzle  expansion. 
Equations  2  and  Equation  3  for  c*  and  Ct  comes  from  Hill  &  Peterson:^ 
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In  Equation  2,  y  is  the  ratio  of  specific  heats,  To  is  the  total  temperature  in  the 
chamber,  R  is  the  universal  gas  constant,  andM  is  the  molecular  weight  of  the 
propellant.  In  Equation  3,  pe,  Pa,  and  po  are  the  exit,  ambient,  and  total  chamber  pressure 
respectively,  Ag  is  the  exit  area  and  A*  is  a  characteristic  area  for  a  Mach  number  of 
one.  Erom  Equation  2,  c*  is  primarily  a  function  of  the  combustion  properties,  and  from 
equation  3,  Ct  is  primarily  a  function  of  nozzle  geometry.  To  demonstrate  the  influence 
of  the  fuel-oxidizer  composition  on  c*.  Table  2  shows  the  values  for  three  fuels  with 
liquid  oxygen. 


Table  2.  Influence  on  c*’ 


Oxidizer 

EO2 

EO2 

EO2 

Fuel 

EH2 

RPl 

CH4 

O/F 

4.83 

2.77 

3.45 

Toi  (K) 

3250 

3700 

3560 

Avs  bulk  density 
(kg/m^) 

320 

1030 

830 

c*  (m/s) 

2386 

1838 

1783 
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In  Table  2  LH2  is  liquid  hydrogen,  RPl  is  a  hydroearbon  fuel  and  CH4  is  liquid 
methane.  The  values  for  Table  2  come  from  Hill  &  Peterson.’  In  each  case  shown  in 
Table  2,  the  0/F  ratio  was  chosen  to  maximize  the  hp  and  the  chamber  pressure  was  6.89 
MPa.  As  apparent  from  the  data,  the  molecular  weights  of  the  propellants  and  the  values 
of  c*  depend  significantly  on  the  fuel-oxidizer  combustion. 

With  Equations  1  through  Equation  3,  a  simple  performance  analysis  can  be 
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conducted.  As  stated  by  Huzle  &  Huang,  the  calculation  of  thrust-chamber  performance 
is  based  on  the  theoretical  propellant  combustion  data  and  the  application  of  certain 
correction  factors.  The  theoretical  propellant  combustion  data  comes  from  thermo¬ 
chemical  computations  equating  the  heat  of  reaction  of  the  propellant  combination  to  the 
rise  in  enthalpy  of  the  combustion  gases.  The  desired  and  actual  performance  of  the 
chamber  calculated  from  these  equations  drives  the  design  of  the  chamber  presented  next. 

2.2.2.  Dimension  of  the  Thrust  Chamber 

The  geometry  of  the  thrust  chamber  is  based  on  the  pressure,  propellant  type, 
propellant  mass  flow  rate,  and  the  oxidizer-to-fuel  ratio  (0/E)  derived  from  the 
performance  analysis.  There  are  several  different  approaches  to  defining  the  geometry  of 
the  combustion  chamber.  Humble  et.al.^  begins  by  finding  the  throat  area.  The  mass 
flow  at  the  throat  must  be  chocked  for  proper  operation,  therefore  the  area  can  be 
determined.  Through  conservation  of  mass,  the  throat  area  (At)  can  be  found  from 
Equation  4; 
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To  insure  long  residence  times  for  the  mixing  and  chemical  reactions  in  the  thrust 
chamber,  the  Mach  number  must  be  nearly  zero.  The  low  Mach  number  implies  that  the 


thrust  chamber  pressure  is  nearly  the  stagnation  pressure.  Therefore,  the  chamber  area 
(Ac)  can  be  found  as  a  multiple  of  the  At,  as  determined  by  the  thermo-chemistry  of  the 
propellants.^  The  result  is  Equation  5; 
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For  Equation  5  the  Mach  number  is  generally  in  the  range  of  0. 1  to  0.6.  Humble 
et.al.^  also  derives  the  chamber  volume  (Vc)  from  a  ratio  of  the  At.  The  chamber  must  be 
large  enough  to  allow  complete  combustion.  However,  the  larger  the  chamber  the  greater 
the  weight.  Equation  6  uses  a  characteristic  length  (E*)  to  aid  in  the  sizing  of  the 
chamber: 
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In  Equation  6,  E*  is  the  thrust  chamber’s  characteristic  length.  Historical  data 
and  gas  dynamics  are  the  basis  for  the  sizing  of  E*.  Small  values  of  E*  imply  a  small 
engine.  The  goal  of  the  rocket  designer  is  to  minimize  the  size  and  mass  of  the  engine. 
Therefore,  E*  must  be  made  as  small  as  possible  while  maintaining  adequate  combustion 
efficiencies. 


2.2.3.  Injectors 

The  job  of  the  injector  is  to  atomize,  mix,  and  ignite  the  fuel  and  oxidizer. 
Ultimately,  the  injector  promotes  the  complete  combustion  of  the  liquid  propellants. 
Additionally,  the  pressure  loss  a  cross  the  injector  is  important  since  it  takes  a  way  from 
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the  pressure  obtained  in  the  ehamber.  To  better  understand  the  meehanisms  involved  in 
the  injeetors,  it  useful  to  see  what  work  has  been  done  in  this  area.  Rahman  &  Santoro^ 
reported  relevant  information  on  predieting  spray  drop-sizes  from  liquid  oxygen  and 
gaseous  hydrogen.  Figure  4  shows  a  sehematie  of  a  eoaxial  injeetor. 


Figure  4.  Shear  Coaxial  Gas/Liquid  Injector 


Atomization  is  of  interest  to  the  propulsion  eommunity  from  the  standpoint  of 
propellant  injeetor  design  for  liquid  roekets  engine  eombustion  ehambers.  Atomization 
from  a  shear  eoaxial  jet  refers  to  the  breakup  of  the  eore  liquid  jet  by  shear  forees  due  to  a 
eo-flowing,  high  veloeity,  annular  gas  jet  surrounding  the  eore  liquid  stream.^  Data 
presented  in  this  paper  did  not  contain  any  experiments  that  simulate  the  unique 
propellant  properties  of  liquid  oxygen  and  gaseous  hydrogen;  such  as  gaseous  hydrogen 
injection  density  and  liquid  oxygen  surface  tension  and  viscosity.  Simulation  of  these 
properties  is  crucial  to  obtaining  reliable  atomization  results  that  are  directly  relevant. 

Lightfoot,  Danczyk  &  Talley^*’  suggested  three  causes  of  atomization  from  wall- 
bounded  fdms:  liquid  turbulence,  stripping  of  waves,  and  stripping/tearing  resulting  from 
gas-phase  vortices.  Turbulent  eddies  within  the  liquid  can  interact  with  the  interface 
causing  it  to  become  roughened  and  eventually  forming  ligaments.  These  ligaments  may 
then  break  down  into  droplets.  Results  cite  several  important  non-dimensional 
parameters  including  ratios  of  film  height  to  hydraulic  diameter,  mean  velocity  to  RMS 
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velocity  fluctuations  as  well  as  a  liquid  Weber  number  based  on  hydraulic  diameter  and 
mean  surface  veIocity/°  Coherent  gas-phase  vortices  may  form  as  a  result  of  injector 
geometry  features  such  as  the  gap  or  the  lip.  If  the  liquid’s  energy  is  sufficiently  larger 
than  the  gas,  then  the  gas-phase  vortex  will  be  displaced;  otherwise,  the  vortex  will  alter 
the  path  of  the  liquid.  In  the  latter  case,  droplets  may  be  formed  when  the  vortex  distorts 
and  tears  liquid  away  from  the  film  or  when  aerodynamic  forces  arising  from  the  new 
shape  of  the  interface  strip  liquid  away.'°  Low  energy  ratio  simulations  show  waves 
form  on  the  interface  downstream  of  the  lip.  These  waves  are  uniform  and  grow  until 
they  reach  a  size  where  the  gas  flow  can  strip  mass  from  their  crests.  The  uniform  nature 
and  growth  of  these  waves  implies  hydrodynamic  instabilities  cause  the  droplet 
formation.  The  lip,  the  spacer  between  the  oxidizer  and  fuel  stream,  seems  to  have  very 
little  effect  on  the  general  atomization  behavior;  a  stronger  impact  is  expected  if  the 
gradient  of  taper  is  small  enough  that  the  flow  remains  attached  to  the  injector.  Current 
understanding  suggest  the  relative  momentum  difference  between  the  liquid  and  gas 
particularly  in  the  axial  direction  plays  a  large  role  in  the  film’s  behavior. At  high 
kinetic  energy  ratios,  gas-phase  structures  control  the  film’s  behavior  and  atomization. 

At  low  kinetic  energy  ratios,  waves  form  on  the  surface  of  the  film  and  are  responsible 
for  atomization.  Experiment  showed  atomization  due  to  turbulence  occurred  at  relatively 
large  gas  velocities;  suggesting  the  gas  turbulence,  not  liquid  turbulence,  was 
responsible.^*’ 

Besides  co-axial  injectors,  another  type  of  injector  that  displays  beneficial 
characteristics  is  the  Pintle  injector.”  Figure  5  shows  a  schematic  of  a  Pintle  injector. 
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fuel 


Figure  5,  Pintle  Injector 


Dressier  and  Bauer  found  a  Pintle  injector  design  could  deliver  high  combustion 
efficiency  and  enables  implementing  some  unique  operating  features,  such  as  deep 
throttling  and  injector  face  shutoff/'  Design  simplicity  makes  it  ideal  for  low  cost 
engines.  Significantly  lower  development  and  qualification  costs  result  because  injectors 
can  be  easily  adjusted  and  optimized  by  changing  only  two  simple  parts.  Additionally, 
there  has  never  been  an  instance  of  combustion  instability  in  a  Pintle  engine  during  any 
ground  or  flight  operations."  Either  fuel  or  oxidizer  can  be  centrally  metered  in  the 
Pintle  design.  The  Pintle  injector’s  flow  field  induces  recirculation  regions. 

Woodward  et.  al.  investigated  combustion  performance  of  coaxial  and  Pintle 
injectors  using  liquid  oxygen  and  ethanol.  Though  not  the  exact  same  as  the  oxygen  and 
hydrogen  used  for  the  DEAN,  the  results  should  show  similar  trends.  The  findings  shows 
marginal  liquid  oxygen  quality  was  seen  to  have  a  significant  influence  on  combustion 
efficiency  for  both  injector  types.  Pintle  injector  designed  specifically  for  orbital 
maneuvering  system  projects  operated  both  in  a  stable  manner  and  with  high 
performance.  Additionally,  changes  in  O/E  ratio  resulted  in  changes  in  c*  efficiency,  a 
point  to  consider  when  optimizing. 
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As  shown  in  this  section,  many  factors  affect  the  design  of  injectors,  however  the 
work  done  in  this  document  will  be  limited  to  the  effect  the  injector  will  have  on  the 
combustion  efficiency  and  pressure  drop  across  the  injector.  Still  this  information  will 
be  useful  in  determining  the  detailed  design  of  the  injectors  later. 

2,3.  Nozzle 

The  key  demands  on  future  launch  systems  are  the  reduction  of  Earth-to-Orbit 
launch  costs  in  conjunction  with  an  increase  in  the  reliability  and  operational  efficiency. 
Launch  systems  operate  in  a  constantly  changing  ambient  environment.  A  nozzle  that 
can  maintain  high  performance  over  a  wide  variety  of  ambient  conditions  could 
dramatically  improve  efficiency  and  decrease  cost.  .  The  DEAN  utilizes  one  such 
concept  known  as  the  aerospike  or  plug  nozzle.  A  booster-stage  would  experience  the 
greatest  change  in  ambient  pressure  and  could  benefit  most  form  the  aerospike  design. 
However,  an  upper-stage  can  benefit  as  well  since  it  can  encounters  a  change  in  ambient 
pressure  depending  on  where  it  ignites.  Eigure  6  show  an  aerospike. 

Flow  Field 
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The  aerospike  nozzle  is  considered  to  have  globally  better  performance  when 
compared  to  a  conventional  bell  nozzle,  since  the  expansion  of  the  jet  is  not  bounded  by  a 
wall  and  therefore  can  adjust  itself  to  the  environment  by  changing  the  outer  jet 
boundary. 

2.3.1.  Aerospike 

Hagemann  et.al.  examined  several  nozzle  concepts  that  could  result  in  superior 
performance  over  conventional  nozzles.  The  advantage  of  aerospike  nozzles  is  it 
demonstrate  altitude  adaptation  up  to  their  geometrical  area  ratio.  This  results  in 
improved  performance  for  the  entire  flight  envelope  over  the  conventional  bell  nozzle. 
Additionally,  for  an  aerospike  nozzle,  at  lower  pressure  ratios  an  open-wake  flow 
established  at  a  pressure  level  practically  equal  to  the  ambient  pressure.  At  a  specific 
pressure  ratio  (PR),  close  to  the  design  PR  of  a  full-length  nozzle,  the  base  flow  suddenly 
changes  its  character  and  turns  over  to  a  closed  form.  A  constant  base  pressure  no 
longer  influenced  by  ambient  pressure  characterizes  this  flow.  Shorter  plug  nozzles  with 
high  truncations  trigger  an  earlier  change  in  wake  flow.  At  transition,  the  pressure  within 
the  wake  approaches  a  value  below  the  ambient  pressure  and  the  full  base  area  induces  a 
negative  thrust.  Beyond  the  transition  point,  pressure  within  the  closed  wake  remains 
constant,  as  ambient  pressure  decreases.  The  base  pressure  is  then  higher  than  the 
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ambient  pressure,  resulting  in  a  positive  thrust  contribution. 

A  more  detailed  analysis  of  the  aerospike  concept  was  done  by  Hagemann, 
Immich  &  Terhardt^"^  using  the  numerical  methods  of  Euler  and  Navier-Stokes.  The 
results  showed  the  altitude  compensation  capabilities  of  the  aerospike  were  indisputable, 
but  they  lose  this  capability  for  pressure  ratios  above  the  design  point. Additional 
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performance  losses  were  induced  due  to  non-isentropic  effects  like  shock  waves.  Figure 
7  shows  the  flow  over  an  aerospike  nozzle.  In  Figure  7,  the  top  picture  depicts  the  flow 
field  at  PR  below  design  (lower  altitude),  the  middle  picture  is  at  the  designed  PR,  and 
the  lower  picture  is  at  a  PR  above  the  design  point  (higher  altitude). 


e^^aitsion  waves 


shear  layer 


Figure  7,  Flow  Phenomena  of  Plug  Nozzle 

As  Figure  7  illustrates,  for  pressure  ratios  lower  than  the  design  PR  of  a  plug 
nozzle  with  a  well-contoured  body,  the  flow  primarily  expands  along  the  plug  body  to  the 
ambient  pressure.  Thus,  only  the  first  part  of  the  nozzle  contour  acts  as  an  expansion 
contour  down  to  the  point  where  the  first  right  running  characteristic  that  feels  the 
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ambient  pressure  meets  the  contour.'"^  At  the  design  pressure  ratio,  the  characteristic  with 
the  design  Mach  number  should  be  a  straight  line  emanating  to  the  tip  of  the  central  plug 
body,  and  the  shear  layer  is  parallel  to  the  centerline.  At  pressure  ratios  above  the  design 
pressure  ratio,  the  wall  pressure  distribution  remains  constant,  and  the  plug  nozzle 
behaves  as  a  conventional  nozzle,  the  loss  of  its  capability  of  further  altitude  adaptation 
included.'"^ 

Truncation  of  the  nozzle  makes  the  aerospike  concept  more  feasible,  but  results  in 
further  performance  losses.  Figure  8  shows  the  flow  over  a  truncated  aerospike.  In 
Figure  8  the  top  drawing  depicts  the  flow  field  at  PR  below  design,  the  middle  drawing  is 
at  the  designed  PR,  and  the  lower  drawing  is  at  a  PR  above  the  design  point. 

interaction  oT coitq»  ^waves 


shear  layer 


Figure  8,  Flow  Phenomena  of  Truncated  Plug  Nozzle 
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As  seen  in  Figure  8,  at  lower  pressure  ratios  an  open  wake  flow  establishes,  with  a 
pressure  level  praeticahy  equal  to  the  ambient  pressure.  At  a  eertain  pressure  ratio  elose 
to  the  design  pressure  ratio  of  the  full-length  plug  nozzle,  the  base  flow  suddenly  ehanges 
its  eharacter  and  turns  over  to  the  closed  form,  characterized  by  a  constant  base  pressure 
that  is  no  longer  influenced  by  the  ambient  pressure. At  the  transition  point,  the 
pressure  within  the  wake  approaches  a  value  that  is  below  ambient  pressure,  and  the  full 
base  area  induces  a  negative  thrust.'"^  Beyond  the  pressure  ratio  at  transition,  the  pressure 
within  the  closed  wake  remains  constant.  A  further  decrease  of  the  ambient  pressure  is 
resulting  therefore  in  a  positive  thrust  contribution  of  the  total  base  area.'"^ 

Ambient  flow  also  slightly  degrades  performance.  Ah  these  losses  add  to  several 
percent,  making  plug  perform  worse  at  high  altitude  when  designed  for  a  lower  altitude 
when  compared  to  similar  bell  nozzles.'"^  To  get  the  most  benefit,  the  design  pressure 
ratio  should  be  chosen  at  as  high  an  altitude  (pressure  ratio)  as  possible. 

Concerning  the  truncation,  Ito,  Fujii  &  Hayashi'^  sate  the  nozzle  performance  was 
believed  not  to  be  strongly  affected  by  cutting  off  the  nozzle  because  base  pressure 
compensated  for  the  loss  of  the  thrust.  For  a  full-length  nozzle,  flow  follows  the  nozzle 
wall  and  smoothly  moves  to  the  wake  region.  For  a  truncated  nozzle,  flow  separates  at 
the  trailing  edge  and  expands.  The  resultant  shear  layer  induces  a  trailing  shock  when  it 
converges  and  turns  the  flow  parallel  to  the  axis.  When  truncated,  the  ramp  area 
decreases  due  to  the  shorter  length.  Therefore,  thrust  from  the  ramp  pressure  reduces. 

The  thrust  generated  by  base  pressure  increases  due  to  the  increased  area.  Therefore,  the 
total  nozzle  thrust  becomes  almost  the  same  for  any  nozzle  truncation. Quantitatively, 
the  thrust  coefficients  of  the  plug  nozzles  have  the  same  trend  as  the  ideal  thrust 
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coefficient  indicating  the  plug  nozzle  operated  at  near  peak  thrust  efficiency  over  a  wide 
pressure  ratio  range.  As  the  altitude  increases,  ambient  pressure  decreases,  therefore  the 
base  pressure  thrust  increases.  Ito  et.al.  have  found  a  contoured  plug  nozzle  had  higher 
performance  compared  to  a  conical  plug  nozzle  over  all  pressure  ratios.  As  the  pressure 
ratio  increases,  the  thrust  difference  between  the  contoured  and  conical  nozzle  increases 
almost  linearly.  The  advantage  of  a  contoured  nozzle  becomes  remarkable  as  the  altitude 
increases.  External  flow  does  not  influence  the  pressure  distribution  on  the  nozzle 
surface  at  high-pressure  ratios.'^ 

In  further  work  done  by  Ito  &  Fujii'^,  they  found  that  in  the  low-pressure  ratio 
regions,  the  stagnation  point  was  located  at  a  distance  from  the  plug  base.  When  the 
pressure  ratio  increases,  the  stagnation  point  suddenly  moves  toward  the  plug  base.  The 
stagnation  point  does  not  move  with  further  pressure  ratio  increases.  The  pressure  ratio 
at  which  an  abrupt  movement  occurs  is  the  same  as  the  pressure  ratio  where  the 
characteristic  transition  occurs.  The  ambient  pressure  influences  base  pressure  as  long  as 
the  flow  initiated  from  the  stagnation  point  due  to  the  envelope  shock-wave  impingement 
reached  the  base.  They  conclude  the  aerospike  nozzle  performance  is  insensitive  to 
length  of  nozzle  because  the  base  pressure  compensates  the  loss  of  the  thrust  force  due  to 
nozzle  truncation.  The  characteristics  of  the  flow  at  the  base  did  not  change  whether 

the  external  flow  was  induced  or  not.  The  ambient  pressure  influences  the  base  pressure 
when  the  pressure  ratio  is  low  and  the  base  pressure  becomes  independent  from  ambient 
when  the  pressure  ratio  is  high.  The  base  produces  positive  thrust  when  the  base  pressure 
is  independent  from  ambient  pressure. 
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In  contrast  to  Ito  &  Fujii  conclusion  on  external  flow  effeets,  Nasuti  &  Onofri'^ 
found  the  numerieal  analyses  they  performed  indieates  the  interaetion  of  the  exhaust  jet 
with  the  external  air  in  truneated  plug  nozzles  may  signifieantly  affeet  the  nozzle 
behavior  in  the  over-expanded  regime.  In  partieular,  even  in  the  ideal  ease  important 
losses  take  plaee  yielding  a  sudden  performanee  drop.  The  presence  of  a  finite-thiekness 
or  thiek  shroud  substantially  ehanges  the  flow  behavior  adding  further  drag  in  itself  and 
yielding  an  overexpansion  at  the  primary  nozzle  lip  (edge  of  the  outer  wall).  Beeause  of 
this  overexpansion,  the  exhaust  jet  finds  a  lower  ambient  pressure  at  the  lip  than  in  the 
still  air  ease,  and  eonsequently  adapts  the  flow  to  a  lower-than-ambient  pressure.  This 
was  a  further  eause  of  drag  appearing  in  part  over  the  plug  wall  and  in  part  over  the  plug 
base.  Moreover,  the  reduetion  of  the  value  of  pressure  ratio  for  transition  from  open  to 
closed  wake  redueed  the  overall  nozzle  performanee.  This  reduction  yielded  a  base  drag 
in  the  range  between  the  aetual  and  still  air  transition  values.  Coneerning  the  mechanism 
of  wake  transition,  in  the  supersonie  ease  the  transition  is  governed  by  the  internal  shoek, 
rather  than  by  the  last  wave  of  the  expansion  fan  at  the  primary  nozzle  lip.^’  On  the 
eontrary,  the  analysis  of  the  flow-field  in  the  subsonie  flight  eonditions  has  shown  only 
slight  ehanges  were  expeeted  in  eomparison  to  the  still  air  ease. 

Ito  &  Fujii  also  performed  a  test  to  determine  the  effeets  of  base  bleed.  They 
found  the  base  pressure  inereased  due  to  an  inerease  in  base  area  and  eompensated  for  the 
total  thrust  loss  due  to  a  deerease  of  the  ramp  pressure  thrust.  As  a  result,  the  total  nozzle 
thrust  beeomes  almost  the  same  for  any  nozzle  truneation.  When  they  introdueed  base 
bleed,  it  expanded  toward  low  pressure  and  interaeted  with  the  main  exhaust  flow  at 
some  distanee  from  base  surfaee.  The  base  bleed  that  promoted  reeireulation  at  the  base 
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region  produced  the  largest  pressure  thrust.'^  Additionally,  base  bleed  producing  no 
divergent  loss  leads  to  the  maximum  total  thrust.  At  low  altitude  (low  PR)  base  pressure 
linearly  decreases  as  the  ambient  pressure  decreases;  showing  the  external  environment 
influences  the  base  region.  Pressure  thrust  produced  by  base  in  this  region  is  small.  At 
high  altitude  (high  PR),  the  base  pressure  is  constant  despite  a  variation  in  altitude 
(decrease  in  ambient  pressure).  The  base-pressure  thrust  increased  in  this  region.  The 
base  pressure  with  and  without  base  bleed  acted  in  a  similar  manner.  Conditions  with 
bleed  had  a  higher-pressure  level  than  without  bleed.  The  thrust  coefficients  with  and 
without  bleed  had  the  same  trends  as  an  ideal  thrust  coefficient  exhibits  indicating  the 
aerospike  operates  at  nearly  peak  thrust  efficiency  over  a  wider  range  than  a  bell  nozzle. 
The  thrust  performance  with  base  bleed  exceeded  the  performance  without  base  bleed 
over  the  whole  altitude  range.  The  base  bleed  orientation  with  the  greatest  performance 
was  at  the  outer  region  of  the  base,  directed  parallel  to  the  nozzle  axis. 

2,4.  Cooling  Jacket 

Cooling  of  the  thrust  chamber  and  nozzle  is  essential  in  the  design  of  any  rocket 
engine.  Many  different  cooling  approaches  exist  and  some  are  outlined  below. 

2.4.1.  Heat  Transfer 

To  describe  the  environment  in  the  thrust  chamber,  heat  transfer  equations  are 
needed.  An  excellent  source  to  provide  an  understanding  of  heat  transfer  is  Incropera 
et.al.  Fundamentals  of  Heat  and  Mass  Transfer.  The  two  major  mechanisms  in  which 
the  heat  is  transferred  are  through  conduction  and  convection.  Conduction  is  the 
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transport  of  energy  in  a  medium  due  to  a  temperature  gradient.'^  A  basie  equation  to 
determine  the  eonduetive  heat  rate  is  shown  in  Equation  7: 

q,=-kA^  {!) 

ax 

In  Equation  7  qx  is  the  heat  rate  with  units  of  watts  (W),  k  is  the  thermal 
conduetivity  with  units  W/(m*K),  A  is  the  area,  and  dX/dx  is  the  derivative  of  the 
temperature.  The  thermal  eonduetivity  is  a  material  eonstant. 

Conveetion  deseribes  the  energy  transfer  between  a  surfaee  and  a  fluid  moving 
over  the  surfaee.'^  A  basie  equation  to  determine  the  eonveetion  heat  flux  is  shown  in 
Equation  8: 

q"=h(T,-T,)  (8) 

In  Equation  8  q”  is  the  heat  flux  with  units  of  W/m  ,  h  is  the  eonveetive  heat 
transfer  eoeffieient  with  units  W/(m  K),  and  Tg  and  To,  are  the  temperature  of  the  surfaee 
and  fluid  respeetively.  Unlike  the  thermal  eonduetivity,  whieh  is  a  material  eonstant,  the 
eonveetive  heat  transfer  eoeffieient  depends  on  numerous  fluid  properties,  surfaee 
conditions  and  flow  conditions. 

Some  of  the  constants  used  in  the  heat  transfer  equations  have  to  be 

experimentally  determined.  Von  Glahn  proposed  correlation  treating  nozzle  heat 

transfer  with  a  cooled  approach  section.  This  concept  applies  to  both  low  and  high 

Reynolds  number  regions.  It  also  applies  to  the  entire  nozzle,  whether  convergent  or 

divergent.  The  proposed  correlation  consists  of  a  fully  turbulent  pipe-flow  heat-transfer 
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equation  modified  by  suitable  nozzle  geometry  parameters. 
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Schacht,  Quentmeyer,  &  Jones^'  conducted  a  gas-side  heat- transfer  experiment 
over  a  wide  range  of  ehamber  pressures  with  an  emphasis  on  the  aceurate  determination 
of  hot  gas-side  heat-transfer  rates  in  the  regimes  of  high  heat  flux.  Equation  9  adequately 
correlates  the  data  at  speeific  locations  in  the  nozzle: 

5t*Pr*°"  =C-Re*/'"  (9) 

In  Equation  9  St*,  Pr*,  and  Re*d  are  Stanton,  Prandtl,  and  Reynolds  number 
respeetively  and  C  is  a  eonstant.  The  eonstant  varies  with  axial  location  and  is  less  than 
0.026  for  all  locations  except  the  chamber.  The  eonstant  at  the  throat  was  42%  less  than 
the  widely  used  value  of  0.026. 

Boldman,  Sehmidt,  &  Gallagher  showed  the  heat  transfer  at  a  given  station  in 
nozzles  generally  exhibits  two  distinct  depressions  from  the  predicted  levels  based  on  a 
turbulent  pipe  flow  type  of  eorrelation.  The  larger  of  these  depressions  oeeurs  at  lower 
Reynolds  numbers  and  was  assumed  the  result  of  laminarization  of  the  initially  turbulent 
boundary  layer.  The  smaller  of  these  depressions,  whieh  oecurred  at  high  Reynolds 
numbers,  was  assumed  the  result  of  reduced  turbulent  transport  assoeiated  with  a 
turbulent  boundary  layer  in  an  aeeelerating  flow. 

23 

An  experimental  investigation  conducted  by  Quentmeyer  &  Roncace 
determined  the  hot-gas-side  heat  transfer  characteristics  for  a  liquid-hydrogen-cooled, 
subscale,  plug-nozzle  roeket  test  apparatus.  The  throat  region  eorrelation  coeffieient  (Cg) 
for  a  Nusselt  number  correlation  of  the  form  in  Equation  10: 

iVM  =  CgRe“*Pr“'  (10) 
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The  averaged  value  is  0.023  for  the  Rigimesh  faeeplate  and  0.026  for  the  platelet 
faeeplate. 

2.4.2.  Cooling  Techniques 

There  are  several  teehniques  used  to  eool  the  ehamber  and  nozzle.  The  most 
eommon  inelude  regenerative  eooling,  dump  eooling,  fdm  eooling,  transpiration  eooling, 
and  ablative  eooling.  Regenerative  eooling  is  most  relevant  to  the  DEAN  and  the 
teehnique  further  studied. 

Preelik  et.  al.^"^  showed  the  development  of  the  wall  heat  flux  is  not  a  quiek 
proeess,  but  spreads  out  from  the  injeetor  down  to  the  eonvergent  nozzle  entranee.  The 
ealorimeter  data  also  highlighted  the  flow  and  burning  oharaeteristies  of  the  present,  eo- 
axial  injeetion  elements.  These  elements  affeet  the  level  of  the  overall  wall  heat  load 
approximately  25%.  Two  prineiples  for  redueing  the  hot  gas  side  heat  transfer  rates,  i.e. 
wall  element  mixture  ratio  trimming  and  gaseous  hydrogen  wall  film  eooling,  were 
investigated  in  more  detail.  The  tests  elearly  indieated  wall  film  eooling  employing  the 
eurrent  slot  design  is  mueh  more  effieient  and  effeetive  when  eompared  to  wall  element 
mixture-ratio  trimming.  With  typieal  film  eooling  mass  flow  rates,  the  wall  heat  loads 
eould  be  deereased  by  20%.  Moreover,  under  the  presenee  of  a  eoolant  film,  the 
individual  flow  oharaoteristies  of  the  different  injeetion  elements  beeome  less  important. 

Kim  et.  al.  eondueted  a  test  of  a  regeneratively  eooled  ehamber  with  film 
eooling.  From  the  experimental  results  with  the  film  eooling  mass  flow  rate  at  10.5%  of 
the  main  fuel  mass  flow  rate,  maximum  heat  flux  at  the  nozzle  throat  was  measured  to  be 
30%  lower  than  without  film  eooling  at  the  nominal  operating  eondition.  The  numerieal 
analysis  resulted  in  this  ease  showing  a  13%  deerease.  The  film  eooling  mass  flow  rate 
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increased  as  the  eharaeteristie  veloeity  deereased.  When  the  mass  flow  rate  of  fllm 
eooling  was  10.5%  of  the  total  fuel  mass  flow  rate,  the  eharaeteristie  veloeity  is  measured 
to  be  1 .2%  lower  than  without  film  eooling  at  the  design  point. 

In  researeh  by  Naragi,  Dunn,  &  Coats,  the  effeetiveness  of  regenerative  dual- 
eireuit  eooling  designs  was  studied.  The  dual-eireuit  showed  redueed  wall  temperatures  at 
the  throat  area  for  the  spaee  shuttle  main  engine  (SSME)  and  for  an  RPl-LOX  engine.  A 
lower  eoolant  pressure  drop  aeeomplished  the  reduetion  in  wall  temperature  and  resulted 
in  lower  eoolant  pumping  power  requirements.  The  overall  result  is  a  higher  performing 
engine  eapable  of  delivering  a  greater  payload  to  orbit. 

Chiaverini,  Sauer,  &  Munson^^  used  a  vortex  eombustion  eold-wall  ehamber 
(VCCW)  to  study  the  heat  transfer  problem.  They  found  the  ehamber  sidewall  heating 
rates  did  not  display  a  signifieant  dependenee  on  ehamber  pressure,  apparently  due  to  the 
similar  effeets  of  elevated  pressure  on  both  thermal  radiation  (aeting  to  heat  the  wall)  and 
eonveetive  eooling  from  the  outer  vortex. 

Immieh  &  Caporieei  eondueted  an  investigations  within  the  FESTIP  teehnology 
program  on  the  high  aspeet  ratio  eoolant  ehannel  and  miero  fin  eooling  ehannel  ehamber. 
These  teehnologies  have  been  seleeted  as  promising  teehnologies,  whieh  ean  be  expeeted 
to  lead  to  eonsiderable  reduetions  of  the  eombustion  ehamber  wall  temperature  by 
improving  the  eoolant  side  heat  transfer  signifieantly.  This  eould  be  aehieved  by  an 
inereased  heat  transfer  surfaee  and  by  taking  partieular  advantage  of  an  enhaneed  fin 
effeet.  The  eoneept  was  to  inerease  the  total  number  of  eooling  ehannels  and  to  introduee 
mierostruetures  at  the  eoolant-side  ehannel  bottom. 
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Schmidt,  Popp,  &  Frohlich^^  stated  an  improvement  of  eoolant  side  heat  transfer 
and  eonsequently  a  deerease  of  the  ehamber  wall  temperature  has  only  very  minor  effeet 
on  the  inerease  of  heat  flux  to  the  eoolant,  beeause  the  driving  temperature  differenee 
eannot  be  ehanged  mueh  by  redueing  the  wall  temperature.  A  doubling  of  the  eoolant- 
side  heat  transfer,  by  higher  eoolant  veloeity,  eoolant  side  surfaee  roughening  or  fins  in 
the  eooling  ehannel  ete.,  would  only  inerease  the  heat  flux  to  the  eoolant  by  about  6-7%. 
Inereases  of  heat  flux  in  the  ehamber  by  lower  chamber  diameter  and  eonsequently 
higher  heat  flux  is  another  possibility.  Redueing  the  ehamber  eontraetion  ratio  eould 
inerease  the  heat  flux  with  a  resulting  10-15  K  (18-27  R)  eoolant  temperature  rise,  but 
only  if  the  ehamber  is  elongated  to  keep  L*  eonstant,  whieh  also  inereases  the  engine 
length.^^  The  assoeiated  inerease  in  eooling  and  eombustion  pressure  drop  would  also 
inerease  the  required  hydrogen  pump  diseharge  pressure.  This  again  would  reduee  the 
overall  benefit  for  the  eyele  to  a  marginal  improvement.  On  the  other  hand,  from  the 
injeetor  side  a  certain  faee  area  is  required,  also  limiting  the  possibility  of  employing 
such  a  measure.  Roughening  of  the  hot  gas  side  surfaee  of  the  ehamber  and  nozzle  would 
also  inerease  the  heat  flux  into  the  eoolant.  This  has  been  demonstrated  by  the  aging 
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effeet  of  ehambers,  running  at  their  thermal  limits.  The  most  straightforward  design 
measure  for  inereased  heat  transfer  is  the  inereased  surfaee  roughness  of  the  hot  gas  side 
surfaee. 

2.4.3.  Aerospike  Cooling 

The  aerospike  adds  additional  ehallenges  when  it  eomes  to  eooling  the  nozzle 
sinee  it  is  surrounded  by  the  exhaust  flow,  whereas  the  Bell  nozzle  is  exposed  to  the 
ambient  surroundings  on  one  side.  In  Sorge,  Carmieino,  &  Noeito  a  eooled  linear  plug 
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nozzle  was  constructed  so  hot  gas  tests  could  be  carried  out  and  performance  measured. 
They  chose  a  two-dimensional  nozzle  shape  and  a  copper  alloy  on  all  the  parts  exposed  to 
hot  gas  flow.  In  addition,  to  promote  the  heat  transfer,  ribs  were  placed  on  the  duct 
surfaces  exposed  to  the  heat  flux  near  the  throat.  The  results  have  showed  the  experiment 
was  able  to  ensure  a  maximum  wall  temperature  on  the  gas  side  of  750  K  (1350  R),  a 
maximum  wall  temperature  on  the  liquid  side  of  450  K  (810  R)  and,  a  water  bulk 
temperature  rise  of  about  298  K  (537  R),  with  a  pressure  drop  per  circuit  equal  to  31 
atm.  Even  though  this  experiment  uses  water,  the  technique  can  be  applied  to  any 
cooling  fluid. 

Kumakawa  et.al.^^  conducted  a  test  with  nickel  plating  on  the  hot-gas  side.  This 
test  showed  the  heat  flux  was  30%  lower  than  without  nickel.  A  cross-flow  condition  of 
water  suppressed  the  burnout  on  the  ribbed  surface  more  than  the  parallel  flow  condition. 
Combustion  performance  and  heat  transfer  characteristics  for  the  truncated  conical  plug 
nozzle  were  quite  similar  to  those  of  conventional  conical  nozzles. 

From  Wang  the  highest  level  of  heating  occurs  near  the  thruster  outlet,  on  the 
ramp  surface.  The  next  highest  level  of  heating  occurs  on  the  plug-base  where  the 
reverse-jet  brings  in  the  hot  plume  gases  torch  the  surface.  The  relief  comes  from 
aspirating  cold  inner-base  airflow  into  the  plug-base  region.  The  airflow  penetrates  as  far 
as  the  base  center.  The  pumping  effect  of  the  engine  plume  causes  this  aspiration.  The 
heating  on  the  plug  sidewall  is  caused  by  the  hot  engine  flow  spilling  off  the  side  of  the 
ramp,  also  known  as  the  plume  spillage  heating. 

Tsutsumi  et.al.  observed  the  maximum  heat-flux  region  is  located  where  the  jets 
from  adjacent  modules  interact.  The  heat  flux  also  increased  at  the  module  exit.  The 
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downstream  part  of  the  module  exit  is  exposed  to  a  lower  heat  load.  The  heat-flux 
distribution  over  the  linear  aerospike  nozzle  is  eorrelated  with  the  viseous-inviseid 
interaction  in  the  near-surface  flow  field.  The  heat  flux  increases  remarkably  at  the  region 
where  the  module  outflows  interact  due  to  the  downwash  of  module  outflow  to  the  nozzle 
surface.  Moreover,  since  the  ambient  air  is  entrained  into  the  vortex  flow  over  the  nozzle 
surface,  a  lower  heat-flux  region  appears. 

2.4.4.  Material 

Selection  of  the  materials  for  the  cooling  jacket  tends  to  drive  many  aspects  of  the 
design  process.  Peer  &  Minick^"^  describes  an  advanced  combustion  chamber  designed 
by  Pratt  &  Whitney.  The  most  significant  feature  is  the  use  of  a  new  copper  alloy 
coupled  with  an  improved  processing  technique.  Their  tubes  can  withstand  repeated 
exposure  to  fabrication  temperatures  in  excess  of  1250  K  (2260  R)  and  still  retain  yield 
strength  five  times  greater  than  other  copper  alloys  used  in  current  rocket  thrust 
chambers. The  tubular  configuration  of  the  chamber  provides  up  to  40  percent  more 
actual  surface  area  due  to  the  circular  tube  crowns  and  therefore,  more  heat  transfer 
capability  with  lower  thermal  strain  (increased  life)  than  smooth  wall  hot-side  fabricated 
channel  configurations^"^.  The  tubular  construction  also  provides  improved  pressure  drop 
characteristics  over  rectangular  channel  designs. 

Schoenman  discusses  the  intensive  effort  to  develop  materials  that  can  operate  at 
higher  temperatures  over  the  last  decade.  A  nominal  target  value  of  2500  K  (4460  R)  has 
been  established  for  the  new  materials  of  construction.  The  iridium/rhenium  layered  wall 
and  the  silicon/carbon  composite  are  the  most  mature  and  can  be  considered  for  flight 
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applications.  Iridium/rhenium  materials  can  operate  in  temperatures  2477  K  (4459  R) 
while  silicon/carbon  can  operate  in  1866  K  (3359  R)  temperatures.^^ 

Kumakawa  et.  al.^^  showed  carbon-carbon  and  silicon/carbon  composites  have 
potential  for  use  as  lightweight  materials  exposed  to  temperatures  in  excess  of  1700  K 
(3060  R).  The  exposed  heat  fluxes  dictate  whether  active  cooling  would  be  required  for 
these  materials.  Aerospike  nozzles  made  of  composite  would  be  exposed  to  severe 
operating  conditions  in  a  hot  gas  flow,  including  hot  spots  caused  by  shock  interactions. 
Thus,  an  actively  cooled  composite  is  necessary  in  this  case  and  considered  a  prime 
candidate  for  the  aerospike-nozzle  wall  structure 

A  significant  challenge  will  be  choosing  materials  to  accommodate  the  liquid 
oxygen  expander  cycle  of  this  concept.  Gloyer,  Knuth,  &  Crawford^^  designed  a  liquid 
oxygen  cooled  gas  generator.  They  found  copper,  nickel  and  some  steels  can  be  used 
adequately  with  oxygen.  These  metals  are  used  since  they  do  not  react  with  oxygen  that 
would  cause  catastrophic  failure.  They  also  concluded  the  design  appears  to  be  limited  to 
a  single  operating  point  design  due  to  the  liquid  oxygen  boiling  sensitivity  to  factors  such 
as  0/F  ratio,  mass  flow,  temperature  and  pressure.  For  wider  operating  ranges,  alternative 
designs  which  avoid  two-phase  flow  would  be  preferred. 

2.4.5.  Aspect-Ratio 

Changing  the  area  of  the  chamber  affects  the  heat  transfer  characteristics.  This 
can  be  done  through  the  aspect-ratio  (AR)  of  the  cooling  circuit  and/or  the  addiction  of 
ribs.  The  heat  energy  requirements  of  the  turbo-pumps  dictate  longer  combustion 
chambers.  Size  limitations  create  the  need  for  a  different  method  to  increase  heat 
extraction.  Increasing  the  area  exposed  to  the  hot-gas  by  using  combustor  ribs  fulfilled 
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this  requirement.  The  ribs  increased  the  total  area  exposed  to  the  hot-gas  by  80%  in  his 
research,  and  thus  enhanced  the  heat  energy  level  imparted  to  the  coolant  working  fluid^^. 

An  investigation  was  conducted  by  Carlile^^  to  determine  the  effectiveness  of 
using  high-aspect-ratio  cooling  passages  to  improve  the  life  and  reduce  the  coolant 
pressure  drop  in  high-pressure  rocket  thrust  chambers.  Figure  9  show  the  difference 
between  a  low  and  high  aspect  ration  channel  (AR). 


Low.kR  Higli.AR 


^  ^  ^  ^  ”1^  ^ 

Heat 

Figure  9,  Channel  Aspect  Ratio 
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Figure  9  shows  that  the  low  AR  channel  is  generally  wider  with  more  surface 
area  along  the  wall  while  the  high  AR  is  taller  with  less  surface  area.  The  coolant 
pressure  drop  for  the  high-aspect-ratio  chamber  was  reduced  in  increments  to  one -half  the 
baseline  chamber  by  reducing  the  coolant  mass  flow.  The  result  still  showed  a  reduction 
in  the  hot-gas-side  wall  temperature.  The  data  indicated  the  hot-gas-side  wall 
temperatures  for  the  high-aspect-ratio  chamber  could  have  been  reduced  substantially 
further  by  using  aspect  ratios  greater  than  5.0. 

Wadel  &  Meyer"^*’  found  by  increasing  the  cooling  channel  surface  area  through 
increasing  aspect  ratios,  heat  from  the  hot-gas-side  wall  is  more  efficiently  transferred  to 
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the  coolant.  The  increased  height  and  number  of  the  ribs  also  enhance  the  heat  transfer 
from  the  chamber  liner  to  the  coolant  (i.e.  enhanced  fin  effect).  Therefore,  it  is  possible  to 
fabricate  chambers  with  sufficiently  greater  total  flow  area  to  reduce  pressure  drop  over  a 
conventional  design,  and  still  gain  an  increase  in  the  heat  transfer  capability. 

In  Neuner  et.al."^^  paper,  laboratory  experiments  with  large  scale  cooling  channel 
models  were  described.  The  purpose  of  the  experiments  was  to  determine  the  impact  of 
high-aspect  ratio  channels  on  curvature  induced  heat  transfer  enhancing  phenomena.  As 
a  general  result,  secondary  flow  structures  were  clearly  identified,  even  in  channels  with 
an  aspect  ratio  of  8.0.  These  vortex  phenomena  have  been  found  to  enhance  the  heat 
transfer  both  in  convex  and  concave  side  heated  bends.  Nevertheless,  they  appear  only  in 
and  not  far  downstream  of  the  corresponding  curvatures. 

The  results  of  a  parametric  study  on  cooling  channel  geometry  showed  as  the 
channel  geometry  changes,  the  coolant  heat  transfer  coefficient  dominates  the  heat 
transfer  rate  as  compared  to  the  area  terms. In  general,  a  small  flow  area  tends  to 
increase  heat  transfer;  however,  the  pressure  loss  across  the  coolant  channel  restricts  the 
extent  the  area  can  be  reduced.  Higher  aspect  ratio  cooling  channels  are  advantageous  in 
balancing  the  pressure  loss  requirements  with  the  heat  transfer  demands.  The  maximum 
allowable  chamber  pressure  is  limited  by  the  survivable  gas  wall  temperature.  With 
respect  to  pressure  loss  and  over  all  engine  mass,  a  stepped-channel  configuration  proves 
to  be  superior  to  an  invariant  channel  while  maintaining  the  thermal  performance  of  the 
regenerative  cooling  jacket.  Evidence  of  this  was  presented  by  Schuff  et.al.  and  is  shown 
in  Table  3. 
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Table  3,  Stepped  Channel  vs.  Constant  Cross-Section  Channel'*^ 


Constant  Cross-  Section 

Stepped 

Pin  (Psi) 

2299 

1838 

Pout  (psi) 

1762 

1680 

APe  (psi) 

537 

158 

ATe  (R) 

766 

737 

Engine  Mass 
(Ibm) 

327 

274 

Table  3  shows  a  APcand  All  decreased  of  340%  and  4%,  respectively,  from  the 
constant  cross-section  to  the  stepped  channel.  Schuff  et.al.  state  near  the  throat  the  gas 
side  heat  transfer  coefficient  increases  significantly  and  as  a  result  the  wall  temperature 
also  increases.  To  adequately  deal  with  the  high  heat  transfer  rate,  the  geometry  remained 
the  same  as  defined  for  the  baseline  case  in  this  region.  In  the  reduced  combustion  gas 
heat  transfer  regions  upstream  and  downstream  of  the  throat  we  increased  the  channel 
cross-sectional  area  to  decrease  the  pressure  drop  and  maintain  a  combustion  gas  side 
wall  temperature  less  than  at  the  throat.  Consequently,  increasing  the  channel  cross- 
sectional  area,  decreased  the  coolant  heat  transfer  coefficient  he,  therefore  the  stepped- 
channel  case  presented  a  tradeoff  in  performance  parameters,  temperature  rise  (AJ^  and 
pressure  drop  (APe)  of  the  coolant  along  the  entire  length  of  the  channels. 
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3.  Methodology 


This  section  will  outline  the  methodology  implemented  to  design  the  thrust 
chamber,  nozzle,  and  cooling  jacket  of  the  DEAN.  To  begin  the  process,  the 
performance  goals  outlined  in  the  introduction  establish  a  baseline  engine  derived  from 
basic  rocket  equations. 

3.1.  General  Sizing 

The  design  started  with  the  determination  of  a  solution  space  for  the  size  of  the 
engine.  The  size  of  the  engine  begins  with  the  choosing  of  initial  parameters;  the  fuel 
and  oxidizer,  the  oxidizer-to-fuel  ratio  (0/F),  the  chamber  pressure  (Pc),  and  expansion 
ratio  (s).  As  stated  in  the  introduction,  the  DEAN  will  be  powered  by  hydrogen  and 
oxygen.  Table  4  shows  the  remaining  initial  values  to  start  the  process. 


Table  4.  Chosen  Parameters 


O/F 

£ 

Pc  (MPa/psi) 

7 

125 

12/1740 

The  values  in  Table  4  are  chose  to  maximize  the  kp  of  the  design.  The  pressure 
and  expansion  ratio  were  chosen  based  on  other  like  designs  such  as  the  RE- 10.  The  exit 
Mach  number  (Me)  was  found  using  e  and  the  ratio  a  specific  heat  (y)  in  Equation  1 1 ; 
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S  - - 

M„ 


r+i 

ly-l 


(11) 


The  value  for  y  is  for  a  eombusted  flow  and  was  assumed  at  1 .2  for  the  initial 
estimates.  The  exit  Maeh  number  along  with  the  ehamber  pressure  and  y  were  used  in 
Equation  12  to  determine  the  exit  pressure: 


Based  on  the  chemieal  reaetion,  the  eharaeteristie  exhaust  veloeity  (e*)  ean  be 
found  using  Equation  13: 


ly-l 


(13) 


In  Equation  13,  the  ehamber  temperature  and  the  gas  eonstant  eome  from  the 
chemieal  kinetics  of  the  hydrogen  and  oxygen.  The  c*  efficiency  (r|c*)  is  set  to  a  realistic 
value  based  on  flow  conditions  and  empirical  evidence.  Erom  the  results  of  Equation  11- 
13,  the  theoretical  kp  can  be  found  using  Equation  14: 


In  Equation  14,  k  is  the  nozzle  efficiency  and  is  set  to  a  reasonable  number  based 
on  other  models.  The  ambient  pressure  is  set  to  zero  for  a  vacuum  providing  a  maximum 
value.  With  the  kp  and  the  thrust,  the  mass  flow  can  be  calculated  from  Equation  1 .  The 
geometry  of  the  combustion  chamber  can  then  be  calculated  as  outlined  is  section  2.2.2. 
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3,2.  Thrust  Chamber 


The  purpose  of  the  thrust  ehamber  is  to  eombine  and  eombust  the  fuel  and 
oxidizer.  Several  variables  are  involved  in  the  design  of  the  thrust  ehamber.  Below  are 
the  input  variables  that  eannot  be  ehanged  (eontrols),  variables  that  ean  be  ehanged  by  the 
designer  to  produeed  the  desired  results  (proeess  variables),  and  the  results  used  to 
determine  the  validity  of  the  design  (outputs). 

3.2.1.  Chamber  Controls 

The  eonstant  variables  in  the  design  of  the  ehamber  are  the  mass  flow  rate  of  the 
hydrogen  and  oxygen.  The  values  for  these  mass  flow  rates  eome  from  a  eombination  of 
the  requirements  of  the  overall  eyele  balanee  and  sizing. 

2.2.2.  Chamber  Process  Variables 

The  ehanging  variables  inelude  several  key  geometrie  aspeets  of  the  ehamber. 

The  primary  aspeets  are  the  radius  and  length  of  the  ehamber.  In  the  NPSS  model,  these 
variables  are  the  radius  and  volume  of  the  ehamber.  Additionally,  the  eontours  of  the 
ehamber  ean  ehange  to  get  desired  results.  These  inelude  the  outer  ehamber  wall  and  the 
inner  spike 

2.2.2.  Chamber  Outputs 

The  output  variables  used  to  evaluate  the  results  are  the  pressure  and  temperature 
of  the  ehamber.  The  pressure  and  temperature  must  be  within  material  limits  so  the 
struetural  integrity  of  the  ehamber  remains  intaet.  However,  with  higher  values,  the 
performanee  of  the  roeket  is  generally  better. 
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3,3,  Nozzle 


The  purpose  of  the  nozzle  is  to  expand  the  flow  and  produee  thrust.  Below  are 
the  eontrols  variables,  proeess  variables,  and  output  results  affeeting  the  design  of  the 
nozzle. 

3.3.1.  Nozzle  Controls 

The  variables  held  eonstant  in  the  design  of  the  nozzle  are  the  nozzle  outputs. 
These  inelude  the  pressure,  temperature,  and  mass  flow. 

3.3.2.  Nozzle  Process  Variables 

The  variables  that  ehange  are  again  the  geometry  of  the  nozzle.  For  the  nozzle, 
this  ineludes  the  throat  area,  s,  and  the  length.  Additionally,  the  operating  altitude  is  set 
by  varying  the  ambient  pressure  input.  The  throat  area  is  eonstrained  by  the  area  of  the 
ehamber  and  the  limit  of  the  Maeh  number  at  the  throat  equal  to  one.  However,  the  area 
of  the  throat  ean  be  ehanged  to  aeeommodate  both  the  ehamber  and  nozzle  but  must 
maintain  the  mass  flow. 

3.3.3.  Nozzle  Outputs 

The  overall  performanee  statisties  eharaeterizing  the  sueeess  of  this  effort  are  the 
Isp  and  the  thrust.  The  desire  is  to  have  the  Isp  and  thrust  elose  to,  if  not  more  than  the 
goal  values  listed  in  Chapter  1 . 

3,4,  Cooling  Jacket 

The  eooling  jaeket  eools  the  walls  of  the  thrust  ehamber  and  transports  the  proper 
energy  to  the  turbines.  Below  are  the  eontrol  variables,  proeess  variables,  and  output 
results  affeeting  the  design  of  the  eooling  jaeket. 
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3.4.1.  Cooling  Jacket  Controls 

The  constant  set  for  the  cooling  jacket  come  from  the  requirements  of  the  turbo¬ 
machinery.  These  include  the  mass  flow  rate  of  the  coolant,  the  inlet  temperature,  and 
the  inlet  and  outlet  pressures. 

3.4.2.  Cooling  Jacket  Process  Variables 

The  process  variables  to  be  examined  in  the  design  space  include  the  geometry  of 
the  cooling  jacket.  In  NPSS,  these  are  the  area,  volume,  and  number  of  channels. 

Another  important  parameter  includes  the  surface  area  of  the  cooling  jacket  that  transfers 
the  heat.  The  geometry  of  the  cooling  jacket  is  dependent  on  the  geometry  of  the 
chamber  and  nozzle.  However,  these  geometries,  such  as  the  cross-sectional  area  of  the 
channels,  can  be  modified  to  meet  the  requirements  of  the  cooling  jacket. 

3.4.3.  Cooling  Jacket  Outputs 

The  analysis  will  focus  on  the  wall  and  outlet  temperatures  of  the  cooling  jacket 
in  evaluating  the  performance.  The  wall  temperature  must  remain  low  enough  to  prevent 
the  wall  from  melting  while  the  outlet  temperature  must  be  high  enough  to  ensure 
sufficient  power  is  transferred  to  the  turbines. 

3.5.  NPSS 

After  defining  the  performance  parameters  of  the  engine,  the  data  was  input  into 
NPSS  to  perform  a  power  balance  and  pressure  budget.  The  power  balance  and  pressure 
budget  confirms  the  validity  of  the  DEAN  design.  Additionally,  the  NPSS  codes  has  a 
high  fidelity  resulting  in  more  accurate  compared  with  the  initial  estimates.  Additionally, 
changes  to  the  design  can  be  made  quickly  with  the  effects  of  the  changes  apparent. 
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NPSS  was  developed  by  NASA  Glenn  Researeh  Center,  in  eonjunction  with  the 
U.S.  aero  propulsion  industry  and  the  Department  of  Defense,  to  grow  teehnologies 
eapable  of  supporting  detailed  aerothermomechanical  eomputer  simulations  of  eomplete 
aircraft  engines.  NPSS  can  realistically  model  the  physical  interactions  taking  place 
throughout  an  engine,  accelerating  the  concept-to-production  development  time  and 
reducing  the  need  for  expensive  full-scale  tests  and  experiments."^^ 

A  rocket  engine  was  built  in  NPSS  using  elements  to  define  the  various 
components.  Figure  10  shows  a  schematic  layout  of  the  NPSS  model  used  to  define  the 
DEAN  that  is  based  on  the  schematic  to  the  DEAN  show  in  Eigure  1 . 


HSHAFT 


B&1  BCV1  6P2 


COi  >8 


Ow*H3 


0ml»4 


OCV2  DudiOl 


Duc<>03  CV< 


'2  OudiOa 


CVi  3  DMcfc04 


OudiOS  C 


DuekOB 


Ducfc07  OCV3  OT  0CV4  0\ 


IbI  TBCVI  T8 


TirttO  Opt  OCV1 


OSHAfT 


T«fAH  HDt  HCV1  HP1  HCV2  HP2 


HCV3  DucfcHI  CVHl  DucU42  CVH2  Du€kK3  CVH3  OudiHA  CVH4  DudiHS  CVH5  DucfcHO  CVHO  DuckH7  CVH7  Du(  ‘H8  CVH8  DuchHB  HCV4  HT  HpV5  KV 


|><^i - H 


■  ■  ■  ■  - 

1  I  I 

L 

Figure  10,  NPSS  Model  Schematic 

As  Eigure  10  illustrates,  each  block  represents  an  element  and  the  elements  are 
linked  to  one  another  through  link  ports.  NPSS  comes  with  defined  elements  for  both  air 
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breathing  and  rocket  engines.  Additionally,  new  elements  can  be  programmed  by  the 
user.  Each  element  has  a  set  of  input/output  variables,  independent/dependent  variables, 
and  equations.  The  user  defines  the  input  variables,  though  not  all  inputs  need  to  be 
defined.  Additionally,  in  the  model,  the  composition  of  the  flow,  such  as  oxygen,  and  the 
thermodynamic  data  pack  are  defined.  The  inputs  defined  by  the  user,  and  the  inputs  of 
the  composition  and  thermodynamic  data,  are  used  by  the  equations,  with  changing  of  the 
independent  and  dependent  variables  as  defined  in  the  element,  to  determine  the  outputs. 
The  elements  are  linked  through  linkPorts.  There  are  link  ports  for  fluid  and  heat  flow. 
Each  element  defines  which  type,  and  how  many  ports  the  element  has. 

The  full  model  used  for  the  DEAN  is  presented  in  Appendix  B.  Eor  the 
combustion  chamber,  labeled  COMB  in  Eigure  10,  was  modeled  using  the  RocketCombl 
element.  This  element  requires  inputs  of  chamber  radius,  chamber  volume,  and  the 
composition  of  the  oxidizer  and  fuel  inlet,  and  outlet.  Additionally,  guesses  for  the 
oxidizer-fuel  ratio  (OER),  chamber  temperature,  chamber  pressure,  and  weight  flow  are 
required.  The  aerospike  nozzle,  labeled  NOZZ  in  Eigure  10,  is  modeled  with  the  element 
RocketNozzle  and  requires  inputs  of  the  throat  area,  area  ratio  (s),  and  ambient  pressure. 

The  cooling  jacket  model  requires  the  use  of  three  elements;  one  for  a  pressure 
drop,  one  for  change  in  energy,  and  one  for  the  transfer  of  thermal  energy,  because  of  the 
way  NPSS  works.  Eigure  1 1  shows  one  of  the  cooling  jackets  of  the  NPSS  model. 
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Figure  11,  NPSS  Cooling  Jacket 


The  Valve04  element  modeled  the  pressure  drop  of  the  cooling  fluid.  The 
elements  DuckH2  in  Figure  11  is  an  example  of  a  Valve04  element  in  the  DEAN  model. 
The  Valve04  element  requires  an  input  of  the  cross-sectional  area  and  friction  factor. 
Furthermore,  a  guess  at  the  weight  flow  is  required.  The  CoolingVolume  element  was 
used  to  model  the  change  in  energy  of  the  flow.  The  elements  CVHl  and  CVH2  in  Figure 
1 1  are  examples  of  this  element  in  the  DEAN  model.  The  CoolingVolume  element 
requires  inputs  of  the  volume,  area,  and  number  of  channels.  In  addition,  a  guess  at  the 
total  pressure  and  total  specifies  enthalpy.  The  Wall2  element  was  used  to  model  the  heat 
transfer  from  the  chamber  to  the  cooling  fluid.  The  elements  Hwalll  and  Hwall?  in 
Figure  1 1  are  examples  of  this  element  in  the  DEAN  model.  This  element  requires  and 
input  of  the  specific  heat  (Cp),  thermal  coefficient  (k),  and  the  density.  Additionally,  a 
guess  of  the  wall  temperature  is  required. 

Other  elements  used  in  the  DEAN  model  are  the  Starter,  Pump,  Turb02,  and  Shaft 
elements.  In  Figure  10,  OTNK  is  an  example  of  a  Starter  element.  These  elements  begin 
and  end  mass  flow  respectively.  These  elements  require  a  temperature  and  pressure  as 
inputs.  The  Pump  element  models  a  fluid  pump  and  is  used  for  such  elements  as  OP  in 
Figure  10.  This  element  requires  and  input  of  weight  flow,  gear  ratio,  efficiency,  and 
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pressure  rise.  To  model  a  turbine  the  Turb02  element  was  used.  One  sueh  element 
shown  in  Figure  10  is  HT.  This  element  requires  an  input  of  weight  flow,  effieieney,  and 
pressure  ratio.  The  Pump  and  Turb02  element  are  linked  by  the  Shaft  element  whieh 
balanees  the  power.  In  Figure  10,  HSHAFT  is  a  Shaft  element.  This  element  only 
requires  and  input  of  the  revolutions  per  minute  (rpm). 

To  allow  mass  flow  from  one  to  another,  elements  are  linked  with 
UnReactedFluidPorts .  These  ports  do  not  require  any  speeifie  imports,  however  the  user 
may  define  a  parameter.  Heat  is  passed  between  elements  using  ThermalPorts.  Inputs 
for  these  ports  inelude  a  radius  of  eurvature,  eross-seetional  flow  area,  and  a  surfaee  area. 

3.6.  TDK 

The  Two-Dimensional  Kineties  (TDK)  eode  was  used  to  establish  the  eontour  of 
the  aerospike  nozzle.  The  results  of  the  this  eode  will  aid  in  the  inputs  for  the  nozzle  in 
NPSS  and  help  eonfirm  the  results  obtained. 

The  TDK  eomputer  program  represents  the  eulmination  of  work  done  by  many 
people  over  the  last  35  plus  years.  A  series  of  eontraets  funded  through  NASA  and  the 
ICRPG  (now  (JANNAF)  resulted  in  TDK  beeoming  the  JANNAF  Standard  Code  for 
predieting  the  nozzle  performanee  in  liquid  propellant  roeket  engines  in  1967."^"^ 
Currently,  TDK  is  property  of  and  under  development  by  Software  and  Engineering 
Assoeiates,  Ine.  (SEA). 

The  TDK  program  is  eomposed  of  seven  modules,  ODE,  SCAP,  ODK,  TRANS, 
MOC,  BEM,  and  MABE.  These  modules  allow  for  a  eomplete  two  dimensional  non- 
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equilibrium  nozzle  performance  calculation  with  a  boundary  layer.  Each  module  requires 
one  set  of  inputs  and  can  be  used  alone,  in  part,  or  in  combination  with  the  other  modules. 

3,7,  Sensitivity  Analysis 

A  sensitivity  analysis  on  the  cooling  jacket  was  performed  to  focus  the  design 
changes  in  the  cooling  jacket  parameters.  To  accomplish  this,  a  simplified  NPSS  model 
was  created  with  a  single  wall  circuit  for  both  hydrogen  and  oxygen  flow.  Figure  12 
show  the  schematic  of  the  sensitivity  model. 
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Figure  12,  Sensitivity  Model 


The  primary  outputs  that  would  be  used  to  judge  the  merit  of  the  design  are  the 
final  temperature  of  the  flow  and  the  temperature  of  the  wall.  The  final  temperature  must 
be  great  enough  to  allow  the  turbine  to  produce  the  required  amount  of  horsepower  to 
drive  the  pumps.  The  wall  temperature  must  be  sufficiently  low  to  insure  structural 
stability  of  the  material.  The  parameters  that  were  controlled  are  the  mass  flow,  change 
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in  pressure,  eross-seetional  area,  and  surface  area.  Equation  15  was  used  to  compute  the 
sensitivity  of  each  parameter; 


Sensitivity  = 


AT 


AS/ 

7q 


(15) 


In  Equation  15,  Q  represents  the  parameter  that  was  changed  and  T  is  the 
temperature.  Table  5  and  Table  6  show  the  initial  values  and  the  results  of  the  sensitivity 
analysis. 


Table  5,  Hydrogen  Sensitivity  Analysis 


Hydrogen 

Original 

New 

Tout 

initial 

Tout 

new 

Twall 

initial 

Twall 

new 

Tout 

sensitivity 

Twall 

sensitivity 

mdot 

10.0 

9.06 

221.6 

223.7 

729.7 

766.2 

-22.1 

-382.5 

AP 

500 

700 

221.6 

221.2 

729.7 

730.2 

-1.18 

1.25 

Cross 

Area 

1 

1.5 

221.6 

220.5 

729.7 

916.9 

-2.28 

374.4 

Surface 

Area 

50 

55 

221.6 

223.9 

729.7 

729.1 

22.4 

-6.51 

Table  6,  Oxygen  Sensitivity  Analysis 


Oxygen 

Original 

New 

Tout 

initial 

Tout 

new 

Twall 

initial 

Twall 

new 

Tout 

sensitivity 

Twall 

sensitivity 

mdot 

10.0 

8.44 

281.1 

285.1 

1656 

2066.1 

-25.9 

-2629.4 

AP 

500 

700 

281.1 

281.2 

1656 

1665.5 

0.313 

23.75 

Cross 

Area 

1 

1.5 

281.1 

280.4 

1656 

1780.6 

-1.38 

249.14 

Surface 

Area 

50 

55 

281.1 

282.8 

1656 

1855.6 

17.7 

1995.6 

The  initial  values  in  Table  5  and  Table  6  were  chosen  as  realistic  round  numbers 
for  ease  of  analysis.  The  highlighted  boxes  are  those  that  affect  the  temperature  of  the 
fluid  and  wall  the  most.  Erom  Table  5,  a  change  in  mass  flow  results  in  a  proportionally 


62 


large  change,  in  the  opposite  direction,  in  both  the  flow  and  wall  temperature. 
Additionally,  a  change  in  cross-sectional  area  results  in  a  substantial  change  in  wall 
temperature  in  the  same  direction.  A  change  in  surface  area  results  in  a  considerable 
change  in  the  flow  temperature  in  the  same  direction. 

From  Table  6,  a  change  in  oxygen  mass  flow  has  a  similar  effect  on  the  flow  and 
wall  temperature  as  the  hydrogen  mass  flow.  Furthermore,  a  change  in  the  surface  area 
has  a  significant  change  in  both  flow  and  wall  temperature  in  the  same  direction. 

Plainly,  the  mass  flow,  cross-sectional  and  surface  area  have  the  largest  effect  on 
the  final  temperature  of  the  flow  and  the  wall  temperature. 

3,8,  Beginning  Frame  Work 

The  inputs  to  the  NPSS  model  require  many  geometric  inputs.  Therefore,  a 
general  frame  for  the  chamber  was  developed  to  begin  the  process  of  defining  the 
geometry  inputs  for  NPSS.  Figure  13  shows  the  chamber  framework. 


Figure  13,  Initial  Contour 


The  contour  illustrated  in  Figure  13  helped  to  visualize  the  design  of  the  chamber. 
With  the  framework  of  the  chamber  and  nozzle;  the  geometry,  such  as  the  radius  and 
length,  could  be  adjusted  based  on  the  needs  of  the  design.  The  changes  would  then  be 
used  to  determine  inputs  into  NPSS  such  as  the  surface  area  for  a  ThermalPort. 
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3.8.1.  Cooling  Channel  Design 

Another  aspect  affected  by  the  contour  of  the  chamber  is  the  design  of  the  cooling 
channels.  The  number  and  size  of  the  cooling  channels  is  constrained  by  the  radius  of  the 
chamber  at  each  point.  Initially  both  milled  chaimel  and  tubular  channels  were 
considered  but  the  milled  channels  proved  to  be  superior  in  performance  and 
manufacturability.  Figure  14  shows  the  cross-section  used  to  define  the  cooling  chaimel 
geometry. 


centerline 


Figure  14,  Cooling  Channel  Cross-Section 


The  channel  depicted  in  Figure  14  represents  half  of  the  actual  channel  where  it  is 
mirrored  about  the  centerline.  The  shaded  region  represents  the  solid  material  while  the 
clear  area  is  where  the  fluid  flows.  The  parameter  ’a’  represents  the  half  spacing  between 
channels,  ‘w’  is  the  half  width,  ‘t’  is  the  thickness  between  the  chamber  wall  and  the 
channel,  ‘h’  is  the  height  of  the  channel.  To  determine  the  number  of  channels  that  will 
fit,  the  circumference  of  the  chamber  is  divided  by  twice  the  sum  of  ‘a’  and  ‘w’,  since 
one  channel  is  represented  by  two  times  their  lengths. 
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Changing  the  thickness  directly  influences  the  hot  and  cold  side  wall-temperature. 
The  ratio  between  ‘a’  and  ‘w’  determines  the  total  surface  area  of  the  channels  thereby 
influencing  the  wall  temperature  and  the  temperature  of  the  fluid  as  indicated  in  the 
sensitivity  analysis.  Changing  the  height  of  the  chamber  affects  the  cross-sectional  area 
of  the  channel  thus  affecting  the  wall  and  fluid  temperature  as  well. 

3.8.2.  Material  Choice 

The  methodology  for  choosing  the  materials  was  based  on  the  properties 
presented  here.  The  material  of  the  chamber  must  deal  with  extreme  thermal  and 
pressure  forces.  To  handle  the  heat,  materials  with  high  thermal  conductivity  and 
elevated  melting  points  were  considered.  These  include  metals  like  copper  and  nickel  as 
well  as  composites  like  silicon  carbide.  Table  7  shows  a  list  of  materials  and  their 
important  properties. 


Table  7.  Material  Thermal  Properties^’ 


Material 

Melting 

Point 

(R) 

50%  of 
Melting 
(R) 

Thermal 

Conductivity 

(Btu/in*R*s) 

Copper 

2444.4 

1222.2 

0.005377 

Beryllium 

2790 

1395 

0.002682 

Chromium 

3812.4 

1906.2 

0.001257 

Cobalt 

3184.2 

1592.1 

0.00133 

Iridium 

4896 

2448 

0.001971 

Molybdenum 

5209.2 

2604.6 

0.001851 

Nickel 

3110.4 

1555.2 

0.001216 

Rhodium 

4024.8 

2012.4 

0.002012 

Silicon 

3033 

1516.5 

0.001985 

Silver 

2223 

1111.5 

0.005753 

Tungsten 

6588 

3294 

0.002333 

Niobium 

4933 

2045 

0.00072 

Silicon  Carbide 

5580 

2790 

0.006571 

Carbon/Carbon 

4091.4 

5045.7 

0.02615 
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The  50%  melting  point  temperature  in  Table  7  will  be  the  values  the  wall 
temperature  is  eompared  to  in  order  to  determine  if  the  material  is  viable  given  that  the 
strength  of  the  material  should  remain  unchanged  up  to  that  point. 

Surrounding  the  cooling  jacket  will  be  a  structural  jacket  to  withstand  the 
pressure.  Materials  for  this  task  will  need  high  tensile  strength  to  handle  the  pressure  and 
low  density  to  minimize  the  weight.  Table  8  displays  materials  considered  for  structural 
support. 


Table  8,  Structural  Materials 


Material 

Ftu  (ksi) 

Density 

(slugs/ft^) 

A1 7075-T6 

80 

5.4 

Copper 

33 

17 

High-strength  Steel 

180 

15.2 

Nickel 

110 

17 

Titanium 

170 

8.7 

Tungsten 

600 

37 

Carbon  Fiber 

30.5 

2.72 

The  values  in  Table  8  are  average  values  and  for  materials  like  carbon  fiber,  are 
highly  dependent  on  the  specific  usages  and  manufacture.  Another  consideration  for  the 
material  choice  is  how  it  reacts  with  oxygen.  According  to  a  NASA  report, materials 
used  in  an  oxygen  environment  should  have  a  low  heat  of  combustion  and  a 
comparatively  low  burn  factor.  Equation  16  shows  how  burn  factor  can  be  defined; 

A//“ 

BF  = - ^  (16) 

a 

In  Equation  16,  AH^’f  is  the  heat  of  oxidation  and  a  is  the  thermal  diffusivity. 

Table  9  shows  to  compare  material  compatibility  with  oxygen. 


66 


Table  9,  Oxygen  Compatibility 


Material 

Heat  of 

Combustion'*^ 

(cal/g) 

Heat  of 

Formation'*^ 

(kcal/gmole) 

Thermal 

Diffusivity*^ 

(m^/s) 

Burn 

Factor 

Copper 

585 

41.8 

117 

0.357 

Aluminum 

7425 

400.4 

97.1 

4.12 

Cobalt 

970 

57.1 

26.6 

2.15 

Molybdenum 

1458 

182.65 

53.7 

3.40 

Nickel 

980 

57.3 

23 

2.49 

Tungsten 

1093 

140.94 

68.3 

2.06 

Titanium 

4710 

225.5 

9.32 

24.19 

Aluminum  and  Titanium  are  metals  known  to  be  highly  reaetive  with  oxygen. 
Copper  and  Niekel  are  generally  unreaetive  with  oxygen.  The  values  in  Table  9  for  these 
materials  are  used  as  baselines  to  evaluate  other  materials.  For  example,  Tungsten’s  burn 
faetor  is  elose  to  copper,  therefore  it  would  be  considered  compatible  with  oxygen.  On 
the  other  hand.  Molybdenum  has  a  burn  factor  closer  to  aluminum  and  consequently 
would  be  considered  to  be  reactive  with  oxygen.  On  the  composite  side,  Schoenman"^^ 
states  that  silicon  carbine  is  compatible  with  use  in  oxygen  environments. 
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4,  Analysis  and  Results 


4.1,  Chapter  Overview 

This  Chapter  presents  the  results  obtained  from  NPSS  and  TDK  for  the  DEAN 
deseribed  in  the  methodology  of  Chapter  three.  The  NPSS  model  of  the  DEAN  eonsists 
of  five  major  components;  the  oxygen  turbo  machinery,  the  oxygen  cooling  jacket,  the 
hydrogen  turbo  machinery,  the  hydrogen  cooling  jacket,  and  the  chamber/nozzle.  The 
turbo  machinery  is  analyzed  in  related  documents  by  Capt  Michael  Arguello'^^  for  the 
hydrogen  side  and  Capt  William  Strain'^^  for  the  oxygen  side.  Presented  below  are  the 
results  pertaining  to  the  cooling  jackets  and  the  chamber/nozzle.  The  code  that  produced 
these  results  is  presented  in  Appendix  C. 

4.2,  Results  of  Simulation  Scenarios 

The  chamber  and  nozzle  contour  was  defined  to  determine  the  geometric  inputs 
for  NPSS  and  visualize  the  deign  of  the  DEAN.  Eigure  15  shows  the  final  contour  of  the 
chamber/nozzle  used  in  the  NPSS  model. 
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Figure  15,  Chamber/Nozzle  Contour  (all  dimension  in  inches) 
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Figure  15  shows  the  ehamber/nozzle  split  into  several  segments.  Splitting  the 
eontour  into  segments  allowed  for  a  profile  of  the  wall  and  fluid  temperature  along  the 
length  of  the  contour.  The  chamber  is  split  into  five  equal  length  sections  with  an 
additional  section  at  the  throat  since  that  point  usually  experiences  the  most  extreme 
temperatures.  The  nozzle  is  split  into  three  sections  with  the  two  equal  length  sections 
being  cooled  and  the  last  section  not  cooled.  The  mid-points  shown  in  Figure  15 
represent  the  points  where  the  geometry  for  the  elements  are  defined  for  each  section. 

4.2.1.  Hydrogen  Cooling  Jacket  Results 

The  hydrogen  cooling-jacket  performs  the  task  of  cooling  the  inner  chamber  wall 
and  part  of  the  nozzle  while  providing  energy  to  the  hydrogen  turbine  to  power  the 
pumps.  The  hydrogen  flow  is  split  50/50  after  the  first  pump  as  shown  in  Figure  10.  The 
result  is  a  mass  flow  of  7.55  Ibm/s  (3.42  kg/s)  for  cooling.  The  cooling  jacket  consists  of 
the  elements  between  the  second  pump  and  turbine.  The  DuckH#  elements  are  Valve04 
elements  and  are  responsible  for  the  pressure  drop  in  the  jacket.  The  pressure  drops  from 
4050  psi  (27.9  MPa)  at  the  exit  of  the  second  pump  to  3625  psi  (25  MPa)  at  the  inlet  of 
the  turbine.  The  friction  factor  (K)  of  the  elements  is  set  to  allow  for  a  constant  pressure 
drop  through  the  jacket. 

The  CoolingVolume  elements  attracted  to  the  Wall2  elements  allow  for  the  heat 
transfer  from  the  chamber  to  the  fluid.  These  elements  consist  of  CVH#  and  Hwall# 
elements.  There  are  a  total  of  314  channels  with  the  dimensions  of  the  channels  changing 
for  each  element  as  revealed  in  Table  10. 
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Table  10,  Hydrogen  Channel  Dimensions 


Element 

r(in) 

a  (in) 

w  (in) 

li(in) 

A(in^) 

AR 

CVHl 

2.5 

0.01 

0.015 

0.02 

0.19 

1.33 

CVH2 

4 

0.01 

0.03 

0.02 

0.38 

0.67 

CVH3 

4.5 

0.01 

0.035 

0.02 

0.44 

0.57 

CVH4 

4 

0.01 

0.03 

0.02 

0.38 

0.67 

CVH5 

2.75 

0.01 

0.0175 

0.03 

0.33 

1.71 

CVH6 

2.25 

0.01 

0.0125 

0.04 

0.31 

3.2 

CVH7 

2 

0.01 

0.01 

0.05 

0.31 

5 

CVH8 

2 

0.01 

0.01 

0.05 

0.31 

5 

The  variables  in  Table  10  eorrespond  to  Figure  14,  r  is  the  radius  and  A  is  the 
total  eross-seetional  area  of  the  eooling  flow.  The  ehange  in  dimension  in  Table  10 
represents  the  ehange  in  priority  of  the  ehannel  along  the  eontour  of  the  ehamber.  At  the 
throat,  CVH3,  the  AR  is  lowest  to  allow  for  maximum  heat  transfer  and  minimum  wall 
temperature.  At  the  end  of  the  flow,  CVH8,  the  eoneern  is  preventing  the  flow  from 
going  supersonie.  A  shook  wave  in  the  ehannel  would  greatly  inorease  the  pressure  loss 
and  may  also  oause  damage  to  the  struoture.  Table  1 1  shows  the  Maoh  number  for  eaoh 
element. 


Table  11,  Hydrogen  Mach  Numbers 


Element 

P  (psi) 

T(R) 

V  (ft/s) 

a  (ft/s) 

Mach  # 

CVHl 

4002 

145.3 

1630 

4556.4 

0.35 

CVH2 

3954 

270.6 

1336 

4338.5 

0.30 

CVH3 

3907 

276.6 

1185 

4333 

0.27 

CVH4 

3860 

415.7 

1972 

4689.1 

0.42 

CVH5 

3813 

477.2 

2584 

4870.7 

0.53 

CVH6 

3766 

524.9 

3029 

5012.1 

0.60 

CVH7 

3719 

567.4 

3278 

5137.1 

0.63 

CVH8 

3672 

609.7 

3532 

5260.7 

0.67 

In  Table  11,  the  pressure,  temperature  and  velooity  are  the  exit  properties  of  eaoh 
element.  The  values  for  speed  of  sound  (a)  oome  from  the  National  Institute  of  Standards 
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and  Technology  (NIST)"^^  database  for  hydrogen  for  the  specific  pressure  and 
temperature.  The  greatest  Mach  number  is  0.67  at  the  end  of  the  cooling  jacket  and  is  not 
high  enough  to  cause  any  issues  associated  with  shock  waves.  The  Mach  number  is  not  a 
issue  for  hydrogen  till  around  0.9,  therefore  there  is  some  leeway  were  this  is  concerned. 

The  heat  flow  and  heating  coefficients  used  to  calculate  the  wall  temperature  are 
outputs  from  the  NPSS  model.  However,  the  wall  temperature  NPSS  calculates  is  only  a 
balance  of  heat  in  and  out  of  the  wall  and  does  not  take  into  account  material  properties. 
To  calculate  the  wall  temperature  for  different  materials  another  calculation  is  required. 
The  calculation  balances  the  three  heat  equations  shown  in  Equation  17  originating  from 
Equations  7  and  Equation  8: 

q  =  h„A„(T„  -T,„)^kA„[ ~ k  /V -ic (T,k  - Tc )  (17) 

V  ' 

In  Equation  17,  q,  hn,  he.  Ah,  Ac,  Th,  and  Tc  come  from  the  NPSS  model.  The 
value  for  k  is  dependent  of  the  material  chosen,  and  t  is  set  to  0.02  in  (0.508  mm)  for  the 
design  of  the  cooling  jacket.  The  results  for  various  materials  are  presented  in  Appendix 
E.  Eigure  16  shows  the  graph  of  the  wall  temperature  for  copper  along  the  axial 
coordinate. 
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In  Figure  16,  the  injector  faceplate  is  at  zero  and  the  beginning  of  the  cooling 
circuit  in  the  nozzle  is  at  31  in  (0.787  m)  (the  hydrogen  travels  from  right  to  left).  The 
throat  is  at  24  in  (0.61  m).  As  expected  the  hottest  wall  temperature  is  at  the  throat  at 
1224  R  (680  K).  The  50%  melting  point  for  copper  is  1222  R  (678  K).  The  calculated 
temperature  of  the  wall  is  within  a  percent  of  the  50%  melting  point  of  copper,  hence  in 
the  current  configuration  the  wall  temperature  maintains  a  reasonable  level  during  normal 
operation. 

The  final  requirement  for  the  cooling  jacket  to  meet  is  the  temperature  of  the  flow 
entering  the  turbine.  Table  12  shows  the  outlet  temperature  of  each  element. 
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Table  12.  Hydrogen  Flow  Temperature 


Element 

T(R) 

CVHl 

145.3 

CVH2 

270.6 

CVH3 

276.6 

CVH4 

415.7 

CVH5 

477.2 

CVH6 

524.9 

CVH7 

567.4 

CVH8 

609.7 

Table  12  shows  the  temperature  into  the  eooling  jaeket  is  145.3  R  (80.7  K)  at  the 
subsequent  rise  in  temperature  to  The  result  is  a  temperature  raise  of  464.3  R  (257.9  K). 
The  temperature  into  the  turbine  is  609.7  R  (338.7  K)  that  allowed  the  turbine  to  produee 
3573  hp,  adequate  to  power  the  hydrogen  pumps. 

4. 2. 2.  Oxygen  Cooling  Jacket  Results 

The  oxygen  eooling-jaeket  is  responsible  for  eooling  the  outer  ehamber  wall  and 
providing  energy  to  the  oxygen  turbine  to  power  the  pump.  Unlike  the  hydrogen,  100% 
of  the  oxygen  flows  through  the  jaeket  as  illustrated  in  Figure  10.  The  result  is  a  mass 
flow  of  106  Ibm/s  (48.08  kg/s)  for  eooling.  The  eooling  jaeket  eonsists  of  the  elements 
between  the  pump  and  turbine.  The  DuckO#  elements  perform  the  same  task  of  pressure 
drop  like  the  DuckH#  elements.  The  pressure  drops  from  4635  psi  (3 1 .9  MPa)  at  the  exit 
of  the  pump  to  3666  psi  (25.2  MPa)  at  the  inlet  of  the  turbine. 

The  elements  of  CVO#  and  Owall#  elements  perform  the  same  as  CVH#  and 
Hwall#  elements.  The  oxygen  eooling-jaeket  has  785  ehannels  with  the  dimensions  of 
the  ehannels  ehanging  for  eaeh  element  as  revealed  in  Table  13. 
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Table  13,  Oxygen  Channel  Dimensions 


Element 

r(in) 

a  (in) 

w  (in) 

li(in) 

A  (in2) 

AR 

CVOl 

5 

0.01 

0.01 

0.03 

0.47 

3 

CV02 

5.5 

0.01 

0.012 

0.04 

0.75 

3.33 

CV03 

6 

0.01 

0.014 

0.04 

0.88 

2.86 

CV04 

6 

0.01 

0.014 

0.04 

0.88 

2.86 

CV05 

6 

0.01 

0.014 

0.04 

0.88 

2.86 

CV06 

6 

0.01 

0.014 

0.05 

1.1 

3.57 

The  variables  in  Table  13  are  the  same  as  Table  10.  The  ehange  in  dimension  in 
Table  13  are  due  to  the  same  faetors  that  influenee  the  ehange  in  AR  seen  in  Table  10. 
Table  14  shows  the  Maeh  number  for  eaeh  element. 


Table  14,  Oxygen  Mach  Numbers 


Element 

P  (psi) 

T(R) 

V  (ft/s) 

a  (ft/s) 

Mach  # 

CVOl 

4497 

179.3 

448.2 

3186.3 

0.14 

CV02 

4359 

303.9 

379.8 

1980.1 

0.19 

CV03 

4221 

379.6 

428.1 

1497.1 

0.28 

CV04 

4084 

452.2 

576.2 

1326.1 

0.43 

CV05 

3947 

531.4 

758.3 

1321.6 

0.57 

CV06 

3809 

616.8 

764.3 

1354.6 

0.56 

In  Table  14,  the  pressure,  temperature  and  veloeity  are  the  exit  properties  of  eaeh 
element.  The  values  for  speed  of  sound  are  for  oxygen  at  the  speeified  pressure  and 
temperature.  The  greatest  Maeh  number  is  0.57  at  CV05  and  is  not  high  enough  to  eause 
any  issues  as  stated  earlier  with  the  hydrogen  side.  The  Maeh  number  beeomes  an  issue 
at  around  0.6  for  oxygen.  Therefore,  this  design  is  at  the  high  end  of  what  is  expeetable. 

The  wall  temperature  for  the  oxygen  side  is  eomputed  using  Equation  17  similar 
to  the  hydrogen  side.  The  results  for  various  materials  are  presented  in  Appendix  E. 
Eigure  17  shows  the  graph  of  the  wall  temperature  for  silieon  earbide  along  the  axial 
eoordinate. 
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In  Figure  17  the  injector  faceplate  is  at  0  inches  and  the  beginning  of  the  cooling 
circuit  in  the  throat  is  at  24  inches  (0.61  m)  (the  oxygen  travels  from  right  to  left). 

Silicon  carbide  is  an  expensive  and  exotic  material,  however,  lower  cost  alternatives  we 
unable  to  meet  the  performance  requirements.  Additionally,  silicon  carbide  is  know  to 
works  well  with  oxygen,  that  was  a  consideration  outlined  in  section  3.8.2.  Silicon 
carbide  has  been  successfully  tested  with  an  actively  cooled  configuration.^”  Therefore,  it 
is  manufacturable  although  it  may  be  expensive.  Another  material  that  might  work  is 
tungsten.  According  to  the  bum  factor  presented  in  Chapter  three,  tungsten  should  be 
compatible  with  oxygen  but  the  manufacture  may  present  problems. 

As  expected  the  hottest  wall  temperature  is  at  the  throat  at  1625  R  (902  K).  The 
50%  melting  point  for  silicon  carbide  is  2790  R  (1550  K).  The  calculated  temperature  of 
the  wall  is  34  %  lower  than  the  50%  melting  point  of  silicon  carbide;  hence,  in  the  current 
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configuration  the  wall  temperature  maintains  a  reasonable  level  during  normal  operation 
and  allows  room  for  growth. 

The  final  requirement  for  the  eooling  jaeket  to  meet  is  the  temperature  of  the  flow 
entering  the  turbine.  Table  15  shows  the  outlet  temperature  of  eaeh  element. 


Table  15.  Oxygen  Flow  Temperature 


Element 

T(R) 

CVOl 

179.3 

CV02 

303.9 

CV03 

379.6 

CV04 

452.2 

CV05 

531.4 

CV06 

616.8 

Table  15  reveals  the  temperature  into  the  eooling  jaeket  is  179.3  R  (99.6  K)  and 
inereases  along  the  jaeket  to  a  final  temperature  of  616.8  R  (342.6  K).  The  result  is  a 
temperature  rise  of  437.4  R  (243  K)  in  the  jaeket.  The  temperature  out  of  the  eooling 
jaeket  allowed  the  turbine  to  produee  2587  hp,  suffieient  to  power  the  oxygen  pump. 

4.2.3.  Chamber/Nozzle  Results 

The  hydrogen  and  oxygen  flow  into  the  ehamber  and  eombust,  and  then  expand 
through  the  nozzle.  The  elements  of  interest  are  the  COMB  and  NOZZ  elements  as 
depleted  in  Figure  10.  The  COMB  element  has  an  input  of  C_02_H2  as  the  eomposition 
to  designate  that  the  reaetion  is  of  hydrogen  and  oxygen  and  the  radius  of  the  ehamber  is 
6  inehes.  The  NOZZ  element  has  inputs  of  15.9  in  for  the  throat  area,  125  for  the 
expansion  ratio,  and  0.1  psi  for  the  ambient  pressure.  Both  elements  have  thermal  ports 
that  are  eonneeted  to  Wall2  elements  for  heat  transfer.  The  inputs  for  radius  of  eurvature. 
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flow  area,  and  surface  area  (areaHx  in  the  code)  come  from  Figure  15.  The  results  for 
these  elements  are  presented  in  Table  16. 


Table  16,  Chamber  Performance 


O/F 

7.03 

Pc 

1739.1 

psi 

12 

MPa 

Te 

6586.3 

R 

3659 

K 

F. 

57231.9 

Ibf 

254 

kN 

Total  mass  flow 

121.1 

Ibm/s 

54.9 

k^/s 

Hydrogen  mass  flow 

15.1 

Ibm/s 

6.8 

kg/s 

Oxygen  mass  flow 

106 

Ibm/s 

48.1 

kg/s 

Isp 

472.7 

sec 

Table  16  show  this  model  meets  the  thrust  and  hp  goals  set  forth  in  Chapter  1. 
Specifically,  the  thrust  goal  is  50,000  Ibf  (222.4  kN),  while  the  result  from  the  DEAN 
NPSS  model  is  57,231  Ibf  (254  kN).  This  is  a  14%  increase  in  thrust  over  the  goal.  The 
Isp  goal  is  464  sec  whereas  the  DEAN  Isp  from  the  NPSS  model  is  472.7  sec.  This  is  a 
1 .9%  increase  in  fip  over  the  goal.  The  increase  of  the  thrust  and  kp  over  the  goal  values 
shows  the  model  performs  better  than  expected.  Additionally,  the  values  for  the  chamber 
are  close  to  the  initial  parameters  shown  in  Table  4.  The  O/E  is  7.03  in  the  model  while 
the  initial  value  was  seven.  The  chamber  pressure  of  1739  psi  (12  MPa)  is  within  one  psi 
of  the  initial  selection.  This  demonstrates  the  initial  assumptions  were  sound  and  the 
NPSS  chamber  model  performs  close  to  what  was  projected. 

4.2.4.  TDK  Nozzle 

TDK  04  determined  the  contour  of  the  nozzle.  The  code  developed  for  the  DEAN 
is  given  in  Appendix  C.  In  developing  the  code,  the  goal  was  to  match  many  of  the 
chamber  properties,  such  as  pressure  and  temperature,  to  the  NPSS  model.  TDK  then 
generated  a  nozzle  based  on  these  properties  and  non-dimensionalized  by  the  throat 
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radius.  The  contour  is  then  incorporated  into  the  design.  The  code  developed  results  in 
the  following  parameters  in  Table  17. 


Table  17,  TDK  Properties 


Pc 

1739 

psi 

12 

Mpa 

Tc 

6586 

R 

3659 

K 

O/F 

7 

Mdot 

121.1 

Ibm/s 

54.9 

k^/s 

Fg 

66167 

Ibf 

294 

kN 

IsD 

546 

s 

Table  17  illustrates  the  chamber  properties  for  NPSS  and  TDK  are  the  same. 
However,  the  hp  and  thrust  calculated  by  TDK  is  much  larger  than  the  NPSS  model. 
Some  difference  is  due  to  the  differences  in  methods  of  TDK  and  NPSS.  Nevertheless, 
the  TDK  hp  is  still  above  the  theoretical  limit  of  500  sec  for  an  O2/H2  engine.  This  is 
explained  by  the  fact  that  the  presented  hp  is  that  of  the  entire  control  volume  TDK  used 
to  make  calculations  and  therefore  is  inflated.  Still,  the  definition  of  the  control  volume 
in  unclear  in  the  output  so  real  fip  remains  unknown.  The  nozzle  contour  was  still  used 
for  NPSS  since  the  mass  flows  matched.  Figure  18  shows  the  non-dimensional  contour 
with  a  temperature  profile. 


Figure  18,  Non-dimensional  Nozzle  Contour  with  Temperature  Profile 
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Figure  18  reveals  the  length  of  the  nozzle  is  about  three  times  the  radius  at  the 
throat.  Additionally,  Figure  18  shows  the  change  in  temperature  along  the  nozzle  contour 
and  is  the  reason  for  part  of  the  nozzle  being  actively  cooled  due  to  the  decrease  in  the 
temperature  of  the  flow. 

Truncation  is  ending  the  nozzle  before  the  flow  is  fully  expanded.  Figure  19 
show  a  truncated  nozzle. 


In  Figure  19,  the  red  line  denotes  a  length  of  75%  of  a  fully  contoured  nozzle,  the 
green  line  50%  and  the  black  like  25%.  To  determine  the  affect  of  truncation,  a  model 
where  the  ratio  of  length  to  throat  radius  varies  was  created.  Table  18  shows  the  relation 
of  nozzle  length  to  kp  and  mass. 


Table  18,  Truncation  of  Nozzle 


L 

Isp  (s) 

AIsp  (s) 

Mass  (Ibm) 

AMass 

100% 

546 

1.05 

75% 

548 

2 

0.96 

8.5% 

50% 

545 

-1 

0.65 

38% 

25% 

537 

-9 

0.38 

64% 
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The  mass  in  Table  18  assumes  niobium,  0.3096  Ibm/in^  (8570kg/m^),  is  used  for 
the  entire  nozzle  with  a  throat  radius  of  one  inch.  Niobium  (columbium)  is  chosen  since 
it  is  the  material  for  the  un-cooled  section  of  the  nozzle  and  results  in  quicker 
calculations.  Table  18  indicated  truncation  of  the  nozzle  by  75%  results  in  a  0.4% 
increase  in  kp  and  a  8.5%  decrease  in  mass.  Why  the  kp  increased  may  be  due  to  the 
same  factors  causing  the  kp  to  be  inflated  over  the  theoretical  maximum.  Further 
reduction  in  length  behaved  as  expected  with  a  0.2%  and  1.6%  decrease  in  Isp  with  a 
38%  and  64%  decrease  in  weight  for  50%  and  25%  of  the  total  length  respectively.  The 
results  indicate  truncation  of  the  nozzle  has  little  impact  on  kp-  However,  as  indicated  in 
Chapter  1,  even  a  small  change  in  kp  has  a  major  impact  on  the  capabilities  of  the  launch 
vehicle.  Therefore,  at  this  point  the  design  will  utilize  the  full-length  nozzle.  Figure  20 
illustrates  the  full-length  nozzle  contour  for  the  DEAN  design. 


The  contour  in  Figure  20  comes  from  the  contour  of  the  TDK  code  presented  in 
Appendix  C.  This  code  matches  the  chamber  pressure,  chamber  temperature,  and  mass 
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flow  between  NPSS  and  TDK  outputs.  The  non-dimensional  eontour  in  this  code  is 
dimensioned  with  the  throat  radius  of  4.46  in.  With  an  outer  radius  of  five  inches,  the 
throat  area  is  15.9  in^,  which  is  the  area  input  in  NPSS. 

4,3,  Investigative  Objectives 

This  investigation  tried  to  meet  three  objectives: 

1 .  Determine  feasibility  of  meeting  the  IHPRPT  Phase  III  orbit  transfer  vehicle 
goals  with  the  DEAN  concept 

2.  Implement  and  improve  upon  a  design  process  focused  on  the  energy 
conversion  section  of  a  rocket  engine  (combustion  chamber,  nozzle) 

3.  Perform  detailed  design  analysis  of  the  energy  transfer  components  (cooling 
jackets)  making  the  DEAN  possible 

The  result  of  the  NPSS  model  showed  that  this  model  is  able  to  meet  the 


IHPRPT  Phase  III  goals  and  detailed  to  analysis  of  the  energy  conversion  and  energy 
transfer  components.  The  result  is  the  conceptual  solid  model  of  the  DEAN  shown  in 
Eigure  21. 


Figure  21,  The  DEAN 
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The  left  portion  of  Figure  21  shows  the  solid  model  while  the  right  portion  depiets 
a  line  drawing  showing  the  inside  of  the  ehamber.  Figure  21  illustrates  the  full-length 
aerospike  nozzle  developed  by  TDK  and  the  ehamber  designed  by  NPSS.  The  outer  wall 
of  the  ehamber  eonsists  of  a  silieon  earbide  eooling  jaeket  with  an  estimated  mass  of  5.96 
Ibm  (2.7  kg).  Around  the  cooling  jaeket  is  a  structural  jacket.  The  structural  jacket  is 
sized  as  a  pressure  vessel  using  Equation  18: 


In  Equation  18,  tw  is  the  wall  thickness,  rc  is  the  radius  of  the  chamber,  Etu  is  the 
ultimate  tensile  strength  of  the  material,  and  Pb  is  the  burst  pressure.  The  bust  pressure  is 
the  chamber  pressure  multiplied  by  a  factor  of  safety,  in  the  case  of  the  DEAN,  the  factor 
of  safety  is  1.5.  Different  materials  were  evaluated  to  minimize  thickness  and  weight. 
Aluminum  7075-T6,  with  an  Em  of  80  ksi  (555.1  MPa),  was  judged  the  best  with  a 
thickness  of  0.2  in  (0.5  cm)  and  a  mass  of  13.2  Ibm  (6  kg).  The  cooled  part  of  the  nozzle 
and  inner  wall  is  composed  of  copper  resulting  in  a  mass  of  about  69  Ibm  (31.3  kg).  The 
un-cooled  portion  of  the  nozzle  is  niobium  and  has  a  mass  of  2.3  Ibm  (1 .05  kg).  The  total 
mass  of  the  chamber/nozzle  is  estimated  at  90.5  Ibm  (41.05  kg).  Combined  with  the 
oxygen  turbo  machinery  at  137.3  Ibm  (62.3  kg)  and  hydrogen  turbo  machinery  of  251.1 
Ibm  (1 13.9  kg)  results  is  an  estimated  total  mass  of  478.9  Ibm  (217.2  kg)  for  the  DEAN 
rocket  engine  concept.  This  results  in  a  thrust  to  weight  ratio  of  1 19: 1 . 


4,4,  Summary 

NPSS  was  used  to  model  the  DEAN.  This  model  resulted  in  several  significant 
results.  The  DEAN  model  produces  57,23 1  Ibf  (254  kN)  of  thrust  with  an  fip  of  472.7 
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sec.  The  hydrogen  cools  the  inner  wall  of  the  chamber  and  part  of  the  nozzle  with  a  mass 
flow  of  7.55  Ibm/sec  (3.42  kg/sec).  The  cooling  jacket  increases  the  flow  temperature  by 
464  R  (257.8  K).  The  walls  are  made  of  copper  and  have  a  maximum  temperature  of 
1224  R  (680  K)  at  the  throat.  The  oxygen  cools  the  outer  wall  of  the  chamber  with  a 
mass  flow  of  106  Ibm/sec  (48.1  kg/sec).  The  cooling  jacket  increases  the  flow 
temperature  by  437.4  R  (243  K).  The  cooling  channels  walls  consist  of  silicon  carbide 
and  have  a  maximum  temperature  of  1625  R  (902.8  K)  at  the  throat.  Due  to  the  accuracy 
of  NPSS,  it  is  expected  that  the  numbers  presented  would  closely  match  the  results  and 
actual  working  engine. 
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5.  Conclusions  and  Recommendations 


5.1.  Chapter  Overview 

The  DEAN  is  a  dual-expander  aerospike  nozzle  design  providing  superior 
performanee  to  the  eurrent  upper  stage  fleet.  The  eonelusion  and  reeommendations  based 
on  the  work  done  in  this  doeument  are  presented  here. 

5.2.  Conclusions  of  Research 

Based  on  the  results  presented  in  Chapter  4,  the  DEAN  would  meet  and  exeeed 
the  IHPRPT  phase  III  goals.  Additionally,  a  dual-expander  eyele  with  separate  turbines 
for  the  fuel  and  oxidizer  flow  is  possible  sinee  the  amount  of  heat  eaptured  from  the 
ehamber  is  suffleient  to  inerease  the  flows  to  the  required  temperature.  Additionally, 
though  the  learning  eurve  is  steep,  NPSS  has  proven  to  be  a  powerful  tool  in  the 
development  of  roeket  engines.  TDK  was  easier  to  learn  in  eomparison,  however,  its 
ability  to  model  aerospike  nozzles  was  left  wanting  due  to  the  eonfusing  kp  result. 

5.3.  Significance  of  Research 

The  work  eontained  within  this  doeument  eould  result  in  signifieant  gains  in 
performanee  for  the  Delta  and  Atlas  platforms.  Based  on  the  eurrent  Centaur  upper  stage, 
the  DEAN  model  would  result  in  a  28.3  s  inerease  in  kp.  This  eould  result  in  as  mueh  as 
a  3,806  Ibm  (1726  kg)  inerease  in  satellite  mass  to  orbit.  In  addition,  the  DEAN  has  the 
potential  to  save  around  $19  million  per  launeh  based  on  the  savings  per  kp  presented  in 
Chapter  1. 
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5,4,  Recommendations  for  Future  Research 


The  work  presented  in  this  doeument  is  for  the  initial  design  eoneepts.  Therefore, 
it  is  reeommended  further  researeh  be  done  based  on  this  model  of  the  DEAN  to 
eomplete  the  design  proeess.  For  one,  the  diserepaney  between  the  kp  from  the  NPSS 
and  TDK  model  should  be  resolved.  Work  was  done  on  an  elements  for  NPSS  that 
allows  NPSS  to  run  TDK.  At  this  point  the  element  has  not  suoeessfully  run.  Getting 
this  element  to  run  should  result  in  a  improve  idea  of  the  kp  of  the  DEAN  design  and 
whether  truneation  of  the  nozzle  is  neeessary. 

Even  though  NPSS  is  a  powerful  tool,  it  has  some  limitation.  For  example,  it 
eannot  show  the  effeets  the  eurvature  of  the  ehamber  will  have  on  the  eooling  jaeket  and 
the  eombustion.  Therefore,  a  detailed  analysis  of  the  flow  within  the  eooling  jaeket  to 
eonfer  the  properties  of  the  flow  obtained  from  NPSS.  This  would  inelude  a  CFD 
analysis  of  the  fluid  flow  to  depiet  the  flow  filed  in  the  ehamber  as  well  as  the  eooling 
jaekets.  The  CFD  model  will  result  in  a  high  fidelity  model  for  the  properties  of  the 
ehamber  and  their  effeet  on  the  eooling  jaeket. 

Furthermore,  in  this  doeument  the  ineorporation  of  the  injeetor  design  only 
resulted  in  the  assumption  of  the  pressure  drop  aeross  the  injeetor  and  effieieney  of 
eombustion.  A  detailed  investigation  on  the  benefits  of  a  Pintle  or  Coaxial  injeetors  into 
the  DEAN  design  should  be  undertaken.  This  will  determine  the  type  of  injeetor  that 
should  be  used  on  the  design  and  a  detailed  design  of  the  injeetors. 

The  work  done  in  this  doeument  only  details  a  design  at  a  single  point.  Therefore, 
further  work  on  the  DEAN  should  involve  on  investigation  on  off-design  eonditions. 

This  is  required  to  outline  the  ability  for  the  model  to  meet  throttling  requirements.  To 
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ensure  the  DEAN  design  is  the  best  it  ean  be,  and  optimization  of  the  eurrent  design 
should  be  undertaken.  The  optimization  would  determine  the  aspeets  of  eaeh  eomponent 
that  results  in  the  overall  best  engine  design.  Aspects  to  consider  are,  but  not  limited  to, 
the  amount  of  truncation  of  the  nozzle,  optimal  0/F  and  s,  and  material  choices  for  the 
walls. 

Another  aspect  of  the  design  to  consider  is  the  possible  use  of  fuels  other  than 
hydrogen.  One  fuel  the  space  command  has  expressed  interest  in  is  the  use  of  methane 
since  it  is  easier  to  store  than  liquid  hydrogen.. 

5,5,  Summary 

In  conclusion,  with  the  incorporation  of  state  of  the  art  turbo  machinery,  and  the 
use  of  the  aerospike  nozzle,  the  proposed  upper  stage  design  of  the  DEAN  achieves  a 
performance  that  could  radically  enhance  the  current  space  launch  fleet.  The  DEAN 
model  resulted  in  57,231  Ibf  (254  kN)  of  thrust,  14%  over  the  goal,  and  an  kp  of  472.7 
sec,  1.9%  increase  over  the  goal.  In  addition,  the  walls  of  the  chamber  and  nozzle  were 
adequately  cooled  to  prevent  the  failure  and  provided  the  necessary  energy  to  the  fuel  and 
oxidizer  turbine. 
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Appendix  A  :  Lessons  Learned 


In  the  execution  of  the  experiment  detailed  in  this  report,  several  lessons  were 
learned.  Most  of  these  lessons  deal  with  the  use  of  NPSS  and  TDK. 

A,1  NPSS  Lessons  Learned 

NPSS  is  a  great  tool,  though  the  learning  curve  is  steep.  First,  it  is  best  not  to  use 
the  graphic-user-interface  (GUI)  when  coding.  The  GUI  is  useful  when  viewing  the 
schematic  of  the  code,  editing,  and  viewing  the  results  of  the  code.  However,  when  using 
it  to  write  the  code  it  often  adds  needless  lines  confusing  more  than  enhancing  the  code. 

A  major  revelation  was  the  elements  are  defined  by  ‘.inf  files  and  were  located  in 
a  folder  within  NPSS.  Before  this  revelation,  the  coefficients  NPSS  uses  were  confusing, 
seemed  unrelated  to  any  known  coefficients,  and  were  not  explained  anywhere  in  the 
documentation  provided  with  NPSS.  However,  in  the  ‘.inf  files,  the  equations  in  each 
element  are  defined  and  the  coefficient  values  were  revealed.  Though  the  user  input 
coefficients  were  convoluted,  the  results  of  the  inputs  were  predictable. 

One  drawback  is  NPSS  only  uses  English  units.  The  use  of  English  units  can 
unduly  confuse  the  problem.  It  would  have  been  useful  to  include  either  SI  elements,  or 
elements  able  to  use  either  unit. 

The  caveat  to  all  this  is  NPSS  allows  for  the  user  to  code  there  own  ‘.inf  files. 
Therefore,  any  complication  can  be  overcome  with  the  user  creating  their  own  elements. 
However,  this  effort  takes  away  from  the  actual  design  since  time  is  spent  fixing  NPSS. 
Additionally,  learning  NPSS  is  difficult  enough,  let  alone  creating  new  elements. 
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A.2  TDK  Lessons  Learned 


TDK  is  a  powerful  program  for  analyzing  nozzles.  One  useful  thing  would  have 
been  graphing  of  the  chamber  as  well  as  the  nozzle.  One  realization  in  using  the  code 
was  TDK  is  very  sensitive  to  the  wall  contour.  Unlike  some  programs  continuing  to 
work  even  if  the  physics  are  impossible,  TDK  aborts  the  code  if  the  wall  fails.  This  can 
be  frustrating  when  trying  to  design  a  nozzle  since  the  contour  is  unknown  often  leaving 
the  researcher  with  a  mystery  as  to  the  fix  to  the  problem. 

The  ability  for  TDK  to  model  aerospike  nozzles  is  relatively  unique.  However, 
the  aerospike  is  a  feature  of  the  scramjet  modeling  capability.  Some  complications  arise 
since  the  scramjet  is  an  aero-device  while  the  aerospike  is  a  rocket  device.  This  creates 
confusion  since  the  scramjet  burns  with  air  while  the  aerospike  does  not  use  air.  It  is 
recommended  future  releases  of  TDK  have  a  separate  aerospike  function. 


88 


Appendix  B  :  NPSS  Code 


This  is  the  NPSS  code  that  defines  the  DEAN  model.  This  code  requires  the  use 
of  modified  elements  developed  by  Tom  Lavelle  at  NASA  Glen.  When  opened  in  NPSS 
these  files  will  generate  data  presented  in  this  report. 


B,1  Final  Model  Code 

1/===^^=^=^^=====^^========^========^========^ 

H  Air  Force  Institute  of  Technology 
//  2950  Hobson  Way,  Bldg  641 
//  Wright  Patterson  AFB,  OH  45433 

II 

H  David  F.  Martin  II,  2Lt.,  USAF 
11============================================ 

//============================================ 

H  Goals 

11============================================ 

real  mdot  =  121.25; 
real  P_c=  1740.5; 
real  0_F  =  7; 
real  A_throat  =  15.9; 
real  epsilon  =  125; 

#include  "Pump.int"; 
include  "Turb02.int"; 

#include  "CoolingVolume.int"; 
real  mdotO  =106; 
real  mdotH  =  15.1; 
real  Wall_ternp  =  900; 

//============================================ 

H  Chamber 

11============================================ 

setDefaultComposition("C_02_H2"); 

RocketCombI  COMB  { 
comp  =  "C_02_H2"; 
radius_tc  =  6.0; 
volume  =  2075.44; 

FuJ.comp  =  "HYDROGEN"; 

FLoxid.comp  =  "OXYGEN"; 

ThermalOutputPort  Hx_zone01  { 
areaFlow  =  32*PI; 
areaHx  =  180.96; 
radCurv  =  parent.radius_tc; 

}; 

ThermalOutputPort  Hx_zone02  { 

areaFlow  =  . COMB. HxzoneOI. areaFlow; 
areaHx  =  . COMB. Hx  zoneOI  .areaHx; 
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radCurv  =  parent.radius_tc; 

}; 

ThermalOutputPort  Hx  zoneOS  { 
areaFlow  =  97.193; 
areaHx  =  .COMB.HxzoneOI  .areaHx; 
radCurv  =  parent.radius_tc; 

}; 

ThermalOutputPort  Hx_zone04  { 
areaFlow  =  89.339; 
areaHx  =  . COMB. HxzoneOI  .areaHx; 
radCurv  =  parent.radius_tc; 

}; 

ThermalOutputPort  Hx_zone05  { 
areaFlow  =  44.76; 
areaHx  =  165.88; 
radCurv  =  5.5; 

}; 

ThermalOutputPort  Hx_zone06  { 

areaFlow  =  A_throat+3; 
areaHx  =  3.14; 
radCurv  =  5; 

}; 

ThermalOutputPort  Hx_zoneH1  { 
areaFlow  =  32*PI; 
areaHx  =  60.32; 
radCurv  =  2; 

}; 

ThermalOutputPort  Hx_zoneH2  { 

areaFlow  =  . COMB. HxzoneHI. areaFlow; 
areaHx  =  . COMB. HxzoneHI. areaHx; 
radCurv  =  . COMB. HxzoneHI. radCurv; 

}; 

ThermalOutputPort  Hx_zoneH3  { 
areaFlow  =  97.19; 
areaHx  =  67.86; 
radCurv  =  2.25; 

}; 

ThermalOutputPort  Hx_zoneH4  { 
areaFlow  =  89.33; 
areaHx  =  82.94; 
radCurv  =  2.75; 

}; 

ThermalOutputPort  Hx_zoneH5  { 
areaFlow  =  44.76; 
areaHx  =  120.64; 
radCurv  =  4; 

}; 

ThermalOutputPort  Hx_zoneH6  { 

areaFlow  =  A_throat+3; 
areaHx  =  2.83; 
radCurv  =  4.5; 

}; 

} 
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RocketNozzle  NOZZ  { 

Ath  =  A_throat; 

AR  =  epsilon; 

Ps  =  0.01; 

ThermalOutputPort  Hx_zoneH7  { 

areaFlow  =  A_throat*2; 
areaHx  =  87.96; 
radCurv  =  4; 

}; 

ThermalOutputPort  Hx_zoneH8  { 

areaFlow  =  A_throat*4; 
areaHx  =  54.98; 
radCurv  =  2.5; 

}; 

} 

//==============================. 

H  Opump 

11==============================. 

setDefaultComposition("OXYGEN"); 
Element  Starter  TankO  {} 

Element  Valve04  OD1  { 

Across  =  2; 

K  =  0.0116; 

} 

Element  CoolingVolume  OCV1  { 

UnReactedFluidlnputPort  Fl_l; 
UnReactedFluidOutputPort  Fl_0; 
volume  =  10; 

Aphys  =  2.112; 

} 

Element  Pump  OP  { 
gearRatIo  =  1; 

W  =  mdotO; 
eff=  0.773; 

PRdes  =  103; 

} 

Element  CoolingVolume  OCV2  { 

UnReactedFluidlnputPort  Fl_l; 
UnReactedFluidOutputPort  Fl_0; 
volume  =  10; 

Aphys  =  1.71; 

} 

//==============================-. 

H  Cooling  Jacketl 

11==============================-. 

Element  Valve04  DuckOI  { 

Across  =  0.47; 

K  =  0.0888; 

} 
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Element  CoolingVolume  CV01  { 

U  n  Reacted  FI  u  id  In  putPort  Fl_l  {OFR  =  1;}; 
UnReactedFluidOutputPort  Fl_0  {OFR  =  1;}; 
ThermalOutputPort  Hx  zonel  {areaHx  =  1.57;}; 
volume  =  0.05; 
n_channels  =  785; 

Aphys  =  0.47; 

} 

Element  Valve04  Duck02  { 

Across  =  .CV01  .Aphys; 

K  =  0.0879; 

} 

Element  CoolingVolume  CV02  { 
UnReactedFluidlnputPort  Fl_l  {OFR  =  1;}; 
UnReactedFluidOutputPort  Fl_0  {OFR  =  1;}; 
ThermalOutputPort  Hx  zonel  {areaHx  =  90.43;}; 
volume  =  3.62; 

n_channels  =  .CV01.n_channels; 

Aphys  =  0.75; 

} 

Element  Valve04  Duck03  { 

Across  =  .CV02.Aphys; 

K  =  0.165; 

} 

Element  CoolingVolume  CV03  { 
UnReactedFluidlnputPort  Fl_l  {OFR  =  1;}; 
UnReactedFluidOutputPort  Fl_0  {OFR  =  1;}; 
ThermalOutputPort  Hx  zonel  {areaHx  =  105.5;}; 
volume  =  4.22; 

n_channels  =  .CV01.n_channels; 

Aphys  =  0.88; 

} 

Element  Valve04  Duck04  { 

Across  =  .CV03.Aphys; 

K  =  0.172; 

} 

Element  CoolingVolume  CV04  { 
UnReactedFluidlnputPort  Fl_l  {OFR  =  1;}; 
UnReactedFluidOutputPort  Fl_0  {OFR  =  1;}; 
ThermalOutputPort  Hx  zonel  {areaHx  =  105.5;}; 
volume  =  4.22; 

n_channels  =  .CV01.n_channels; 

Aphys  =  0.88; 

} 

Element  Valve04  Duck05  { 

Across  =  .CV04  .Aphys; 
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K  =  0.127; 

} 

Element  CoolingVolume  CV05  { 

U  n  Reacted  FI  u  id  In  putPort  Fl_l  {OFR  =  1;}; 
UnReactedFluidOutputPort  Fl_0  {OFR  =  1;}; 
ThermalOutputPort  Hx  zonel  {areaHx  =  105.5;}; 
volume  =  4.22; 

n_channels  =  .CV01.n_channels; 

Aphys  =  0.88; 

} 

Element  Valve04  Duck06  { 

Across  =  .CV05.Aphys; 

K  =  0.0968; 

} 

Element  CoolingVolume  CV06  { 
UnReactedFluidlnputPort  Fl_l  {OFR  =  1;}; 
UnReactedFluidOutputPort  Fl_0  {OFR  =  1;}; 
ThermalOutputPort  Hx  zonel  {areaHx  =  105.5;}; 
volume  =  5.28; 

n_channels  =  .CV01.n_channels; 

Aphys  =  1.10; 

} 

Element  Valve04  Duck07  { 

Across  =  .CV06.Aphys; 

K  =  0.125; 

} 

Element  Wall2  OwalH  { 

Thermal InputPort  CoolSidel; 

ThermalInputPort  HotSidel; 

} 

Element  Wall2  Owall2  { 

ThermalInputPort  CoolSidel; 

ThermalInputPort  HotSidel; 

} 

Element  Wall2  Owall3  { 

ThermalInputPort  CoolSidel; 

ThermalInputPort  HotSidel; 

} 

Element  Wall2  Owall4  { 

ThermalInputPort  CoolSidel; 

ThermalInputPort  HotSidel; 

} 

Element  Wall2  Owall5  { 

ThermalInputPort  CoolSidel; 

ThermalInputPort  HotSidel; 

} 

Element  Wall2  Owall6  { 

ThermalInputPort  CoolSidel; 

ThermalInputPort  HotSidel; 
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} 

//===================================== 

Element  CoolingVolume  OCV3  { 

U  n  Reacted  FI  u  id  In  putPort  Fl_l; 
UnReactedFluidOutputPort  Fl_01,  Fl_02; 
volume  =  10; 

Aphys  =  26.42; 

} 

1/^=====^^^^^=========^==^======^====== 

H  OT  bypass 

11^=====^^^^==========^=========^^=^=== 

Element  Valve04  TB1  { 

Across  =  2; 

K=  194.5; 

} 

Element  CoolingVolume  TBCV1  { 

UnReactedFluidlnputPort  Fl_l; 
UnReactedFluidOutputPort  Fl_0; 
volume  =  10; 

Aphys  =  2; 

} 

Element  Valve04  TB2  { 

Across  =  2; 

K=  194.5; 

} 

1/=======^============^========^=^===== 

Element  Turb02  OT  { 

Wflow  =  95; 

FI_O.Aphys  =  14.85; 
eff  =  .949; 

PR=  1.8; 

} 

Element  CoolingVolume  OCV4  { 
UnReactedFluidlnputPort  FI_I1  {OFR  =  1;}; 
UnReactedFluidlnputPort  FI_I2  {OFR  =  1;}; 
UnReactedFluidOutputPort  FI  O  {OFR  =  1;}; 
volume  =  10; 

Aphys  =  14.85; 

} 

Element  Valve04  OV  { 

Across  =  2; 

K  =  0.5255; 

} 

Element  Shaft  OSHAFT  { 

ShaftInputPort  Sh_l1,  Sh_l2; 

Nmech  =  32000; 

} 

1/================^===^=============^== 
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//  Hpump 

//===================================== 

setDefaultCompositionC'HYDROGEN"); 

Element  Starter  TankH  {} 

Element  Valve04  HD1  { 

Across  =  2; 

K  =  0.0336; 

} 

Element  CoolingVolume  HCV1  { 

U  n  Reacted  FI  u  id  In  putPort  Fl_l; 
UnReactedFluidOutputPort  Fl_0; 
volume  =  10; 

Aphys  =  2; 

} 

Element  Pump  HP1  { 
gearRatIo  =  1; 

W  =  mdotH; 
eff=  .8; 

PRdes  =  45; 

} 

Element  CoolingVolume  HCV2  { 

UnReactedFluidlnputPort  Fl_l; 
UnReactedFluidOutputPort  Fl_01,  Fl_02; 
volume  =  10; 

Aphys  =  2; 

} 

1/=========^^============^=======^=^=== 

H  Bypass 

11=====^===^^=========^==========^===== 

Element  Valve04  BD1  { 

Across  =  2; 

K  =  4; 

} 

Element  CoolingVolume  BCV1  { 

UnReactedFluidlnputPort  Fl_l; 
UnReactedFluidOutputPort  Fl_0; 
volume  =  10; 

Aphys  =  2; 

} 

Element  Valve04  BD2  { 

Across  =  2; 

K  =  4; 

} 

//===================================== 

Element  Pump  HP2  { 
gearRatIo  =  1; 

W  =  mdotH/2; 
eff=  .83; 
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PRdes  =  2; 


} 

Element  CoolingVolume  HCV3  { 

U  n  Reacted  FI  u  id  In  putPort  Fl_l; 
UnReactedFluidOutputPort  Fl_0; 
volume  =  10; 

Aphys  =  2; 


//=========================== 

//  Cooling  Jacket2 

1/^=======^^==========^^===== 

Element  Valve04  DuckHI  { 

Across  =  0.19; 

K  =  0.0633; 

} 

Element  CoolingVolume  CVH1  { 
UnReactedFluidlnputPort  FIJ  {}; 
UnReactedFluidOutputPort  Fl_0  {}; 
ThermalOutputPort  Hx_zone1  { 
areaHx  =  32.97; 

}; 

volume  =  .66; 
n_channels  =  314; 

Aphys  =  0.19; 

} 

Element  Valve04  DuckH2  { 

Across  =  .CVH1. Aphys; 

K  =  0.0474; 

} 

Element  CoolingVolume  CVH2  { 
UnReactedFluidlnputPort  FIJ  {}; 
UnReactedFluidOutputPort  Fl_0  {}; 
ThermalOutputPort  Hx_zone1  { 
areaHx  =  65.94; 

}; 

volume  =  1 .32; 
n_channels  =  314; 

Aphys  =  0.38; 

} 

Element  Valve04  DuckH3  { 

Across  =  .CVH2.Aphys; 

K  =  0.114; 

} 

Element  CoolingVolume  CVH3  { 
UnReactedFluidlnputPort  FIJ  {}; 
UnReactedFluidOutputPort  Fl_0  {}; 
ThermalOutputPort  Hx_zone1  { 
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areaHx  =  2.2; 

}; 

volume  =  0.04; 
n_channels  =  314; 

Aphys  =  0.44; 

} 

Element  Valve04  DuckH4  { 

Across  =  .CVHS.Aphys; 

K  =  0.149; 

} 

Element  CoolingVolume  CVH4  { 

U  n  Reacted  FI  u  id  In  putPort  FIJ  {}; 
UnReactedFluidOutputPort  Fl_0  {}; 
ThermalOutputPort  Hx_zone1  { 
areaHx  =  90.43; 

}; 

volume  =  1 .81 ; 
n_channels  =  314; 

Aphys  =  0.38; 

} 

Element  Valve04  DuckHS  { 

Across  =  .CVH4.Aphys; 

K  =  0.0773; 

} 

Element  CoolingVolume  CVH5  { 
UnReactedFluidlnputPort  FIJ  {}; 
UnReactedFluidOutputPort  Fl_0  {}; 
ThermalOutputPort  Hx_zone1  { 
areaHx  =  52.74; 

}; 

volume  =  1 .58; 
n_channels  =  314; 

Aphys  =  0.33; 

} 

Element  Valve04  DuckH6  { 

Across  =  .CVH5.Aphys; 

K  =  0.0511; 

} 

Element  CoolingVolume  CVH6  { 
UnReactedFluidlnputPort  FIJ  {}; 
UnReactedFluidOutputPort  Fl_0  {}; 
ThermalOutputPort  Hx_zone1  { 
areaHx  =  37.68; 

}; 

volume  =  1.51; 
n_channels  =  314; 

Aphys  =  0.31 ; 

} 
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Element  Valve04  DuckH7  { 

Across  =  .CVH6.Aphys; 

K  =  0.041; 

} 

Element  CoolingVolume  CVH7  { 

U  n  Reacted  FI  u  id  In  putPort  FIJ  {}; 
UnReactedFluidOutputPort  Fl_0  {}; 
ThermalOutputPort  Hx_zone1  { 
areaHx  =  30.14; 

}; 

volume  =  1.51; 
n_channels  =  314; 

Aphys  =  0.31; 

} 

Element  Valve04  DuckHS  { 

Across  =  .CVH7.Aphys; 

K  =  0.0379; 

} 

Element  CoolingVolume  CVH8  { 
UnReactedFluidlnputPort  FIJ  {}; 
UnReactedFluidOutputPort  Fl_0  {}; 
ThermalOutputPort  Hx_zone1  { 

areaHx  =  .CVH7.Hx_zone1. areaHx; 

}; 

volume  =  .CVH7. volume; 
n_channels  =  314; 

Aphys  =  .CVH7.Aphys; 

} 

Element  Valve04  DuckH9  { 

Across  =  .CVHS.Aphys; 

K  =  0.0351; 

} 

Element  Wall2  Hwalll  { 

Thermal InputPort  CoolSidel; 
ThermalInputPort  HotSidel; 

} 

Element  Wall2  Hwall2  { 

ThermalInputPort  CoolSidel; 
ThermalInputPort  HotSidel; 

} 

Element  Wall2  Hwall3  { 

ThermalInputPort  CoolSidel; 
ThermalInputPort  HotSidel; 

} 

Element  Wall2  Hwall4  { 

ThermalInputPort  CoolSidel; 
ThermalInputPort  HotSidel; 

} 
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Element  Wall2  Hwall5{ 

Thermal InputPort  CoolSidel; 
ThermalInputPort  HotSidel; 

} 

Element  Wall2  Hwall6  { 

ThermalInputPort  CoolSidel; 
ThermalInputPort  HotSidel; 

} 

Element  Wall2  Hwall7  { 

ThermalInputPort  CoolSidel; 
ThermalInputPort  HotSidel; 

} 

Element  Wall2  Hwall8  { 

ThermalInputPort  CoolSidel; 
ThermalInputPort  HotSidel; 

} 

//================================== 

Element  CoolingVolume  HCV4  { 

UnReactedFluidlnputPort  Fl_l; 
UnReactedFluidOutputPort  Fl_0; 
volume  =  10; 

Aphys  =  2; 

} 

Element  Turb02  HT  { 

Wflow  =  mdotH/2; 

FI_O.Aphys  =  2; 
eff=  .9; 

PR  =  1.85; 

} 

Element  CoolingVolume  HCV5  { 

UnReactedFluidlnputPort  FI_I1,  FI_I2; 
UnReactedFluidOutputPort  Fl_0; 
volume  =  10; 

Aphys  =  2; 

} 

Element  Valve04  HV  { 

Across  =  2; 

K  =  2; 

} 

Element  Shaft  HSHAFT  { 

ShaftInputPort  ShJI,  Sh_l2,  Sh_l3; 
Nmech  =  110000; 

} 

//================================== 

//links 

//================================== 

linkPorts("COMB.FI_tc",  "NOZZ.FI_l",  "Cl"); 

linkPorts("TankO.FI_0",  "ODI.FIJ",  "01"); 
linkPorts("OD1.FI_0",  "OCVI.FIJ",  "02"); 
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linkPorts("OCV1.FI_0",  "OP.FI_l",  "03"); 
linkPorts("OP.FI_0",  "0CV2.FI_I",  "04"); 

linkPorts("OCV2.FI_0",  "DuckOlFIJ",  "05"); 
linkPorts(  "Duck01.FI_0",  "CV01.FI_I",  "06" ); 
linkPorts(  "CV01.FI_0",  "Duck02.FI_l",  "07" ); 
linkPorts(  "Duck02.FI_0",  "CV02.FI_I",  "08" ); 
linkPorts(  "CV02.FI_0",  "Duck03.FI_l",  "09" ); 
linkPorts(  "Duck03.FI_0",  "CV03.FI_I",  "01 0" ); 
linkPorts(  "CV03.FI_0",  "Duck04.FI_l",  "Oil" ); 
linkPorts(  "Duck04.FI_0",  "CV04.FI_I",  "012" ); 
linkPorts(  "CV04.FI_0",  "Duck05.FI_l",  "013" ); 
linkPorts(  "Duck05.FI_0",  "CV05.FI_I",  "014" ); 
linkPorts(  "CV05.FI_0",  "Duck06.FI_l",  "015" ); 
linkPorts(  "Duck06.FI_0",  "CV06.FI_I",  "016" ); 
linkPorts(  "CV06.FI_0",  "Duck07.FI_l",  "017" ); 
linkPorts(  "Duck07.FI_0",  "0CV3.FI_I",  "018" ); 

linkPorts("OCV3.FI_02",  "TBI.FIJ",  "T1"); 
linkPorts("TB1.FI_0",  "TBCVI.FIJ",  "T2"); 
linkPorts("TBCV1.FI_0",  "TB2.FI_I",  "T3"); 
linkPorts("TB2.FI_0",  "0CV4.FI_I2",  "T4"); 

linkPorts("OCV3.FI_01",  "OT.FI_l",  "019"); 
linkPorts("OT.FI_0",  "0CV4.FI_I1",  "020"); 
linkPorts("OCV4.FI_0",  "OV.FIJ",  "021"); 
linkPorts("OV.FI_0",  "COMB.FI_oxid",  "022"); 

linkPorts("OP.Sh_0",  "OSHAFT.ShJI",  "SOI"); 
linkPorts("OT.Sh_0",  "0SHAFT.Sh_l2",  "S02"); 

linkPorts("CV01.Hx_zone1","Owall1.CoolSide1", 
linkPorts("COMB.Hx_zone06","Owall1  .HotSidel ' 
linkPorts("CV02.Hx_zone1","Owall2.CoolSide1", 
linkPorts("COMB.Hx_zone05","Owall2. HotSidel' 
linkPorts("CV03.Hx_zone1","Owall3.CoolSide1", 
linkPorts("COMB.Hx_zone04","Owall3. HotSidel' 
linkPorts(''CV04.Hx_zone1'',''Owall4.CoolSide1'', 
linkPorts(''COMB.Hx_zone03'',''Owall4. HotSidel' 
linkPorts(''CV05.Hx_zone1'',''Owall5.CoolSide1'', 
linkPorts(''COMB.Hx_zone02'',''Owall5.HotSide1' 
linkPorts(''CV06.Hx_zone1'',''Owall6.CoolSide1'', 
linkPorts(''COMB.Hx_zone01 ", ''Owall6. HotSidel ' 

linkPortsC'TankH.FLO",  "HDI .FIJ",  "FI"); 
linkPorts(''HD1.FI_0'',  "HCVIFIJ",  ''F2''); 
linkPorts(''HCV1.FI_0'',  "HPIFIJ",  ''F3''); 
linkPorts(''HP1.FI_0'',  ''HCV2.FI_I'',  ''F4''); 
linkPorts(''HCV2.FI_01'',  ''HP2.FI_I'',  ''F5''); 

linkPorts(''HCV2.FI_02'',  "BDIFIJ",  "Bl"); 
linkPorts(''BD1.FI_0'',  "BCVIFIJ",  ''B2''); 
linkPorts(''BCV1.FI_0'',  ''BD2.FI_I'',  ''B3''); 
linkPorts(''BD2.FI_0'',  ''HCV5.FI_I2'',  ''B4''); 


"Wl"); 

,  ''W2''); 
''W3''); 

,  ''W4''); 
''W5''); 

,  ''W6''); 
"W7''); 

,  ''W8''); 
"W9''); 

,  "WIO"); 

"WII"); 

,  ''W12''); 
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linkPorts("HP2.FI_0",  "HCVS.FIJ",  "F6"); 
linkPorts("HCV3.FI_0",  "DuckHlFIJ",  "F8"); 

linkPorts(  "DuckH1.FI_0",  "CVHI.FIJ",  "H2" ); 
linkPorts(  "CVH1.FI_0",  "DuckH2.FI_l",  "H3" ); 
linkPorts(  "DuckH2.FI_0",  "CVH2.FI_I",  "H4" ); 
linkPorts(  "CVH2.FI_0",  "DuckH3.FI_l",  "H5" ); 
linkPorts(  "DuckH3.FI_0",  "CVH3.FI_I",  "H6" ); 
linkPorts(  "CVH3.FI_0",  "DuckH4.FI_l",  "H7" ); 
linkPorts(  "DuckH4.FI_0",  "CVH4.FI_I",  "H8" ); 
linkPorts(  "CVH4.FI_0",  "DuckHS.FIJ",  "H9" ); 
linkPorts(  "DuckH5.FI_0",  "CVHS.FIJ",  "H10" ); 
linkPorts(  "CVH5.FI_0",  "DuckH6.FI_l",  "H11" ); 
linkPorts(  "DuckH6.FI_0",  "CVH6.FI_I",  "H12" ); 
linkPorts(  "CVH6.FI_0",  "DuckH7.FI_l",  "H13" ); 
linkPorts(  "DuckH7.FI_0",  "CVH7.FI_I",  "H14" ); 
linkPorts(  "CVH7.FI_0",  "DuckH8.FI_l",  "H15" ); 
linkPorts(  "DuckH8.FI_0",  "CVH8.FI_I",  "H16" ); 
linkPorts(  "CVH8.FI_0",  "DuckH9.FI_l",  "H17" ); 
linkPorts(  "DuckH9.FI_0",  "HCV4.FI_I",  "H18" ); 

linkPorts("HCV4.FI_0",  "HT.FI_I",  "F12"); 
linkPorts("HT.FI_0",  "HCV5.FI_I1",  "F13"); 
linkPorts("HCV5.FI_0",  "HV.FIJ",  "F15"); 
linkPorts("HV.FI_0",  "COMB.FuJ",  "F16"); 

linkPorts("HP1.Sh_0",  "HSHAFT.ShJI",  "SH1"); 
linkPorts("HT.Sh_0",  "HSHAFT.Sh_l2",  "SH2"); 
linkPorts("HP2.Sh_0",  "HSHAFT.Sh_l3",  "SH3"); 

linkPorts("CVH1.Hx_zone1","Hwall1.CoolSide1",  "W13"); 
linkPorts("NOZZ.Hx_zoneH8","Hwall1.HotSide1",  "W14"); 
linkPorts("CVH2.Hx_zone1","Hwall2.CoolSide1",  "W15"); 
linkPorts("NOZZ.Hx_zoneH7","Hwall2.HotSide1",  "W16"); 
linkPorts("CVH3.Hx_zone1","Hwall3.CoolSide1",  "W17"); 
linkPorts("COMB.Hx_zoneH6","Hwall3.HotSide1",  "W18"); 
linkPorts("CVH4.Hx_zone1","Hwall4.CoolSide1",  "W19"); 
linkPorts("COMB.Hx_zoneH5","Hwall4.HotSide1 ",  "W20"); 
linkPorts("CVH5.Hx_zone1","Hwall5.CoolSide1",  "W21"); 
linkPorts("COMB.Hx_zoneH4","Hwall5.HotSide1 ",  "W22"); 
linkPorts("CVH6.Hx_zone1","Hwall6.CoolSide1",  "W23"); 
linkPorts("COMB.Hx_zoneH3","Hwall6.HotSide1 ",  "W24"); 
linkPorts("CVH7.Hx_zone1","Hwall7.CoolSide1",  "W25"); 
linkPorts("COMB.Hx_zoneH2","Hwall7.HotSide1 ",  "W26"); 
linkPorts("CVH8.Hx_zone1","Hwall8.CoolSide1",  "W27"); 
linkPorts("COMB.Hx_zoneH1","Hwall8.HotSide1",  "W28"); 
//=========================================== 

//Guesses 

//=========================================== 

COMB.OFR  =  0_F  /  (1  +  0_F); 

COMB.TtJc  =  6550; 

COMB.PtJc  =  P_c; 

COMB.FIJc.W  =  mdot; 
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TankO.Pt  =  46; 
TankO.Tt=  150; 

OD1.W  =  mdotO; 
OV.W  =  mdotO; 
DuckOl.W  =  mdotO; 
Duck02.W  =  mdotO; 
DuckOS.W  =  mdotO; 
Duck04.W  =  mdotO; 
DuckOS.W  =  mdotO; 
DuckOS.W  =  mdotO; 
Duck07.W  =  mdotO; 
TB1.W  =  mdotO-95; 
TB2.W  =  mdotO-95; 

OCV1.Pt  =  45; 
OCV1.ht  =  66; 
OCV2.Pt  =  4635; 
OCV2.ht  =  66; 
CV01.Pt  =  4497; 
CV01.ht  =  66; 
CV02.Pt  =  4359; 
CV02.ht  =  66; 
CV03.Pt  =  4221; 
CV03.ht  =  66; 
CV04.Pt  =  4083; 
CV04.ht  =  66; 
CV05.Pt  =  3945; 
CVOS.ht  =  66; 
CV06.Pt  =  3807; 
CVOS.ht  =  66; 
OCV3.Pt  =  3663; 
OCV3.ht  =  66; 
TBCV1.Pt  =  2849; 
TBCV1.ht  =  66; 
OCV4.Pt  =  2035; 
OCV4.ht  =  66; 

Owalll.T  =  Wall_temp; 
Owall2.T  =  Wall_temp; 
Owall3.T  =  Wall_temp; 
Owall4.T  =  Wall_temp; 
OwallS.T  =  Wall_temp; 
OwallS.T  =  Wall_temp; 

TankH.Pt  =  46; 
TankH.Tt  =  40; 

HD1.W  =  mdotH; 
HV.W  =  mdotH; 

BD1.W  =  mdotH/2; 
BD2.W  =  mdotH/2; 
DuckH1.W  =  mdotH/2; 
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DuckH2.W  =  mdotH/2; 

DuckHS.W  =  mdotH/2; 

DuckH4.W  =  mdotH/2; 

DuckHS.W  =  mdotH/2; 

DuckHS.W  =  mdotH/2; 

DuckH7.W  =  mdotH/2; 

DuckHS.W  =  mdotH/2; 

DuckHS.W  =  mdotH/2; 

HCV1.Pt  =  45; 

HCV1.ht  =  66; 

HCV2.Pt  =  2035; 

HCV2.ht  =  66; 

BCV1.Pt  =  2008; 

BCV1.ht  =  66; 

HCV3.Pt  =  4080; 

HCV3.ht  =  66; 

CVH1.Pt  =  4034; 

CVH1.ht  =  66; 

CVH2.Pt  =  3988; 

CVH2.ht  =  66; 

CVH3.Pt  =  3942; 

CVH3.ht  =  66; 

CVH4.Pt  =  3896; 

CVH4.ht  =  66; 

CVH5.Pt  =  3850; 

CVHS.ht  =  66; 

CVH6.Pt  =  3804; 

CVHS.ht  =  66; 

CVH7.Pt  =  3758; 

CVH7.ht  =  66; 

CVH8.Pt  =  3712; 

CVHS.ht  =  66; 

HCV4.Pt  =  3663; 

HCV4.ht  =  66; 

HCV5.Pt=  1980; 

HCVS.ht  =  66; 

HwalH.T  =  Wall_temp; 

Hwall2.T  =  Wall_temp; 

Hwall3.T  =  Wall_temp; 

Hwall4.T  =  Wall_temp; 

HwallS.T  =  Wall_temp; 

HwallS.T  =  Wall_temp; 

Hwall7.T  =  Wall_temp; 

HwallS.T  =  Wall_temp; 

//Solver  will  adjust  these  variables 
Independent  H_valve  { 
varName  =  "HPI.eff"; 
autoSetup  =  TRUE; 

} 

Independent  0_valve  { 
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varName  =  "OP.eff'; 
autoSetup  =  TRUE; 

} 

//===================================================================== 

H  Solver 

//====================================================================. 

solver.solutionMode  =  "STEADY_STATE"; 
presolverSequence  =  {}; 

solverSequence  =  {"NOZZ",  "TankO",  "OD1",  "OP",  "DuckOI", 

"Duck02",  "DuckOS",  "Duck04",  "DuckOS",  "Duck06",  "Duck07",  "TB1",  "TB2",  "OT",  "OV", 
"TankH",  "HD1",  "HP1",  "BD1",  "BD2",  "HP2",  "DuckHI",  "DuckH2",  "DuckH3", 

"DuckH4",  "DuckHS",  "DuckH6",  "DuckH7",  "DuckHS",  "DuckH9",  "HT",  "HV", 

"COMB",  "OCV1",  "OCV2",  "CV01",  "CV02",  "CV03",  "CV04",  "CV05",  "CV06",  "OCV3", 
"TBCV1", 

"OCV4",  "OSHAFT",  "OwalU",  "Owall2",  "Owall3",  "Owall4",  "OwallS",  "Owall6", 

"HCV1",  "HCV2",  "BCV1",  "HCV3",  "CVH1",  "CVH2",  "CVH3",  "CVH4",  "CVH5", 

"CVH6",  "CVH7",  "CVH8",  "HCV4",  "HCV5",  "HSHAFT","Hwall1",  "Hwall2",  "Hwall3", 
"Hwall4",  "HwallS",  "HwallG",  "Hwall7",  "HwallS"} 

DataViewer  CaseRowViewer  Chamber  { 
titleBody  = 
titleVars  =  {} 

variableList  =  { "COMB.PtJc",  "COMB.TtJc",  "COMB.OF",  "NOZZ.Fg", 

"NOZZ.Isp",  "NOZZ.W" }; 

} 

DataViewer  CaseRowViewer  CoolO  { 
titleBody  = 
titleVars  =  {} 

variableList  =  { "CV01.FI_l.Tt",  "CV01.FI_O.Tt",  "CV02.FI_0.Tt", 

"CV03.FI_0.Tt",  "CV04.FI_0.Tt",  "CV05.FI_0.Tt",  "CV06.FI_0.Tt"}; 

} 

DataViewer  CaseRowViewer  CoolH  { 
titleBody  = 
titleVars  =  {} 

variableList  =  { "CVH1.FI_l.Tt",  "CVH1  .FI_O.Tt",  "CVH2.FI_0.Tt", 

"CVH3.FI_0.Tt",  "CVH4.FI_0.Tt",  "CVH5.FI_O.Tt",  "CVH6.FI_0.Tt",  "CVH7.FI_O.Tt", 
"CVH8.FI_0.Tt"}; 

} 

postsolverSequence  =  { "Chamber",  "CoolO",  "CoolH"}; 
autoSolverSetupO; 

solver.maxiterations  =  5000; 
solver.maxJacobians  =  5000; 

run(); 

Chamber.displayO; 

CoolO.displayO; 

CoolH. display(); 

cout «  "Oxygen  =============================="«endl; 

cout  «  "mass  Flow  "  «  CVOl.Wavg  «  "  Ibm/sec"  «  endl; 

cout  «  "max  temp  "  «  OwalH  .T  «  "  R, "  «  OwallG.T  «  "  R"«  endl; 


104 


cout  «  "velocity6  "  «  CV06.FI_0.V  «  "  ft/s"  «  endl; 
cout  «  "velocityS  "  «  CV05.FI_0.V  «  "  ft/s"  «  endl; 
cout  «  "velocity4  "  «  CV04.FI_0.V  «  "  ft/s"  «  endl; 

cout  «  "Sum  Q  "  «  OwalU  .Qout+Owall2.Qout+Owall3.Qout+Owall4.Qout+Owall5.Qout 
+Owall6.Qout «  "  Btu/sec"  «  endl; 

cout  «  "Total  SA "  «  CV01.Hx_zone1.areaHx+CV02.Hx_zone1.areaHx+ 
CV03.Hx_zone1  .areaHx+CV04.Hx_zone1  .areaHx+CV05.Hx_zone1  .areaHx+ 
CV06.Hx_zone1.areaHx  «  "  in''2"  «  endl  «  endl; 
cout  «  "OP  massFlow  "  «  OP.W  «  "  Ibm/s,  OP  PR  "  «  OP. PR  «  endl; 
cout  «  "ONmech  "  «  OSHAFT.Nmech  «  "  rpm"  «  endl; 

cout  «  "OT  massFlow  "  «  OT.Wflow  «  "  Ibm/s, "  «  "  OT  PR  "  «  OT.PR  «  endl; 
cout  «  "OT  power "  «  OT.pwr «  "  hp,  OP  power "  «  OP.pwr  «  "  hp"  «  endl; 
cout  «  "Hygrogen  ============================="«  endl; 

cout  «  "mass  Flow  "  «  CVHI.Wavg  «  "  Ibm/sec"  «  endl; 
cout  «  "max  temp  "  «  Hwall3.T  «  "  R"  «  endl; 
cout  «  "velocityS  "  «  CVH8.FI_O.V  «  "  ft/s"  «  endl; 
cout  «  "velocity?  "  «  CVH7.FI_O.V  «  "  ft/s"  «  endl; 
cout  «  "velocityS  "  «  CVH6.FI_O.V  «  "  ft/s"  «  endl; 

cout  «  "Sum  Q  "  «  Hwall1.Qout+Hwall2.Qout+Hwall3.Qout+Hwall4.Qout+Hwall5.Qout 
+Hwall6.Qout+Hwall7.Qout+Hwall8.Qout  «  "  Btu/sec"  «  endl; 
cout  «  "Total  SA "  «  CVH1.Hx_zone1.areaHx+CVH2.Hx_zone1.areaHx+ 
CVH3.Hx_zone1  .areaHx+CVH4.Hx_zone1  .areaHx+CVH5.Hx_zone1  .areaHx+ 
CVH6.Hx_zone1  .areaHx+CVH7.Hx_zone1  .areaHx+CVH8.Hx_zone1  .areaHx  « 

"  in''2"  «  endl  «  endl; 

cout  «  "HP1  massFlow  "  «  HP1  .W  «  "  Ibm/s,  HP1  PR  "  «  HP1  .PR  «  endl; 
cout  «  "HP2  massFlow  "  «  HP2.W  «  "  Ibm/s,  HP2  PR  "  «  HP2.PR  «  endl; 
cout  «  "HNmech  "  «  HSHAFT.Nmech  «  "  rpm"  «  endl; 
cout  «  "HT  massFlow  "  «  HT.Wflow  «  "  Ibm/s  "  «  "  HT  PR  "  «  HT.PR  «  endl; 
cout  «  "HT  pwr "  «  HT.pwr «  "  hp,  HP1  pwr "  «  HP1  .pwr  «  "  hp,  HP2  pwr "  « 
HP2.pwr «  "  hp"  «  endl; 


B,2  Included  Cooling  Volume  Element 

#ifndef  _COOLINGVOLUME_ 

#define  _COOLINGVOLUME_ 

yy************************************************************************* 

//  *  NASA  Glenn  Research  Center 
//*  21000  Brookpark  Rd 
//‘Cleveland,  OH  44135 
II* 

II  ************************************************************************* 

//#include  "/NPSS/dev/Rockets/Common/RocketIncludes.npss" 
#include  <lnterplncludes.ncp> 

class  CoolingVolume  extends  Element{ 

// - 

//  *******  DOCUMENTATION  ******* 

// - 

I* 

title  = 

M 
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COOLING  VOLUME  ELEMENT 


Any  number  of -->  I  | 

incoming  flows  |  | 

I  Cooling  Volume  |  -->  Any  number  of 

Any  number  of -->  I  |  exiting  flows 

Thermal InputPorts  |  | 


*/’ 

description  =  isA()  +  "  performs  mass  and  energy  storage  calculations.  Also  allows  for  heat 
transfer."; 

usageNotes  = 

M 

-  The  cooling  volume  element  accounts  for  mass  and  energy  storage  at  a 
volume  location  in  space. 

-  Any  number  of  Fluidinput,  FluidOutput,  and  Thermal  ports  can  be  requested 
by  the  user  at  run  time.  These  ports  are  then  connected  to  the  model  using 
the  usual  linkPorts  command. 

-  The  user  may  also  specify  and  external  heat  transfer  that  is  applied  to 
the  element.  Note  that  this  energy  will  come  from  outside  the  simulation 
and  will  break  continuity  of  energy  for  the  overall  system. 

-  The  purpose  of  this  element  is  to  calculate  mass  and  energy  derivatives 
that  are  used  by  the  solver  to  balance  the  conditions  and  this  location 

in  space. 

-  The  element  has  two  default  solver  independents  and  two  default  solver 
states. 

-  In  steady-state  mode  the  solver  will  vary  enthalpy  and  pressure 

until  the  mass  and  energy  derivatives  are  zero  (mass  in  =  mass  out,  energy 
in  =  energy  out). 

-  In  transient  mode,  the  pressure  and  enthalpy  will  be  varied  until  the 
calculated  internal  energy  and  density  equal  the  integrated  energy  and 
density. 

-  To  nullify  the  transient  effects  of  this  element  in  an  otherwise  transient 
solution,  set  the  solution  mode  of  the  element  to  steady  state.  Ex: 
elementName.solutionMode  =  'STEADY-STATE'.  Note  that  this  element  does 
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have  solutionMode  as  a  variable.  However,  giving  the  above  command  will 
set  all  of  the  pertinent  solver  objects  to  'STEADY-STATE'. 


// - 

//  *******  SETUP  VARIABLES  ******** 


real  Aphys{ 

lOstatus  =  INPUT;  units  =  INCH2; 

description  =  "Representative  cross  sectional  area  (used  to  calculate  statics)"; 
trigger  =  1 ; 

} 

string  comp{ 

description  =  "Composition  of  the  volume"; 
trigger  =  1 ; 

} 

real  drhotqdt{ 

value  =  0.0;  lOstatus  =  OUTPUT;  units  =  LBM_PER_FT3_SEC; 
description  =  "Time  derivative  of  total  density"; 

} 

real  dutqdt{ 

value  =  0.0;  lOstatus  =  OUTPUT;  units  =  BTU_PER_LBM_SEC; 
description  =  "Time  derivative  of  total  specific  internal  energy"; 

} 

intHset  { 

value  =  FALSE;  lOstatus  =  INPUT;  units  =  NONE; 

description  =  "Used  internally  to  determine  if  enthalpy  has  been  set  (User  Should  Not  Touch)"; 
hide(1); 

} 

real  ht{ 

lOstatus  =  OUTPUT;  units  =  BTU_PER_LBM; 
description  =  "Total  specific  enthalpy"; 
trigger  =  1 ; 

} 

real  Htln{ 

value  =  0.0;  lOstatus  =  OUTPUT;  units  =  BTU_PER_SEC; 
description  =  "Total  energy  flowing  in  to  the  volume  through  the  ports"; 

} 

real  HtOut{ 

value  =  0.0;  lOstatus  =  OUTPUT;  units  =  BTU_PER_SEC; 
description  =  "Total  energy  flowing  out  of  the  volume  through  the  ports"; 

} 

string  inportList[]; 
inportList  { 

description  ="list  of  all  fluid  input  ports"; 
ptrType  =  "UnReactedFluidlnputPort"; 

} 

int  n_channels  { 

value  =  0.0;  lOstatus  =  INPUT;  units  =  NONE; 
description  =  "Number  of  channels  in  the  cooling  volume"; 

} 

string  outportListQ; 
outportList  { 
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description  ="list  of  all  fluid  output  ports"; 
ptrType  =  "UnReactedFluidOutputPort"; 

} 

int  Pset  { 

value  =  FALSE;  lOstatus  =  INPUT;  units  =  NONE; 

description  =  "Used  internally  to  determine  if  pressure  has  been  set  (User  Should  Not  Touch)"; 
hide(1); 

} 

real  Pt{ 

lOstatus  =  OUTPUT;  units  =  PSIA; 
description  =  "Total  pressure"; 
trigger  =  1 ; 

} 

real  Qnet{ 

value  =  0.0;  lOstatus  =  OUTPUT;  units  =  BTU_PER_SEC; 
description  =  "Net  heat  tranfers  out  of  the  volume  through  thermal  ports"; 

} 

real  Qext{ 

value  =  0.0;  lOstatus  =  INPUT;  units  =  BTU  PER  SEC; 

description  =  "External  heat  transfer  supplied  by  user  (note  that  this  value  breaks  conservation 
of  energy"; 

} 

real  s_Qnet  { 

value  =  1.0;  lOstatus  =  INPUT;  units  =  NONE; 
description  ="Heat  transfer  scale  factor  (on  Q_net  only)"; 

} 

real  rhot{ 

value  =  0.0;  lOstatus  =  INPUT;  units  =  LBM  PER  FTS; 
description  =  "Density  based  on  total  conditions"; 

} 

real  ut{ 

value  =  0.0;  lOstatus  =  INPUT;  units  =  BTU  PER  LBM; 
description  =  "Total  specific  internal  energy"; 

} 

real  volume{ 

value  =  0.0;  lOstatus  =  INPUT;  units  =  INCH3; 
description  =  "Volume"; 

} 

real  Wavg{ 

value  =  0.0;  lOstatus  =  OUTPUT;  units  =  LBM  PER  SEC; 
description  =  "Average  mass  flowing  through  each  port"; 

} 

real  Win{ 

value  =  0.0;  lOstatus  =  OUTPUT;  units  =  LBM  PER  SEC; 
description  =  "Total  mass  flowing  into  the  volume"; 

} 

real  Wout{ 

value  =  0.0;  lOstatus  =  OUTPUT;  units  =  LBM  PER  SEC; 
description  =  "Total  mass  flowing  out  of  the  volume"; 

} 


// - 

//  ******  SETUP  PORTS,  FLOW  STATIONS,  SOCKETS,  TABLES  ******** 
// - 
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//  FLUID  PORTS 


UnReactedFlowStation  FI_vol{ 

description  =  "Internal  station  used  to  calculate  volume  properties"; 

} 

//  FUEL  PORTS 

//  BLEED  PORTS 

//  THERMAL  PORTS 
string  HxPorts[]; 

HxPorts  { 

description  =  "Array  of  Thermal  ports"; 
ptrType  =  "ThermalOutputPort"; 

} 

//  FLOW  STATIONS 
//  SOCKETS 

// - 

//  *******  INTERNAL  SOLVER  SETUP  ******* 


// - 

//  ******  ADD  SOLVER  INDEPENDENTS  &  DEPENDENT  ****** 

// - 

Independent  ind_Pt  { 
varName  =  "Pt"; 
autoSetup  =  TRUE; 

description  =  "Varies  the  volume  total  pressure"; 

} 

Independent  ind_ht  { 
varName  =  "ht"; 
autoSetup  =  TRUE; 

description  =  "Varies  the  volume  total  specific  enthalpy"; 
indepRef  =  "100."; 

} 

Integrator  integ_rho  { 
stateName  =  "rhot"; 
derivativeName  =  "drhotqdt"; 
eq_rhs  =  "Win"; 
eqjhs  =  "Wout"; 
autoSetup  =  TRUE; 

description  =  "Balances  the  volume  mass  storage"; 

} 


Integrator  integ_U  { 
stateName  =  "ut"; 
derivativeName  =  "dutqdt"; 
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eq_rhs  =  "Htin  +  Qnet  +  Qext"; 
eqjhs  =  "HtOut"; 
autoSetup  =  TRUE; 

//eq_Ref  =  "100."; 

description  =  "Balances  the  volume  energy  storage"; 

} 


// - 

//  *******  VARIABLE  CHANGED  METHODOLOGY  ******* 

// - 

// - 

//  *******  OPTION  VARIABLE  SETUP  ******* 

// - 

// - 

//  *******  PERFORM  ENGINEERING  CALCULATIONS  ******* 
// - 


// - 

//  calculate  function 

// - 

void  calculateO  { 
int  i; 

// - 

//  set  the  working  area 

// - 

FI_vol.Aphys  =  Aphys; 

// - 

//  Heat  Transfer  Calculations 
// - 

// - 

//  Determine  the  average  flow  rate 

// - 

Wavg  =  0.0; 
int  iport; 

for  ( iport  =  0;  iport  <  inportList.entries();  iport  =  iport  +  1  )  { 
Wavg  =  Wavg  +  inportList[iport]->W; 

} 

for  ( iport  =  0;  iport  <  outportList.entries();  iport  =  iport  +  1  )  { 
Wavg  =  Wavg  +  outportList[iport]->W; 

} 

Wavg  =  Wavg  /  2.; 

real  Wc  =  Wavg; 
real  Ac  =  Aphys; 

// - 

//  Compute  the  hydralic  diameter 
// - 


no 


real  D_hyd  =  2.  *(  Ac  /  PI  /  n_channels  )**0.5; 

Qnet  =  0; 

// - 

//  Loop  through  the  thermal  ports  calculating  the  heat  flux 

// - 

for  ( iport  =  0;  iport  <  HxPorts.entries();  iport  =  iport  +  1  )  { 

// - 

//  determine  the  heat  transfer  coefficient  for  the  heat 
//  flow  through  Hx_zone1 

// - 

real  he  =  s_Qnet  *  0.023  *(  abs(  Wc  )/  Ac  )**0.8 
*  Fl_vol.kt**0.6  *(  FLvol.Cpt/  FLvol.mut  )**0.4 
/  D_hyd**0.2; 


// - 

//  Compute  the  heat  transfer  rates 

//  Note  that  the  area  here  is  total  for  all  coolant  channels 

// - 

real  Q  =  he  *  HxPorts[iport]->areaHx  *  (  HxPorts[iport]->MassTemp  -  FI_vol.Tt ); 

// - 

//  Set  the  coefficients  and  rates  in  the  corresponding  thermal 
//  ports  connected  to  the  hot-wall  and  the  outer  wall 

// - 

HxPorts[iport]->HeatTransferCoef  =  he; 

HxPorts[iport]->HeatTransferRate  =  Q; 

// - 

//  add  the  heat  transfer  rates  to  the  volume  overall  Qhxdt 

// - 

Qnet  =  Qnet  +  Q; 


} 

// - 

//  Create  some  working  variables  and  initialize  values  for 
//  summations 

// - 

real  Wt; 

Win  =  0.0; 

Wout  =  0.0; 

Htin  =  0.0; 

HtQut  =  0.0; 

// - 

//  Grab  the  total  density  and  internal  energy  from  the  working 
//  station 

//  These  values  are  calculated  by  the  station  when  Pt  and  ht  are 
//  set 

// - 

rhot  =  FI_vol.rhot; 
ut  =  FI_vol.ut; 
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// - 

//  Sum  up  energy  and  flow  over  the  incoming  ports 

// - 

//  int  i; 

for  ( i=0;  i  <  inportList.entries()  ;  i  =  i  +  1  )  { 

Wt  =  inportList[i]->W; 
if  (  Wt  >=  0.0  ){ 

Win  =  Win  +  Wt; 

Htin  =  Htin  +  wt  *  inportList[i]->ht; 

} 

else  { 

Wout  =  Wout  -  Wt; 

HtOut  =  HtOut  -  Wt  *  ht; 

} 

} 

// - 

//  Sum  up  energy  and  flow  over  the  exiting  ports 

// - 

for  ( i=0;  i  <  outportList.entries();  i  =  i  +  1  )  { 

Wt  =  outportList[i]->W; 
if  (  Wt  >=  0.0  ){ 

Wout  =  Wout  +  Wt; 

HtOut  =  HtOut  +  Wt  *  ht; 

} 

else  { 

Win  =  Win  -  Wt; 

Htin  =  Htin  -  Wt  *  outportList[i]->ht; 


} 

} 

// - 

//  calculate  the  state-derivatives 

// - 

drhotqdt  =  (  Win  -  Wout )/  volume*12*12*12; 

dutqdt  =  (  Htin  +  Qnet  +  Qext  -  HtOut  -  ut  *  (  Win  -  Wout ))  /  ( rhot  *  volume/12/12/12  ); 

} 


// - 

//  VERIFY  function  (same  as  for  basic  Volume  class) 

// - 

int  verifyO 

{ 

// - 

//  set  the  composition  in  all  the  stations 

// - 

int  i; 

for  ( i  =  0;  i  <  inportList.entries();  i  =  i  +  1  )  { 
inportList[i]->comp  =  FI_vol.comp; 

} 

for  ( i  =  0;  i  <  outportList.entries();  i  =  i  +  1  )  { 
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outportList[i]->comp  =  FI_vol.comp; 


// - 

//  check  to  be  sure  the  volume  has  a  reasonable  value 


if  (volume  <=  0.0000  )  { 

cerr  «  "verify  failed  because  volume  has  not  been  set"  «  endl; 
return  FALSE; 


return  TRUE; 

} 


// - 

//  variable  changed  function 

// - 

void  variableChanged(  string  name,  any  oldValue  ){ 
int  i; 

// - 

//  Put  the  enthalpy  and  pressure,  which  are  probably  set  in  the 
//  solver  to  all  the  stations  and  ports 

// - 

if  (  name  ==  "comp" )  { 


FI_vol.comp  =  comp; 

//int  i; 

for  ( i  =  0;  i  <  inportList.entries();  i  =  i  +  1  )  { 
inportList[i]->comp  =  comp; 

} 

for  ( i  =  0;  i  <  outportList.entries();  i  =  i  +  1  )  { 
outportList[i]->comp  =  comp; 

} 

} 

} 


// - 

//  *******  PREPASS  FUNCTION  ******* 

// - 

void  prePassO  { 

//  Set  the  conditions  in  the  volume  and  all  the  ports.  The 
//  static  conditions  in  the  ports  will  depend  on  their  flow 
//  and  area, 
int  i; 

FI_vol.setTotal_hP(ht,  Pt); 

for  (i=0;  i<inportList.entries();  i  =  i  +  1)  { 
inportList[i]->setTotal_hP(ht,  Pt); } 
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for  (i=0;  i<outportList.entries();  i  =  i  +  1)  { 
outportList[i]->setTotal_hP(ht,  Pt); } 

} 


// - 

//  ************  POSTCREATE  FUNCTION  ********** 


void  postcreate(  string  name  )  { 

// - 

//  allow  for  the  creation  of  ports  at  run  time 

// - 

if  (  name->haslnterface(  "ThermalOutputPort" ))  { 
HxPorts.append(name); 

} 

if  (  name->haslnterface(  "UnReactedFluidlnputPort" ))  { 
inportList.append(name); 

} 

if  (  name->haslnterface(  "UnReactedFluidOutputPort" ))  { 
outportList.append(  name ); 

} 

} 


} 

//dumpMaps("maps.out"); 

#endif 


B,3  Included  Pump  Element 

#ifndef_PUMP_ 

#define  PUMP 


yy************************************************************************* 

//  *  NASA  Glenn  Research  Center 
//*  21000  Brookpark  Rd 
//‘Cleveland,  OH  44135 
11* 

II  ************************************************************************* 

//#include  "/NPSS/dev/Rockets/Common/RocketIncludes.npss" 
#include  <lnterplncludes.ncp> 

class  Pump  extends  Element{ 

// - 

//  *******  DOCUMENTATION  ******* 

// - 

/* 

title  = 

"  PUMP  ELEMENT 
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FIJ  -^1  Pump 


I— >  Fl_0 
I— >  Sh_0 


V 

socket:  S  Map 
socketType:  PUMP  MAP 
returns:  head,  trq 


description  =  isA()  +  "  calculates  the  performance  of  a  rocket  pump  element."; 
usageNotes  = 


-  This  element  works  by  determining  the  exit  conditions  based  on  either 
user  supplied  or  subelement  calculated  values  of  torque  and  head. 

-  The  element  takes  weight  flow  as  input  and  determines  an  exit  pressure. 

-  There  is  a  default  solver  independent  and  dependent  available  that  will  vary  the 
weight  flow  until  the  calculated  exit  pressure  matches  the  exit  pressure 

seen  at  the  port. 

-  The  user  must  supply  an  initial  guess  for  the  weight  flow. 

-  The  inertia  is  kept  in  the  mechanical  port.  The  user  should  set  the  value 
directly  there  (Sh_0. inertia). 


// - 

//  *******  SETUP  VARIABLES  ******** 

// - 

real  dht  { 

value  =  0.0;  lOstatus  =  OUTPUT;  units  =  BTU  PER  LBM; 
description  =  "Actual  specific  enthalpy  change"; 

} 

real  dhtideal  { 

value  =  0.0;  lOstatus  =  OUTPUT;  units  =  BTU  PER  LBM; 
description  =  "Ideal  specific  enthalpy  change"; 

} 

real  eff  { 

value  =  0.0;  lOstatus  =  OUTPUT;  units  =  NONE; 
description  =  "Compressor  efficiency"; 

} 

real  gearRatio  { 

value  =  1.0;  lOstatus  =  OUTPUT;  units  =  NONE; 
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description  =  "Gear  ratio  on  mechanical  connection"; 

} 

real  head  { 

value  =  0.0;  lOstatus  =  OUTPUT;  units  =  INCH; 
description  =  "Pump  head  diameter"; 

} 

real  Nmech  { 

value  =  0.0;  lOstatus  =  OUTPUT;  units  =  RPM; 
description  =  "Rotational  speed  (after  gear  ratio  applied)"; 

} 

real  PR  { 

value  =  0.0;  lOstatus  =  OUTPUT;  units  =  NONE; 

description  =  "Pressure  ratio  calculated  from  pump  characteristics"; 

} 

real  PRg{ 

value  =  0.0;  lOstatus  =  OUTPUT;  units  =  NONE; 

description  =  "Guess  for  pressure  ratio  calculated  from  fluid  ports"; 

} 

real  pwr  { 

value  =  0.0;  lOstatus  =  OUTPUT;  units  =  HORSEPOWER; 
description  =  "Compressor  power"; 

} 

real  trq  { 

value  =  0.0;  lOstatus  =  OUTPUT;  units  =  FT  LBF; 
description  =  "Compressor  torque"; 

} 

real  W  { 

value  =  0.0;  lOstatus  =  INPUT;  units  =  LBM  PER  SEC; 
description  =  "Compressor  weight  flow"; 

} 

real  PRdes; 

// - 

//  ******  SETUP  PORTS,  FLOW  STATIONS,  SOCKETS,  TABLES  ******** 
// - 

//  FLUID  PORTS 

UnReactedFluidlnputPort  Fl_l{ 
description  =  "Incoming  flow"; 

}; 


UnReactedFluidOutputPort  Fl_0{ 
description  =  "Exiting  flow"; 

}; 

UnReactedFlowStation  FI_Otemp{ 
description  =  "Exiting  flow"; 

}; 


//  FUEL  PORTS 
//  BLEED  PORTS 
//  THERMAL  PORTS 
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//  MECHANICAL  PORTS 


ShaftOutputPort  Sh_0{ 
description  =  "Mechanical  connection"; 

} 

//  FLOW  STATIONS 

UnReactedFlowStation  FI_Oi{ 
description  =  "Ideal  exit  conditions"; 

} 

//  SOCKETS 

Socket  S  map  { 
allowedValues  =  { "trq",  "head" } 
socketType=  "PUMP_MAP"; 
description  =  "Pump  performance  map"; 

} 


// - 

//  ******  ADD  SOLVER  INDEPENDENTS  &  DEPENDENT  ****** 

// - 

Independent  ind_W{ 
varName  =  "W"; 
autoSetup  =  TRUE; 
description  =  "Varies  the  weight  flow"; 

} 

Dependent  dep_PR{ 
eqjhs  =  "PR"; 
eq_rhs  =  "PRdes"; 
autoSetup  =  TRUE; 

description  =  "Compares  the  pressure  ratios  calculated  from  the  ports  and  the  map"; 

} 

// - 

//  *******  VARIABLE  CHANGED  METHODOLOGY  ******* 

// - 

// - 

//  *******  OPTION  VARIABLE  SETUP  ******* 

// - 

Option  switch  Des{ 

allowedValues  =  {  DESIGN,  OFFDESIGN  }  //default  is  DESIGN 
description  =  "Design  mode  switch  indicator  [DESIGN  /  OFFDESIGN]"; 
trigger  =  0; 

rewritableValues  =  FALSE;  //  Enables  converter  optimization. 

} 

// - 

//  *******  PERFORM  ENGINEERING  CALCULATIONS  ******* 

// - 

void  calculateO  { 
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// - 

//  Set  the  weight  flows  in  the  ports 

// - 

FI_I.W  =  W; 

FI_O.W  =  W; 

// - 

//  Determine  the  mechanical  speed  from  the  shaft  and  gear  ratio 

// - 

Nmech  =  gearRatio  *  Sh_O.Nmech; 

PR  =  FI_O.Pt/FI_I.Pt; 

real  Sout  =  FIJ.s; 

//real  PtOut  =  PR  *  Ptin; 

FI_Otemp.setTotalSP(  Sout,  FI_O.Pt ); 
real  htIdealOut  =  FI_Otemp.ht; 

// - 

//  Set  the  exit  conditions  (done  to  get  ht  right) 

// - 

FI_0.setTotal_hP(  FIJ.ht  +  (htIdealOut  -  Fl_l.ht)/eff,  FI_O.Pt ); 

dht  =  (  Fl_0.ht-Fl_l.ht ); 

// - 

//  Determine  the  power  for  output 

// - 

pwr  =  -dht  *  W  *  C_BTUtoFT_LBF  /  C_HPtoFT_LBF_PER_SEC; 

//cout  «  "pwr  is  "  «  pwr «  endl; 

trq  =  C_HP_PER_RPMtoFT_LBF  *  pwr  /  Nmech; 

// - 

//  Set  values  in  the  ports 

// - 

Sh_0.trq  =  trq; 

}//End  calculate 

}//End  element 

#endif 

BA  Included  Turbine  Element 

#ifndef_TURB02_ 

#define  TURB02 


yy************************************************************************* 

//  *  NASA  Glenn  Research  Center 
//*  21000  Brookpark  Rd 
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//‘Cleveland,  OH  44135 
II* 

II  ************************************************************************* 

#include  <lnterplncludes.ncp> 

class  Turb02  extends  Element{ 

// - 

//  *******  DOCUMENTATION  ******* 

// - 

I* 

title  = 

"  TURB02  ELEMENT 


I  |->  FLO 

FIJ— >|  Turb02  | 

I  I— >  Sh_0 


V 

socket:  S  map 
socketType:  TURB02_MAP 
return:  eff,  W 


*/ 

description  =  isA()  +  "  calculates  the  performance  of  a  turbine  that  is  being  driven  by  a  single 
constituent  fluid"; 

usageNotes  = 

M 

-  This  element  works  by  determining  the  exit  conditions  based  on  either 
user  supplied  or  subelement  calculated  values  of  efficiency  and  wieght 
flow. 

-  The  inertia  is  keep  in  the  mechanical  port.  The  user  should  set  the  value 
direcIty  there  (Sh_0. inertia). 

-  There  are  no  solver  independents  and  dependents  directly  associatted 
with  this  element. 


// - 

II  *******  SETUP  VARIABLES  ******** 

// - 

real  dht  { 

value  =  0.0;  lOstatus  =  OUTPUT;  units  =  BTU  PER  LBM; 
description  =  "Change  in  enthalpy"; 
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} 

real  dhtideal  { 

value  =  0.0;  lOstatus  =  OUTPUT;  units  =  BTU  PER  LBM; 
description  =  "Ideal  enthalpy  change"; 

} 

real  eff  { 

value  =  0.0;  lOstatus  =  OUTPUT;  units  =  NONE; 
description  =  "Efficiency"; 

} 

real  gearRatio  { 

value  =  1.0;  lOstatus  =  OUTPUT;  units  =  NONE; 
description  =  "Gear  ratio  on  mechanical  connection"; 

} 

real  Nmech  { 

value  =  0.0;  lOstatus  =  OUTPUT;  units  =  RPM; 
description  =  "Rotational  speed  (after  gear  ratio  applied)"; 

} 

real  Nrad  { 

value  =  0.0;  lOstatus  =  OUTPUT;  units  =  RAD  PER  SEC; 
description  =  "Rotational  speed  in  radians  per  sec"; 

} 

real  PR  { 

value  =  0.0;  lOstatus  =  OUTPUT;  units  =  NONE; 
description  =  "Pressure  ratio  (seen  from  ports)"; 

} 

real  pwr  { 

value  =  0.0;  lOstatus  =  OUTPUT;  units  =  HORSEPOWER; 
description  =  "Overall  power"; 

} 

real  trq  { 

value  =  0.0;  lOstatus  =  OUTPUT;  units  =  FT  LBF; 
description  =  "Torque  on  the  shaft"; 

} 

real  Wflow  { 

value  =  0.0;  lOstatus  =  INPUT;  units  =  LBM  PER  SEC; 
description  =  "Weight  flow"; 

} 

// - 

//  ******  SETUP  PORTS,  FLOW  STATIONS,  SOCKETS,  TABLES  ******** 


//  FLUID  PORTS 

UnReactedFluidlnputPort  Fl_l{ 
description  =  "Incoming  flow"; 

}; 


UnReactedFluidOutputPort  Fl_0{ 
description  =  "Exiting  flow"; 

}; 


//  FUEL  PORTS 
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//  BLEED  PORTS 


//  THERMAL  PORTS 

//  MECHANICAL  PORTS 

ShaftOutputPort  Sh_0{ 
description  =  "Mechanical  connection"; 

} 

//  FLOW  STATIONS 
UnReactedFlowStation  FI_Oi{ 
description  =  "Ideal  exit  conditions"; 

} 

//  SOCKETS 
Socket  S  map  { 

allowedValues  =  { "eff" ,  "Wflow" } 
description  =  "Turbine  map"; 
//socketType  =  TURB02_MAP; 

} 

Independent  ind_W{ 
varName  =  "Wflow"; 

//autoSetup  =  TRUE; 

description  =  "Varies  the  weight  flow"; 

} 


// - 

//  *******  VARIABLE  CHANGED  METHODOLOGY  ******* 

// - 

// - 

//  *******  OPTION  VARIABLE  SETUP  ******* 

// - 

Option  switch  Des{ 

allowedValues  =  {  DESIGN,  OFFDESIGN  }  //default  is  DESIGN 
description  =  "Design  mode  switch  indicator  [DESIGN  /  OFFDESIGN]"; 
trigger  =  0; 

rewritableValues  =  FALSE;  //  Enables  converter  optimization. 

} 

// - 

//  *******  PERFORM  ENGINEERING  CALCULATIONS  ******* 

// - 

void  calculateO  { 


// - 

//  Determine  the  mechanical  speed  from  the  shaft  and  gear  ratio 

// - 

Nmech  =  gearRatio  *  Sh_O.Nmech; 
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// - 

//  Switch  the  mechanical  speed  to  radians  per  sec 

// - 

Nrad  =  Nmech*2*PI/60; 

// - 

//  Determine  PR 

// - 

PR  =  FI_I.Pt/FI_O.Pt; 

// - 

//  determine  the  ideal  enthalpy  change 

// - 

FI_Oi.copyFlow(  "FIJ" ); 

FI_Oi.setTotalSP(  FIJ.s,  FI_O.Pt ); 
dhtldeal  =  FI_l.ht-FI_Oi.ht; 


// - 

//  Calculate  turbine  torque 

// - 

dht  =  dhtideal  *  eff; 

FI_0.setTotal_hP(  FIJ.ht  -  dht,  FI_O.Pt ); 

trq  =  C_BTUtoFT_LBF  *  Wflow  *  dht  /  Nrad; 

pwr  =  dht  *  Wflow  *  C_BTUtoFT_LBF  /  C_HPtoFT_LBF_PER_SEC; 

// - 

//  set  the  conditions  in  the  shaft  port 

// - 

Sh_0.trq  =  trq; 

// - 

//  set  the  flows 

// - 

FI_I.W  =  Wflow; 

FI_O.W  =  Wflow; 

}//End  calculate 

}//End  element 

#endif 
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Appendix  C  :  TDK  Code 


This  is  TDK  input  and  output  developed  for  use  with  the  DEAN  NPSS  model  file. 


C.l  TDK  Input  for  DEAN  Upper  Stage  Engine 

TITLE  DEAN_coiitourlOO% 

DATA 

$DATA 

SCRJET  =  1,  PLANAR  =  F, 

ODE  =  1,  ODK=l,  TDK=1,  MABL  =  2,  IMABL  =  0,  MABLK  =  T, 

RSI  =  2.084,  RWTD  =  0.05,  RZERO  =  -  I, 

IWALL=  2,  THETA  =  25,  RMAX  =  -O.OI,  ZMAX  =  3.2, 

LWALL  =  I,  THALW  =  2.5,  XLW(I)  =  0,  YLW(I)  =  -LI I,  XCOWL  =  O.I,  ICOWL  =  I, 
JET  =  T,  PINF  =  0.1,  MINF  =  16.4,  GINF  =  1.4, 

SEND 

REACTANTS 

H2.  100.  -2I54.L  20.27F  .0709 

0  2.  100.  -3I02.L  90.560  1.149 

NAMELISTS 

$ODE 

RKT  =  T,  PSIA  =  T,  XP  =  I,  DELH  =  0, 

OF  =  T,  OFSKED  =  7,  PSJ  =  1739,  TSJ  =  6586,  MSJ  =  3.0, 

SEND 

REACTIONS 

H  -r  H  =  H2  ,MI,  A  =  6.4EI7,  N  =  1.0,  B  =  0.0,  (AR)  BAULCH  72  (A)  30U 
H  -r  OH  =  H20  ,M2,  A  =  8.4E2I,  N  =  2.0,  B  =  0.0,  (AR)  BAULCH  72  (A)  lOU 
O  -r  O  =  02  ,M3,  A  =  I.9EI3,  N  =  0.0,  B  =-I.79,(AR)  BAULCH  76  (A)  lOU 

O  -r  H  =  OH  ,M7,  A  =3.62EI8,  N  =  1.0,  B  =  0.0,  (AR)  JENSEN  78  (B)  30U 
END  TBR  REAX 

02-rH  =  0-r0H  ,  A  =  2.2EI4, N  =  0.0,  B  =16.8,  BAULCH  72  (A)  1. 5U 
H2-rO  =  H-rOH  ,  A  =  I.8EI0, N  = -I.,  B  =  8.9,  BAULCH  72  (A)  1. 5U 
H2  -r  OH  =  H20  -r  H  ,  A  =  2.2EI3,  N  =  0.0,  B  =5.15,  BAULCH  72  (A)  2U 
OH  -r  OH  =  H20  -r  O  ,  A  =  6.3EI2,  N  =  0.0,  B  =1.09,  BAULCH  72  (A)  3U 
LAST  REAX 

THIRD  BODY  REAX  RATE  RATIOS 

MI  =  25*H,4*H2,I0*H2O,25*O,25*OH,I.5*O2, 

M2  =  I2.5*H,5*H2,I7*H20,I2.5*0,I2.5*0H,6*02, 

M3  =  I2.5*H,5*H2,5*H20,I2.5*0,I2.5*0H,II*02, 

M7  =  I2.5*H,5*H2,5*H20,I2.5*0,I2.5*0H,5*02, 

LAST  CARD 
SODK 
EP  =  125, 

SEND 

STRANS 

SEND 

SMOC 

SEND 

SMABL 

ADBATC=  0,  NTQW=2,  TQW=2*2000,  XTQW=-I.E6,I.E6, 

DXI=I0E-4,  NDXI=50,  NYI=I05,  DCIMAX  =  .005,  DXLIM(I)  =  .04,.75, 
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INNER  =  F, 

OFC=l,  DISTRB^O,  XC0=-2.16,  XCE=5.13, 

SEND 

$MABL 

ADBATC=  0,  NTQW=2,  TQW=2*2000,  XTQW=-1.E6,1.E6, 

DXI=10E-4,  NDXI=50,  NYI=105,  DCIMAX  =  .005,  DXLIM(l)  =  .04, .75, 

INNER  =  T, 

OFC=2,  DISTRB^O,  XC0=-2.I6,  XCE=5.I3, 

SEND 

OTITLE  DEAN_coiitourI00% 

ODATA 

SDATA 

SCRJET  =  1,  PLANAR  =  F, 

ODE  =  1,  ODK=l,  TDK=1,  MABL  =  2,  IMABL  =  0,  MABLK  =  T, 

RSI  =  2.084,  RWTD  =  0.05,  RZERO  =  -  I, 

IWALL=  2,  THETA  =  25,  RMAX  =  -0.0 1,  ZMAX  =  3.2, 

LWALL  =  1,  THALW  =  2.5,  XLW(I)  =  0,  YLW(I)  =  -1. 1 1,  XCOWL  =  O.I,  ICOWL  =  I, 
JET  =  T,  PINF  =  0.1,  MINF  =  16.4,  GINF  =  1.4, 

SEND 

C.2  TDK  Performance  Summary  for  the  DEAN 

1  TDK  PERFORMANCE  SUMMARY  :  DEAN_contourl00% 

**********FIRST  TDK  SOLUTION  :  WITH  RAMP  AND  COWL  BOUNDARY 

REAL  WALL  CONTOUR  1  ZONES  EXIT  PLANE 

FIRST 

TDK/MABL  SOLUTION 


THRUST  CHAMBER  OPERATING  CONDITIONS 


CHAMBER  PRESS  [PSIA] 
CHAMBER  TEMP  [R] 
MIXTURE  RATIO  [-] 

H  (OXID)  [CAL/MOLE] 

H  (FUEL)  [CAL/MOLE] 

HCHAM  (ODE)  [BTU/LB] 
DELH  (AVERAGE)  [BTU/LB] 
DELHI  (AVE)  [BTU/LB] 


1739.000 

6586.000 

7.000000 

0.000000 

0.000000 

-394.6774 

0.000000 

0.000000 


THRUST  CHAMBER  GEOMETRY 


ECRAT  [-]  4.000000 

RI  [-]  O.lOOOOOOE-01 

THETAI  [DEGREES]  30.00000 

RWTU  [-]  1.000000 

RSI  [INCHES]  2.084000 

RWTD  [-]  0.5000000E-01 

NIT  [-]  248.0000 

THE  [DEGREES]  13.06846 

THETA  [DEGREES]  25.00000 


EP  (NOZZLE)  [-] 


1.000010 
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ASURF  [IN*  *2] 


41.76751 


EXIT  FLOW  PROPERTIES 

P  (AXIS, EXIT)  [PSIA]  324.6112 

P  (WALL, EXIT)  [PSIA]  688.3699 

T(WALL,EXIT)  [R]  6363.542 

V(WALL,EXIT)  [FT/SEC]  14975.46 

MA  (WALL,EXIT)  [-]  2.963 189 

ONE-DIMENSIONAL  FLOW  PERFORMANCE 

ISP  (ODE)  [SECONDS]  509.8608 

ISP(ODK)  [SECONDS]  571.7090 

ISP  (ODF)  [SECONDS]  509.8608 

CSTAR  (ODE)  [FT/SEC]  0.0000 

CSTAR  (ODK)  [FT/SEC]  1462.3 

CSTAR  (ODF)  [FT/SEC]  0.0000 

TWO-DIMENSIONAL  FLOW  PERFORMANCE 

CD  [-]  0.000000 

CF(TDK)  [-]  2.831902 

CSTAR  (TDK)  [FT/SEC]  6302.546 

THRUST  (TDK)  [POUNDS]  67192.85 

WDOT(TDK)  [LB/SEC]  121.1250 

ISP  (TDK)  [SECONDS]  554.7396 

BOUNDARY  LAYER  PARAMETERS 

DFOPT(MABLK)  [POUNDS]  1176.185 

DF  (MABLK)  [POUNDS]  1025.826 

DISP  (MABLK)  [SECONDS]  8.469148 

THETA  (EXIT)  [INCH]  0.8708488E-03 

DEL*  (EXIT)  [INCH]  0.980 1783E-03 

DEL*  (THROAT)  [INCH]  0.7956905E-03 

EP  (REGEN)  [-]  1.222427 

SQDOT  (REGEN)  [BTU/SEC]  6711.148 

SQDOT  (LOSS)  [BTU/SEC]  0.000000 

SUMQDOT  [BTU/SEC]  6711.148 

DH(SUMQDOT)  [BTU/LBM]  55.40678 

THRUST  CHAMBER  PERFORMANCE 

THRUST  (TC)  [POUNDS]  66167.02 

CF  (TC)  [-]  2.790799 

WDOT(TC)  [LB/SEC]  121.1250 

ISP  (TC)  [SECONDS]  546.2704 

1  TDK  SCRAMJET  SUMMARY  :  DEAN_contourl00% 

**********FIRST  TDK  SOLUTION  :  WITH  RAMP  AND  COWL  BOUNDARY 

THRUST  ON  THROAT  PLANE  (FXBE)  =  61740.66  (LBF/FT) 

THRUST  ON  EXHAUST  NOZZLE  (FX  WALL)  =  4867.923  (LBF/FT) 
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THRUST  ON  COWL  (FX  COWL)  =  -441.5600  (LBF/FT) 

TOTAL  AXIAL  THRUST  (FX  TOTAL)=  66167.02  (LBF/FT) 
RADIAL  THRUST  ON  THROAT  (FY  BE)  =  0.000000  (LBF/FT) 

RADIAL  THRUST  ON  NOZZLE  (FY  WALL)  =  0.000000  (LBF/FT) 

RADIAL  THRUST  ON  COWL  (FY  COWL)  =  0.000000  (LBF/FT) 

TOTAL  RADIAL  THRUST  (FY  TOTAL)=  0.000000  (LBF/FT) 

COWL  LENGTH  (XCOWL)  =  0.1736667E-01  (FEET) 

COWL  LENGTH  (XCOWL)  =  0.1000000  (NONDIM) 


XBAR  WALL  (XB  WALL)  =  0.000000 

XBAR  WALL  (XB  WALL)  =  0.000000 

WALL  TORQUE  (FY  WALL)*(XBAR  WALL)  = 
YBAR  WALL  (YB  WALL)  =  0.000000 

YBAR  WALL  (YB  WALL)  =  0.000000 

WALL  TORQUE  (FX  WALL)*(YBAR  WALL)  = 


(FEET) 

(NONDIM) 

0.000000  (FT*LB/FT) 
(FEET) 

(NONDIM) 

0.000000  (FT*LB/FT) 


XBAR  BRNR  (XB  BRNR)  =  0.000000 

XBAR  BRNR  (XB  BRNR)  =  0.000000 

BRNR  TORQUE  (FY  BRNR)*(YBAR  BRNR)  = 
YBAR  BRNR  (YB  BRNR)  =  0.000000 

YBAR  BRNR  (YB  BRNR)  =  0.000000 

BRNR  TORQUE  (FX  BRNR)*(YBAR  BRNR)  = 


(FEET) 

(NONDIM) 

0.000000  (FT*LB/FT) 
(FEET) 

(NONDIM) 

0.000000  (FT*LB/FT) 


XBAR  COWL  (XB  COWL)  =  0.000000 

XBAR  COWL  (XB  COWL)  =  0.000000 

COWL  TORQUE  (FY  COWL)*(XBAR  COWL)  = 
YBAR  COWL  (YB  COWL)  =  0.000000 

YBAR  COWL  (YB  COWL)  =  0.000000 

COWL  TORQUE  (FX  COWL)*(YBAR  COWL)  = 


(FEET) 

(NONDIM) 

0.000000  (FT*LB/FT) 
(FEET) 

(NONDIM) 

0.000000  (FT*LB/FT) 


TOTAL  TORQUE  ABOUT  THE  EXIT  PLANE  =  0.000000  (FT*LB/FT) 


PERCENTAGES  OF  THRUST  COMPONENTS  (RELATIVE  TO  AXIAL  THRUST) 
THRUST  ON  THROAT  PLANE  (FX  BE)  =  93.3 1  % 

THRUST  ON  EXHAUST  NOZZLE  (FX  WALL)  =  7.36  % 

THRUST  ON  COWL  (FX  COWL)  =  -0.67  % 

RADIAL  THRUST  ON  THROAT  (FY  BE)  =  0.00% 

RADIAL  THRUST  ON  NOZZLE  (FY  WALL)  =  0.00% 

RADIAL  THRUST  ON  COWL  (FY  COWL)  =  0.00  % 

TOTAL  RADIAL  THRUST  (FY  TOTAL)=  0.00  % 

THRUST  VECTOR  ANGLE  =  0.00  DEG 
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Appendix  D  :  Isp  Calculation 


This  appendix  shows  a  simple  ealculation  to  demonstrate  the  impact  or  hp  hinted 
to  in  Chapter  1.  Equation  19  is  the  change  in  velocity  absent  drag  and  gravity. 

M 

M  =  (19) 

In  Equation  19,  Au  is  the  change  in  velocity,  go  is  acceleration  due  to  gravity,  Isp 
is  the  specific  impulse.  Mo  is  the  initial  mass,  and  Mb  is  the  burnout  mass  or  the  mass  at 
the  end  of  the  thrust  period.  Using  the  gross  mass  of  the  Centaur  IIA  19,073  kg  as  the 
initial  mass^',  the  empty  mass  of  2293  kg  as  the  burn  out  mass^^  and  the  stated  kp  of  449 
s^',  result  is  in  a  Au  of: 

Am  =  9.8 1  ^  *  4495  *  In  ^g235  — 

2293kg  s 

Know  keeping  the  same  Au  goal,  same  initial  mass,  but  increasing  the  Isp  by  one 
second  yields  a  bum  out  mass  of: 

9235— =  9.81-^*4505  *  In 2354kg 

5  5 

In  this  analysis  it  is  assumed  the  propellant  mass  is  the  initial  mass  minus  the  bum 
out  mass.  Eor  the  lower  fip  the  propellant  mass  is  16,780  kg  and  for  the  higher  Isp  the 
propellant  mass  in  16,719  kg.  This  less  fuel  in  required  to  perform  the  same  Au 
maneuver  with  an  increase  in  fip.  The  result  is  a  savings  of  61  kg  (134.5  Ibm).  This 
weight  savings  can  mean  less  fuel,  as  shown,  or  larger  satellites. 
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Appendix  E  :  Wall  Temperatures 


This  section  displays  the  results  for  different  materials  and  their  affect  on  the  wall 
temperature.  Table  19  shows  tested  materials  for  the  hydrogen  wall, 


Table  19  Hydrogen  Wall  Materials 


CVHl 

CVH2 

CVH3 

CVH4 

CVH5 

CVH6 

CVH7 

CVH8 

hn 

0.00437 

0.007436 

0.011791 

0.006182 

0.003611 

0.003411 

0.003332 

0.003317 

Ah  (in^) 

54.98 

87.96 

2.83 

120.64 

82.94 

67.86 

60.32 

60.32 

Th(R) 

6586.32 

6586.32 

6586.32 

6586.32 

6586.32 

6586.32 

6586.32 

6586.32 

he 

0.197015 

0.131318 

0.115463 

0.131837 

0.147505 

0.154929 

0.154406 

0.154432 

Ac  (in^) 

32.97 

65.94 

2.2 

90.43 

52.74 

37.68 

30.14 

30.14 

Tc(R) 

145.391 

270.699 

276.61 

415.765 

477.206 

524.952 

567.469 

609.732 

Q  (Btu/s) 

1439.05 

3745.58 

185.042 

4328.38 

1761.69 

1349.22 

1159.36 

1146.54 

t(in) 

0.02 

0.02 

0.02 

0.02 

0.02 

0.02 

0.02 

0.02 

Copper 

k 

(Btu/sRin) 

0.005377 

0.005377 

0.005377 

0.005377 

0.005377 

0.005377 

0.005377 

0.005377 

T^e(R) 

371.6508 

703.1089 

980.9731 

771.1716 

700.7632 

753.2799 

813.6587 

853.1702 

TwH  (  R  ) 

468.999 

861.4854 

1224.16 

904.613 

779.7623 

827.2277 

885.1436 

923.8645 

50%  Tmelt 

1222.2 

1222.2 

1222.2 

1222.2 

1222.2 

1222.2 

1222.2 

1222.2 

Silicon 

Carbide 

k 

(Btu/sRin) 

0.006571 

0.006571 

0.006571 

0.006571 

0.006571 

0.006571 

0.006571 

0.006571 

Twe  (  R  ) 

372.2814 

705.1291 

986.1016 

772.5984 

701.2952 

753.792 

814.1962 

701.2952 

TwH  (  R  ) 

451.948 

834.7393 

1185.118 

881.8026 

765.9455 

814.3086 

872.697 

765.9455 

50%  Tmelt 

2790 

2790 

2790 

2790 

2790 

2790 

2790 

2790 

Tungsten 

k 

(Btu/sRin) 

0.002333 

0.002333 

0.002333 

0.002333 

0.002333 

0.002333 

0.002333 

0.002333 

Twe  (  R  ) 

367.121 

688.5986 

944.1369 

760.923 

696.9423 

749.6009 

809.7983 

696.9423 

TwH  (  R  ) 

591.4695 

1053.593 

1504.586 

1068.452 

879.0035 

920.021 

974.5421 

879.0035 

50%  Tmelt 

3294 

3294 

3294 

3294 

3294 

3294 

3294 

3294 

Nickel 

k 

(Btu/sRin) 

0.001216 

0.001216 

0.001216 

0.001216 

0.001216 

0.001216 

0.001216 

0.001216 

Twe  (  R  ) 

359.7719 

665.0572 

884.3742 

744.2957 

690.7434 

743.6322 

803.5352 

690.7434 

TwH  (  R  ) 

790.1648 

1365.266 

1959.546 

1334.263 

1040.012 

1070.568 

1119.582 

1040.012 

50%  Tmelt 

1555.2 

1555.2 

1555.2 

1555.2 

1555.2 

1555.2 

1555.2 

1555.2 
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In  Table  19,  the  important  values  are  the  hot-side  wall  temperature  (Twh)  and  the 
50%  melting  temperature.  As  expeeted,  the  extreme  Twh  is  for  the  CVH3  element, 
whieh  represents  the  throat.  Table  19  indieates  that  all  but  the  Niekel  would  meet  the 
wall  temperature  eriteria.  Moreover,  Table  19  shows  how  the  temperature  of  the  wall  is 
affeeted  by  the  thermal  eonduetivity  of  the  material  (k).  The  oxygen  side  is  presented  in 
Table  20. 


Table  20  Oxygen  Wall  Materials 


CVOl 

CV02 

CV03 

CV04 

CV05 

CV06 

bn 

0.011217 

0.005725 

0.003248 

0.003003 

0.002891 

0.002833 

Ah  (in^) 

3.14 

165.88 

180.96 

108.96 

180.96 

180.96 

Th(R) 

6586.32 

6586.32 

6586.32 

6586.32 

6586.32 

6586.32 

he 

0.087317 

0.062727 

0.053795 

0.049671 

0.046118 

0.03672 

Ac  (in^) 

1.57 

90.43 

105.5 

105.5 

105.5 

105.5 

Tc(R) 

179.391 

303.911 

379.602 

452.208 

531.483 

616.862 

Q  (Btu/s) 

177.326 

5107.5 

3305.34 

3019.91 

2859.56 

2702.71 

t(in) 

0.02 

0.02 

0.02 

0.02 

0.02 

0.02 

Copper 

k 

(Btu/sRin) 

0.005377 

0.005377 

0.005377 

0.005377 

0.005377 

0.005377 

Twc(R) 

1446.081 

1188.386 

955.7072 

806.6633 

1113.574 

1308.098 

TwH  (  R  ) 

1656.12 

1302.903 

1023.642 

909.7455 

1172.346 

1363.647 

50%  Tmelt 

1222.2 

1222.2 

1222.2 

1222.2 

1222.2 

1222.2 

Silicon 

Carbide 

k 

(Btu/sRin) 

0.006571 

0.006571 

0.006571 

0.006571 

0.006571 

0.006571 

Twc(R) 

1453.879 

1191.369 

956.8652 

807.7637 

1114.61 

1309.278 

TwH  (  R  ) 

1625.768 

1285.086 

1012.461 

892.1228 

1162.707 

1354.737 

50%  Tmelt 

2790 

2790 

2790 

2790 

2790 

2790 

Tungsten 

k 

(Btu/sRin) 

0.002333 

0.002333 

0.002333 

0.002333 

0.002333 

0.002333 

Twe(R) 

1390.07 

1166.962 

947.3899 

798.7596 

1106.13 

1299.628 

TwH  (  R  ) 

1874.125 

1430.879 

1103.952 

1036.323 

1241.577 

1427.645 

50%  Tmelt 

3294 

3294 

3294 

3294 

3294 

3294 
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Table  20  cont. 


Nickel 

k 

(Btu/sRin) 

0.001216 

0.001216 

0.001216 

0.001216 

0.001216 

0.001216 

Twc(R) 

1299.199 

1132.204 

933.8958 

785.9367 

1094.054 

1285.886 

TwH  (  R  ) 

2227.816 

1638.505 

1234.246 

1241.681 

1353.897 

1531.476 

50% 

Tmelt 

1555.2 

1555.2 

1555.2 

1555.2 

1555.2 

1555.2 

Cobalt 

k 

(Btu/sRin) 

0.00133 

0.00133 

0.00133 

0.00133 

0.00133 

0.00133 

Twc(R) 

1299.199 

1132.204 

933.8958 

785.9367 

1094.054 

1285.886 

TwH  (  R  ) 

2148.247 

1595.122 

1208.51 

1202.63 

1331.632 

1510.432 

50%  Tmelt 

1592.1 

1592.1 

1592.1 

1592.1 

1592.1 

1592.1 

Table  20  shows  how  Twh  at  the  throat,  CVOl,  exeited  the  50%  melting  point  for 
not  only  eopper,  but  niekel  and  eobalt  as  well.  Of  the  materials  shown  only  silieon 
earbide  and  tungsten  meet  the  requirements  due  to  their  high  melting  points. 
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